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ABSTRACT 


A machine computational program has been developed under NASA Contract 
MS 3- 13317 to calculate temperature and reactant concentrations as functions of 
time and axial position in typical catalyzed hydrogen-oxygen reaction chambers. The 
program is based upon a transient model of the reactor system which describes the 
behavior of reactors operated under conditions of continuous flow. 

The computer program developed from this model is described in detail in this 
computer manual. The manual contains operating instructions for this program as 
well as descriptions of input and output formats, including all output messages. 
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FOREWORD 


This work was performed by United Aircraft Research Laboratories for the 
National Aeronautics and Space Administration under Contract NAS 3-13317 initiated 
January l4, 1970- 

Included among those who cooperated in performance of the work under NAS 3-13317 
were Dr. A. S. Kesten, Program Manager, Dr. W. G. Burwell, Chief, Kinetics and 
Thermal Sciences Section, Mr. D. B. Smith, Project Analyst, and Mrs. E. J. Smith, 
Computer Analyst. 

This work was conducted under program management of the NASA Lewis Research 
Center and the Technical Manager was Mr. P. N. Herr, NASA Lewis Research Center, 
Cleveland, Ohio. 
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Transient Model of Hydrogen/Oxygen Reactor 
Contract MS 3-13317. 


SUMMARY 


The Research Laboratories of United Aircraft Corporation, under Contract 
MAS 3-13317 "with the Rational Aeronautics and Space Administration, have performed 
an analytical study of the transient behavior of a hydrogen/oxygen catalytic igni- 
tion system. This study has included the development of a computer program which 
is used to calculate temperature and species concentration profiles as functions of 
time in typical reaction chamber configurations. This manual describes the com- 
puter program. A description of input and output for the transient program is 
included in the discussion together with examples for typical data cases. In addi- 
tion, a short write-up of the subroutines contained in the deck is included. 
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INTRODUCTION 


Under Contract MS 3-13317? a comprehensive analytical program was formulated 
with the objectives of (l) .developing a transient analysis of a continuous flow 
hydrogen-oxygen catalytic ignition system to permit the prediction of temperatures? 
pressures? and species concentrations as functions of time as well as axial position 
in typical reactor systems? (2) developing a computer program based on this analysis? 
and (3) performing calculations using this computer program to demonstrate the 
effects of various system parameters on the transient performance of the reactor? 
and thus to define ignition delay times. As part of the contract effort? attention 
•has been directed toward preparing a manual describing to potential users the 
operation of the computer program. The manual includes a detailed discussion of the 
program. 
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DISCUSSION OF THE H 2 /0 2 TRANSIENT COMPUTER PROGRAM 


The equations representing "the transient model of a H 2 /O 2 catalytic reactor 
have been programmed for the UNIVAC 1108 digital computer using FORTRAN V.* This 
computer program is discussed below. Included in this discussion are input and 
output descriptions and a list and short description of the various routines, in the 
program. 


Input Description 

The following is a description of the necessary input for the Hg/C^ transient 
computer program. To complete each data case run with this program, a minimum of 
two separate runs must be made; a third type of run may also be necessary. The 
three types of runs are: (l) initial non-diffusion- controlled run (always necessary) 

(2) continuation of a non- diffusion- controlled run (sometimes needed), and (3) 
diffusion- controlled (concluding) run (always necessary). The non-diffusion- 
controlled runs (types 1 and 2) are based upon a very general reasonably complete 
model of the reactor system. It would be best to run this type of run all the time, 
but the large amount of computer time needed necessitates restricting its use to 
the lower temperature calculations and then switching to the much faster diffusion- 
controlled run which is accurate at the higher temperatures. Splitting each case 
up into two or more runs this way also allows for the use of fewer axial stations 
in the non-diffusion-controlled runs, which reduces computer time for these slower 
runs . 


Initial Non-Diffusion-Controlled Run 


This type of run is the first one made when running a transient Hg/ 0 2 case. A 
general input format is given in Table I. The coding of a sample data case is shorn 
in Fig. 1, and a listing of the input data punch cards corresponding to this sample 
data case is shown in Fig. 2. Below is a detailed description of the input; the 
numbers correspond to the card numbers (first column) of Table I. 

1. The first card contains the title of the data case plus two logical indica- 
tors, H2LEAD and DEBUG. 

Title: The title serves to give individual data case identification. 

It consists of any alpha numeric information desired and should 
be punched in card columns I- 78 . 


H2EEAD: This is a logical unit which should indicate whether the catalyst 

in the reactor is preconditioned with hydrogen (hydrogen run 
through the reactor at its inlet concentration before the oxygen 
is turned on) . If it is preconditioned, the letter "T" should 

^Discussion of these equations and their derivation can be found in NASA Contract 
Report NASA CR-120799, "Final Report - Transient Model of Hydrogen /Oxygen Reactor", 
February, 1971. 
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appear in card column 79- If it is not, the letter "F" should be 
punched. 

DEBUG: This logical -unit indicates whether the detailed printing is 

desired. . A large amount of print will be generated under this 
option. For example, a data case where N0FZ=l6 and NPASS=5 (see 
description of these variables on card 2) will produce over 150 
pages of detailed print. This printing will sometimes be neces- 
sary, however, since s umma ry printing only gives the temperatures 
and concentrations at the end of every twenty time steps. To get 
the detailed print, the letter "T" should be punched in card column 
80 ; if the letter "F" is punched, only the summary print will 
result . 

2. The second card contains five variables, NOFZ, NPASS , INFERP, MODNO , PUNCH, 

and ITSTEP. 

NOFZ: This is the number of axial positions (Z's) at which analysis will 

be made. These axial positions will be input on cards 6a, 6b, etc. 
The maximum number of Z's allowed by the program is 75- However, 
computer run time for the non-diffusion-controlled runs is 
affected drastically by the number of Z 's. For these runs, it is 
best to use as few as possible (10-15 were used in test cases), 
especially for lower pressures. (Using 10 axial positions, one 
data case run at P=15 psia took over 30 minutes of computer time 
for the non-diffusion-controlled run. A comparable case run at 
P=300 psia took only four minutes.) 

NPASS: This is the total number of time "passes" to be taken -for this run. 

Each "pass" consists of "INEERP" time steps (see below), and during 
each pass, a constant time step size is used. For the non- 
diffusion-controlled runs, it is best to use NPASS=10 or less 
although the maximum number allowed by the program is 40. If the 
case completes the non-diffusion-controlled run before the 10 passes 
are through, the program will take care of itself (i.e., punch out 
the appropriate cards for the diffusion-controlled run (run type 3) 
and end) . If more than 10 passes are needed (as in the low pres- 
sure cases) the program will punch out the cards for run type 2 
(continuation of non-diffusion- controlled run) and end. For the 
two cases mentioned under the description of NOFZ, the low pressure 
case which took over 30 minutes to complete needed over 17 passes, 
whereas the high pressure case needed fewer than 2 passes which 
were completed in under 4 minutes. 

If it is decided to input the time increments (see ITSTEP) NPASS 
must equal the number of input DELT's. However for non-diffusion- 
controlled runs, inputting the DELT's is not advised. 
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INEERP: This is the number of time increments per pass. The number used 

in the test cases was "20" which seemed to work out fairly well. 
Since summary plotting and printing is done only at the end of 
each pass, at times it might be desirable to use a smaller number 
for this .variable. 

MODUO : This number indicates which axial position values will be plotted 

on the summary printer-plots. If MODNO=1, every point will be 
plotted, if it equals 2, every-other point will be plotted, etc. 

PUNCH: This variable indicates which of various read/write punch card 

options will be used. For an initial non-diffusion- controlled 
run, this variable should = 1 . This indicates cards will be 
punched at the end of the run (either when NPASS is reached or 
when the non-diffusion-controlled part of the run is completed) . 

ITSTEP: This value indicates whether time steps will be input or calculated 

by the program. For non-diffusion-controlled runs, ITSTEP should 
equal 0 which indicates the program will calculate the time steps. 

3. The third card contains the 8 values for G-, TS, T(l), P(l), H(l), CIH2(l), 

CI02(1) and CIHE(l). 

G: This is the inlet (at Z=0) mass flow rate in Ib/ft^-sec. It must 

be greater than zero. 

TS: This is the freezing temperature of water in deg R. 

T(l): This is the inlet interstitial temperature in deg R. 

P(l): This is the steady-state inlet chamber pressure in psia. 

H(l): This is the inlet enthalpy in BTU/lb, (reference number, usually 

taken as 0) . 

CIH2(l): This is the inlet concentration of hydrogen in Ib/ft^. 

CI02(l): This is the inlet concentration of oxygen in lb/ft3. 

CIHE(l): This is the inlet concentration of helium in lb/ft^. 

4. The fourth card contains the 8 values, TI, PI, CH2I, C02I, CH20I, CHEI, PTI 

and SDTIME. 

TI: This is the initial (at time = 0) interstitial temperature through- 

out the reactor in deg R. 
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PI: This is the initial chamber pressure throughout the reactor in 

psia. 

CH2I: This is the initial concentration of hydrogen throughout the 

reactor in lb/ft3. 

C02I: This is the initial concentration of oxygen throughout the reactor 

in lb/ft 3 . 

CH20I: This is the initial concentration of -water throughout the reactor 

in lb/ft 3 . 

CHEI: This is the initial concentration of helium throughout the reactor 

in lb/ft^. 

PTI: This is the initial catalyst bed temperature throughout the reactor 

in deg R. 

SDTIME: This is the time in seconds at which the oxygen flow is shut off. 

5. The fifth card contains the 8 values AW, CW, MW, DC, TA, HA, HA1 and HA2. 

AW: This is the total surface area of the chamber plus nozzle walls in 

ft 2 . 

CW: This is the specific heat of the chamber walls in BTU /lb -deg R. 

MW: This is the thermal mass of the chamber plus nozzle walls in lbs. 

DC: This is the diameter of the reaction chamber in ft. 

TA: This is the temperature of the surrounding atmosphere in deg R. 

HA: This is the heat transfer coefficient for forced convection between 

the chamber and the surrounding- atmosphere in BTU/ft^-sec-deg R. 

HA1: This is the heat transfer coefficient for natural convection between 

the chamber and the surrounding atmosphere in BTU/ft 2 -sec-deg . 

HA2: This is the radiative heat transfer coefficient between the chamber 

and the surrounding atmosphere in BTU/ft^-sec-deg R . 

6. Card six should only be included for those cases where the catalyst is not 

preconditioned with hydrogen (H2LEAD = F on card l) . This card contains the 

6 values HP, TSS, MBARSS, AC, VC and PREST. 
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PF: This is the feed pressure in psia. 

TSS: This is steady-state hed exit temperature in deg R. This number 

can he obtained by running the H 2 /O 2 steady- state computer program 
(see UARL Report H910721 . . . MSA CR- 72567 ). 

MBARSS: This is the steady-state bed exit average molecular weight. It can 

also be obtained from running the steady-state computer program. 

O 

AC: This is the cross-sectional area of the reaction chamber in ft . 

VC: This is the volume of the reactor up to the nozzle throat exclu- 

sive of the volume occupied by catalyst particles in ft3, 

PREST: This is the pressure of the surrounding atmosphere in psia. 

7, Cards 7a, 7b, etc. contain the axial station values (Z's) in feet at which 
calculations will be made. There should be NOFZ (see card 2) values of Z 
input at the rate of 8 per card. For most cases, it is recommended that 
more Z’s be placed near the reactor inlet. The first Z must equal 0; the 
second Z must be very close to 0; the last Z equals the bed exit. 

8. Cards 8a, 8b, etc. contain the input time increments (BELT'S ) in seconds. 
These cards should be included only if ITSTEP on card 2 equals 1. For 
non- diffusion- controlled runs, it is advisable to allow the program to cal- 
culate DEBT’S (lTSTEP=0) . 


Cards 9 through 11 contain AVSZ(l), the interpolation table used to obtain the 
catalyst particle radius at any point along the reactor bed. Subroutine 
UNBAR, an interpolation routine developed at UARL, is used to obtain an appro- 
priate particle radius. A, for a given axial station along the bed. For this 
table, there should be a total of NZ (see third entry on card 9)values of ZandA 

9. Card 9 contains the four table descriptors needed by UNBAR. The first 
descriptor signifies the table number. For this program, it should equal 
0.0. The second descriptor tells at what location in the array the table 
starts; the tables in this program are read in such that this number equals 
1.0. The third descriptor is the number of independent variables in the 
table (in this case the number of Z's). This number is called NZ. The 
fourth descriptor for a univariate table such as this one should equal 0.0. 

10. Cards 10a, 10b, etc. contain the monatonically increasing Z values in feet. 
Enough cards should be used to contain NZ values of Z at the rate of 10 per 
card. For example, if NZ=12, 12 values of Z should be input using two cards, 
with ten values on the first card and the two remaining values on the second 
card. 
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11. Cards 11a, lib, etc. contain the A's (in feet) which correspond to the Z's 
listed on cards 10a, 10b, etc. Enough cards should be used to contain HZ 
values of A at the rate of 10 per card. 

Cards 12 through 14 contain AFVSZ(l), the interpolation table used to obtain 
the total external catalyst particle surface area per unit volume of bed 
(AP) . These AP values are obtained as functions of Z using UNBAR as in the 
AVSZ table discussed above. For this table there should be a total of NZ 
values of Z and AP. 

12. Card 12 contains the table descriptors and should be exactly the same as 
card 9. 

13. Cards 13a, 13b, etc. contain the Z values and should be exactly the same 
as cards 10a, 10b, etc. 

14. Cards l4a, l4b, etc. contain the AP values in ft -1 which correspond to the 
Z's listed in cards 13a, 13b, etc. Enough cards should be used to contain 
HZ values of AP at the rate of ten per card. 

Cards 15 through 17 contain DELVSZ (i ) , the interpolation table used to 
obtain the interparticle void fraction (DELTA). These DELTA values are 
obtained as functions of Z using UNBAR as in the AVSZ table discussed 
above. For this table there should be a total of IE values of Z and DELTA. 

15. Card 15 contains the table descriptors and should be exactly the same as 
cards 9* 

16. Cards l6a, l6b, etc. contain the Z values and should be exactly the same as 
cards 10a, 10b, etc. 

17. Cards 17a , 17b, etc. contain the DELTA values which correspond to the Z's 
listed on cards l6a, l6b, etc. Enough cards should be used to contain NZ 
values of DELTA at the rate of 10 per card. 

18. Card 18 contains the values CIH2(l) and CI02(l). These two values are the 
inlet concentrations of hydrogen and oxygen, respectively, after "shutdown” 
(that is, after the oxygen flow is turned off). 

Continuation of a Non-Diffusion-Controlled Run 


If the initial non-diffusion-controlled run completes the input NPASS number of 
passes without printing out "GO TO ALL DIFFUSION-CONTROLLED RUN", a continuation run 
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will be- necessary. Also if the continuation run ends without this printed message, 
a continuation of the continuation run should be made, etc, A general input format 
for this type of run is given in Table II. The listing of the punch cards for a 
continuation run is given in Tigs. 3a - 3e- Below is a detailed description of the 
input; the numbers correspond to the card numbers (first column) of Table II. 

1. The first card contains, again, the title H2LEAD and DEBUG. This card 
should be exactly the same for the continuation runs as it was for the 
initial non-diffusion-controlled run. 

2. The second card again contains HOFZ , EPASS, H'JIERP, MODUO, POUCH and ITSTEP. 

HOFZ: This value must be the same as it was for the initial run. 

HPASS: This is the number of time "passes" to be taken. It may be 

changed, but the maximum still equals 40. 

INFERP: This is the number of time increments per "pass". It may be 

changed. 

MODUO: This number indicates which axial position values will be plotted 

as before. It may be changed. 

PUUCH: This indicator specifies which read/write punch card option will 

be used. For continuation runs, it must equal 3 » This indicates 
punch cards mil be read in with the data case and also punched 
out at the end of the run. 

ITSTEP: As in the initial run, this indicates whether time increments will 

be input or internally calculated. Again, as in the initial run, 
it should equal 0 (internally calculated) . 

3-17- Cards 3 through 17 should be exactly the same as they were for the initial 
run. 

Cards 18- 32 are obtained from the cards punched out at the end of' an initial 
non-dlffusion-controlled run or at the end of a previous continuation run 
if this- is the 2nd or more continuation. These cards should only be included 
for continuation runs . 

18. Cards l8a, 18b, etc. contain EOFZ values of Z(l) . . . (axial positions). 

19. Cards 19a, 19b, etc. contain HOFZ values of PT(l) . . . . (particle surface 
temperature at each axial station at end of previous run) . 
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20. Cards 20a, 20b, etc. contain WOFZ values of T(l) . . . (interstitial 
temperature at each axial station at end of previous run) . 

21. Cards 21a, 21b, etc. contain WOFZ values of H(l) . . . (enthalpy at each 
axial station at end of previous run) . 

22. Cards 22a, 22b, etc. contain WOFZ values of P(l) . . . (chamber pressure 
at each axial station at end of previous run) . 

23. Cards 23a, 23b, etc. contain WOFZ values of CIH2(l) . . . (interstitial 
concentration of hydrogen at each axial station at end of previous run). 

24. Cards 24a, 24b, etc. contain WOKS values of CI02(l) . . . (interstitial 
concentration of oxygen at each axial station at end of previous run) . 

25. Cards 25a, 25b, etc. contain WOFZ values of CIH20V(l) . . . (interstitial 
concentration of water at each axial station at end of previous run) . 

26. Cards 26a, 26b, etc. contain WOFZ values of CIHE(l) . . . (interstitial 
concentration of helium at each axial station at end of previous run) . 

27. Cards 27a, 27b, etc. contain WOFZ values of PPT(l) . . . (particle surface 
temperature at each axial station one time increment previous to end of 
previous run) . 

28. Cards 28a, 28b, etc. contain WOFZ x WPP values of TPP(l, J). . . 
(temperature, profile within catalyst particle at each axial station at end 
of previous run) . WPP is the number of radial positions within the particle 
considered in the analysis (set at 30 in the program). 

29. Cards 29a, 29b, etc. contain WOFZ x WPP values of C02(l, j). . . (oxygen 
concentration profile within catalyst particle at each axial station at end 
of previous run) . 

30. Cards 30a, 30b, etc. contain WOFZ x WPP values of CH2(l, J) . . . (hydrogen 
concentration profile with catalyst particle at each axial station at end of 
previous run) . 

31. Cards 31a, 31b, etc. contain WOFZ values of MHJT(l) . . . (diffusion- 
controlled indicator for each axial station at end of previous run) . 

32. Card 32 contains the values at the end of the previous run for TIME,' TW, and 
DELTT. 
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TIME: is the cumulative reactor run time in seconds 

TW: is the wall temperature in deg R 

DELTT: is the time increment 

33. Card 33 for the continuation run is exactly the same as card 18 of the 
initial non- diffusion- controlled run. 

Diffusion-Controlled Run 

This type of run is the last one to be made when running a transient H2/O2 ease. 
When "GO TO ALL DIFFUSION- CONTROLLED RUN" is printed out at the end of an initial 
or continued non-diffusion-controlled run, this third type of run should follow. 
Compared to the first two types, this kind of run takes relatively little computer 
time. Depending on the number of axial stations being analyzed, a diffusion- con- 
trolled run, on the UARL UNIVAC 1108, took between three and five minutes for all 
cases run (between 20 and 40 Z’s were used). 

A general input format is given in Table III. A listing of the input data punch 
cards for a run of this type is shorn in Fig. 4. Below is a detailed description of 
the input; the numbers correspond to the card numbers (first column) of Table III. 

1. The first card again contains the title of the data case plus the logical 
indicators H2LEAD and DEBUG. 

Title: This should be the same as it was for the two previous run types. 

H2LEAD: This also should be the same as it was for the non- diffusion- 

controlled runs. 

DEBUG: This indicator should be set equal to "F" for the diffusion 

controlled ruins; o the raise over 1000 pages of output may be pro- 
duced. For the diffusion- controlled runs, the summary printing 
gives sufficient information. 

2. The second card contains the five values of NOFZ, NPASS, INEERP, M0DN0, 

PUNCH and ITSTEP. 

NOFZ: This is the number of axial positions (Z's) at which analysis will 

be made in the diffusion-controlled run. The maximum number 
allowed is 75. Since this type of run takes relatively little 
run time no matter how many Z's are used, it is advised to use the 
near the maximum (the test cases were rum using 36 Z's). 
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NPASS: This is the total number of time "passes” to be taken for this 

run. For the dif fusion-controlled runs, it is advised to input 
the time increments (see discussion of ITSTEP). If this is the 
case, NPASS must equal the number of time increments input on 
cards 8aj 8b, etc. ; the maximum number allowed by the program is 

40. 

INPERP: This is the number of time increments per pass. It may change 

from previous run value. 

MODNO: This number indicates which axial position values are to be 

summary plotted. It may be changed from previous run value. 

PUNCH: This variable indicates which of various read/write punch card 

options will be used. For diffusion-controlled runs, this num- 
ber must equal 4 . This indicates a special set of data cards will 
be read in (see discussion of cards 18-25). 

ITSTEP: This value indicates whether time steps will be input or calculated 

by the program (ITSTEP=1 means input . . . ITSTEP=0 means calculated). 
For diffusion-controlled runs, ITSTEP should equal 1 and the time 
increments (DELT's) should be input on cards 8a, 8b, etc. As the 
case approaches steady-state, larger DELT's might be desired (see 
sample input DELT's in Fig. 3). If ITSTEP is left equal to 0, all 
DELT's will be calculated equal to .001 sec for the diffusion- 
controlled runs. 

3-6. Cards 3-6 should be exactly the same as they were for the non- diffusion- 
controlled runs . 

7. Cards 7 a » 7b, etc., contain the axial station values (Z's) in feet at which 
calculations will be made by the program. There should be NOFZ (see card 
2) values of Z input at the rate of 8 per card. Note that this will be a 
different set of Z's than those used in the non-diffusion-controlled 

runs . . . there should be many more. However, as -before, for most cases 
it is best to place more Z's near the reactor inlet. Also, the first Z 
must still equal 0, the second Z must be very close to 0, and the last Z 
must be at the bed exit. 

8. Cards 8a, 8b, etc., contain the input time increments (DELT's) in seconds. 

These cards must be included if ITSTEP=1. It is recommended that these be 
input for diffusion-controlled runs. 

9-17* Cards 9-17 should be exactly the same as they were for non-diffusion- 
controlled runs . 
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Cards 18-25 are obtained from the cards punched out at the end of the non- 
diffusion- controlled runs (when the message "GO TO AIL DIFFUSION-COI'JTHOT .LED 
MM" is printed) . These cards should only be included for diffusion-controlled 
runs. 


18. Card 18 contains the number for NAXT. This value is the number of Z's used 
in the preceding non- diffusion- controlled runs. 

19. Cards 19a, 19b, etc., contain ZT(l) . . . (the axial stations used in the 
non-diffusion-controlled runs). Enough cards are used to punch out NAXT 
values of ZT at the rate of 8 per card. 

20. Cards 20a, 20b, etc., contain PTT(l) . . . (the particle surface temperature 
at each axial station at the end of the non-diffusion-controlled run). 

NAXT values of PTT are punched out. 

21. Cards 21a, 21b, etc., contain C02T(l) . . . (the oxygen concentration at 
each Z at the end of the non-diffusion-controlled run) . NAXT values of 
C02T are punched. 

22. Cards 22a, 22b, etc., contain CH2T(l) . . . (the hydrogen concentration at 
each Z at the end of the non-diffusion-controlled run) . NAXT values of 
CH2T arc punched. 

23. Cards 23a, 23b, etc., contain CH20T(l) . . . (the water concentration at 
each Z at the end of the non-diffusion-controlled run). NAXT values of 
CH20T are punched. 

24. Cards 24a, 24b, etc., contain TTl(l) . . . (the interstitial temperature at 
each Z from the end of the non- diffusion- controlled run). NAXT values of 
TTI are punched. 

25. Card 25 contains TIME and TW, the cumulative reactor run time in seconds 

and the wall temperature in deg E from the end of the non-diffusion-controlled 
run. 


26. Card 26' of the diffusion- controlled run is exactly the same as card 18 of the 
initial non- diffusion- controlled run. 


13 



TABLE I 

INPUT FORMAT 

TRANSIENT H 2 /0 2 COMPUTER PROGRAM: INITIAL NON-DIFFUSION-CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN Y 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

1 

1 

13A6, 2L1 

1-78 

Title 






T9 

H2LEAD 

- - 

Logical indicator 




80 

DEBUG 

- — 

Logical indicator 

2 

1 

6i4* 

1-4 

NOFZ 


Number of axial stations 




5-8 

NPASS 

- - 

Number of "passes" 




9-12 

INPERP 

- - 

Number of increments per "pass" 




13-16 

MODNO 

- - 

Axial position plotting indicator 




17-20 

PUNCH 

- - 

Punch card indicator 




21-24 

ITSTEP 

- “ 

Input time step indicator 

3 

1 

8EIO.5 

1-10 

G 

G 

Mass flow rate 




11-20 

TS 

- _ 

Freezing temperature of water 




21-30 

T(l) 

T± 

Inlet interstitial temperature 




31-4o 

P(l) 

P 

Steady-state inlet chamber pressure 




41-50 

H(l) 

h 

Inlet enthalpy 




51-60 

CIH2(1) 

Ci H 2 

Inlet concentration of hydrogen 




61-70 

cioe(i) 

c i°2 

Inlet concentration of oxygen 


- 


71-80 

CIHE(l) 

Ci HE 


4 

I 

8E10 . 5 

1-10 

TI 

*i 

Initial interstitial temperature 




11-20 

PI 

P 

Initial chamber pressure 




21-30 

CH2I 

C/S 

Initial concentration of hydrogen 




31-40 

C02I 

C ± °2 

Initial concentration of oxygen 


*right adjusted 









TABLE I (Cant'd) 

INPUT FORMAT 


TRANSIENT H 2 /0 2 COMPUTER PROGRAM: INITIAL NON -DIFFUSION- CONTROLLED RUN 


CARD 

NO. 

NO, OF 
CARDS 

FORTRAN T 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

4 

1 

8eio . 5 

41-50 

CH20I 

C, H 20 

Initial concentration of water 

(Cont’< 

0 


51-60 

CHEI 

. Ci He 

Initial concentration of helium 




61-70 

PTI 

(T P ) S 

Initial catalyst Led temperature 




71-80 

SDTIME 

- - 

Shutdown time 

5 

1 

8E10.5 

l- 

■10 

AW 

A W 

Chamber wall total surface area 




11- 

•20 

CW 

c w 

Chamber wall specific heat 




21- 

•30 

MW 

m w 

Chamber wall thermal mass 




31-40 

DC 

dc 

Reaction chamber diameter 




4l- 

-50 ■ 

TA 

Z'a 

Ambient temperature 




51-60 

HA 

h a 

Forced convection heat transfer coefficient 




61-70 

HA1 

K 

Natural convection heat transfer coefficient 




71-80 

HA2 

V, 11 

h a 

Radiative heat transfer coefficient 

6* 

1 

6E10.5 

l- 

-10 

PF 

P F 

Feed pressure 




11- 

-20 

TSS 

Ti ss 

Steady-state exit interstitial temperature 




21-30 

MBARSS 

M ss 

Steady-state exit average molecular weight 




31-40 

AC 

A C 

Cross-sectional area of reaction chamber 




41-50 

VC 

v c 

Volume of reactor minus catalyst volume 


• 


51-60 

PREST 


Ambient pressure 

7a 

** 

8E10.5 

1- 

■10 

z(i) 

z 

Axial station 

7h 



11- 

-20 




etc. 











71- 

■80 





*This card is included only if H2LEAD (Card l) is equal to "F" 

**Enough cards should he used to contain (NOFZ) values of z at the rate of 8 per card 




TABLE I (Cant'd) 

INPUT FORMAT 

TRANSIENT H 2 /02 COMPUTER PROGRAM: INITIAL NON-DIFFUSION -CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

fortran y 

FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

8 a* 

** 

8E10.5 

1 - 

-10 

BELT (I ) 

At 

Time increments 

8 h 



11. 

-20 




etc. 











71-80 




9 

1 

4E8.it 

1-8 

0. 


Table descriptor 




9-16 

1 . 

- - 

Table descriptor 




17-24 

NZ. 

- - 

Table descriptor 




25 . 

-32 

0. 

» — 

Table descriptor 

10 a 

*** 

10E8.4 

1-8 

Z(I) 

z 

Axial stations 

lOh 



9- 

•16 




etc. 











73- 

-80 




11 a 

V- V v 
A A A 

10E8.4 

1-8 

A(I) 

a 

Catalyst particle radii 

lib 



9-16 




etc. 











73-80 




12 

1 

4e8.4 

1-8 

0. 

M M 

Table descriptors 




9-16 

1. 

- - 

• 




.17-24 

NZ. 

- _ 





25-32 

0. 

- - 



*These cards are included only if ITSTEP (Card 2) is equal to "1” 

**Enough cards should he used to contain (NPASS) values of At at the rate of 8 per card 
***Enough cards should he used to contain (NZ) values of z, a, A p , or 8 at the rate of 10 per card 




% 


TABLE I (Cant’d) 

INPUT FORMAT 


TRANSIENT H 2 /o 2 COMPUTER PROGRAM: INITIAL NON-DIEFUSION-CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

fortran y 

FORMAT 

COLUMNS 

USED 

symbol or 

DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

13a 

* 

10E8.4 

1-8 

Z(I) 

z 

Axial stations 

13b 



9- 

-16 




etc. 











73- 

-80 




l4a 

* 

10E8.4 

l- 

-8 

AP(l ) 

Ap 

Total external catalyst particle surface areas 

l4b 



9- 

-l 6 



per unit volume of bed 

etc. 







• 




rn 

C- 

•80 



. 

15 

1 

4e8.4 

1-8 

0 . 


Table descriptors 




9-16 

1 . ■ 

- - 

” 




17-24 

NZ. 

_ - 





25- 

■32 

0 . 

- ~ 


l 6 a 

• * 

10E8.4 

1 - 8 ’ 

Z(I). 

z 

Axial stations 

l 6 h 



9-16 




etc . 









' 


73-80 




■ 

IT a 

*■ 

10E8.4 

1-8 

DELTA (I ) 

8 . 

Interparticle void fractions 

17 b 



9- 

l 6 




etc. 











73- 

80 





^Enough cards should he used to contain (HZ) values of z, a , A p; ' or 8 at the rate of 10 per card 








































. TABLE I (Gont'd) 

INPUT FORMAT 

TRANSIENT H 2 /0 2 COMPUTER PROGRAM: INITIAL NON-DIFFUSION-CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN Y 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

l8 

I 

2E10.5 

1-10 

11-20 


0^2 

Inlet concentration of hydrogen after "shutdown" 
Inlet concentration of oxygen after "shutdovna" 









TABLE II 

INPUT FORMAT 

TRANSIENT H 2 /0 2 COMPUTER PROGRAM: CONTINUATION OF NON-DIFFUSION-CONTROLIFD RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN T 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

1-17 

Cards 1-17 follow exact lj 
in values on card 2 in t! 

the same fo 
is detailed d 

rmat as cards 1-: 
iscussion of coni 

L7 in Table I. Note, however, the differences 
;inuat ion-run input. 

Care 

Ls I 8 - 3 : 

. are cards 

punched out by the program at the end of an initial or previous continuation run. 

18 a 

l8b 

etc. 

* 

8 EIO. 5 . 

1-10 

11-20 

I 

71-80 

Z(I) 

i Z 

Axial stations 

19 a 

19 b 

etc. 

* 

8E10 . 5 

If 

PT(l) 

(T ) 
v p's 

Particle surface temperatures 

20a 

20b 

etc. 

* 

8 EIO .5 

tf 

t(t) 

T. 

i 

Interstitial temperatures 

21a 

21b 

etc. 

* 

8E10 , 5 

If 

h(i) 

h 

Enthalpies 


■^Enough cards should "be used to contain (NOFZ) values at the rate of 8 per card. 




TABLE II (Cant'd) 

INPUT FORMAT 

TRANSIENT H 2 /0 2 COMPUTER PROGRAM : CONTINUATION OF NON -DIFFUSION-CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN T 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

22 a 
22 t 
etc . 

* 

8E10.5 

1-10 

11-20 

71-80 

P(I) 

P 

Chamber pressures 

23a 

23 b 

etc. 

* 

8E10.5 

tt 

CIH2(I) 


Interstitial hydrogen concentrations 

24a 
24b 
etc . 

* 

8E10.5 

t? 

CI02(I) 

=i° 2 

Interstitial oxygen concentrations 

25a 
25 b 
etc . 

* 

8E10.5 

IT 

CIH20V ( I ) 

Ci H 2 ° V 

Interstitial water concentration 

26 a 

26 b 

etc. 

'* 

8E10.5 

ft 

cihe(i) 

=i He 

Interstitial helium concentration 

27 a 

27 b 

etc. 

* 

8E10.5 

tt 

ppt(i) 

<V« 

Particle surface temperatures at previous time 


*Enough cards should "be used to contain (NOFZ) values at the rate of 8 per card. 



















TABLE II (Cant’d) 

INPUT FORMAT 

TRANSIENT H 2 /c> 2 COMPUTER PROGRAM: CONTINUATION OF N ON -DIFFUSION-CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN Y 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

28 a 

** 

8E10.5 

1- 

■10 

TPP(I,J) 

T(x,t) 

Particle temperatures 

28 b 



11. 

-20 




etc . 



1— 1 

0 

°o 




29 a 

** 

8E10.5 

IT 

C02(I,J) 

-P 

*\ 

of 

Particle oxygen concentration 

29 b 








etc. 








30 a 

30 b 

** 

8E10.5 

It 

CH2(I,J) 


Particle hydrogen concentration 

etc. 








31 a 

* 

4012** *** 

1- 

-2 

MINT (I) 


Diffusion-controlled indicators 

31b 



3rb 




etc. 











79-80 




32 

1 

3E10.5 

1- 

-10 

TIME 

t 

Cumulative time 




11- 

-20 

TW 

T w 

Wall temperature 




21- 

•30 

DELTT 

At 

Time increment 

33 

1 

2E10.5 

1- 

•10 

CIE2(1) 

Ci H 2 

Inlet concentration of hydrogen after "shutdown” 




11- 

-20 

CI02(l) 

Ci°2 

Inlet concentration of oxygen after "shutdown" 


**Enough cards should he used to contain (NOFZ x NPP) values at the rate of 8 per card (NPP values for first Z, then NPP 
values for second Z . . . etc. . . to NPP values for (NOFZ) Z). 

***Right adjusted. 


TABLE III 

INPUT FORMAT 


TRANSIENT 1 * 2/02 COMPUTER PROGRAM: DIFFUSION-CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN Y 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

1-17 

Cards 
in va 

1-17 folic 
lues on care 

m exactly 
Is 1, 2, 7 

the same format as cards 1-17 in Table I, Note, however, the differences 
and 8 in the detailed discussion of diffusion-controlled input. 

Cards 18-25 are cards punched out by the program, at the end of the final non-diffusion-controlled run. 

18 

1 

12* 

1-2 

NAXT 

- - 

Number of Z’s in non-d.c. runs 

19 


8E10.5 

1-10 

11-20 

71-80 

ZT(I) 

z 

Axial stations in non-d.c. runs 

20a 
20b 
etc . 

** 

8E10.5 

It 

ptt(i) 

<y s 

Particle surface temperatures at end of non-d.c. 
runs 

21a 

21b 

etc. 

** 

8E10.5 

tt 

C02T(l ) 

Ci °2 

Interstitial concentrations of oxygen at end of 
non-d.c, runs 

22a 
22b 
etc . 

#* 

8E10.5 

11 

CH2T(I) 


Interstitial concentrations of hydrogen at end 
of non-d.c. runs 


*Right adjusted. 

**Enough cards should he used to contain (nAxt) values at the rate of 8 per card 




TABLE III (Coat’d) 

INPUT FORMAT 


TRANSIENT H 2 /O 2 COMPUTER PROGRAM: DIPFUS ION -CONTROLLED RUN 


CARD 

NO. 

NO. OF 
CARDS 

FORTRAN Y 
FORMAT 

COLUMNS 

USED 

SYMBOL OR 
DESCRIPTION 

CORRESPONDING 
SYMBOL USED 
IN EQUATIONS 

NOMENCLATURE 

23 a. 

23 h 

etc. 

* 

8E10.5 

1-10 

11-20 

niso 

CH20T(I) 

C H 2° 

''l 

Interstitial concentrations of water at end of 
non-d.c. runs 

24a 

24b 

etc . 

* 

8E10.5 

It 

TTI(I) 

T i 

Interstitial temperatures at end of non-d.c. runs 

25 

1 

2E10.5 

1-10 

11-20 

TIME 

TW 

t 

T W 

Cumulative time at end of non-d.c. runs 
Wall temperature at end of non-d.c. runs 

26 

1 

2E10.5 

1-10 

11-20 

CIH2 ( 1 ) 

cioe(i) 

c H 2 
c i°2 

Inlet concentration of hydrogen after "shutdown" 
Inlet concentration of oxygen after "shutdown" 


■*Enough cards should he used to contain (NAXT) values at the rate of 8 per card 




CODING OF A SAMPLE DATA CASE 


TRANSIENT H2/02 COMPUTER PROGRAM: Initial Non-Diffusion-Controlled Run 



K9 10962- 13 1 FIG. 1 


LISTING OF INPUT DATA PUNCH CARDS:. Sample Input 

TRANSIENT H 2 /0 2 COMPUTER PROGRAM 
(Initial Non— Diffusion— Controlled Run) 


Card No. 


TRANSIENT H2/02 CASE 13 *-+ 

(P=3Q0 

♦ G— 1.05 

* TF-360 * 

0/F=l * 25-30 MSH>TT 

1 

.11... 10 

20 . i.. _ 

. .1 ,._.o 

__ __ _ b __ 




2 

1.05 

492, 

360. 

300. 

0. 

,1472 

.1472 0. 

3 

360, 

300. 

.1472 

0. 

0, 

0. 

360. 0.5 

4 

.0047 

0.12 

,0066 

,0358 

360. 

0 • 

.00011 0. 

5 

0. 

.00001 

,001 

.002 

.004 

.006 

.008 .010 

7a 

.012 

.015 

,020 





7b 

0. , 

1. . . 

. 4 *. — .. — 0., 

„ _ 




9 

0. 

.01 

.03 ,05 





10 

.001 

.001 

.001 ,001 





11 

0. 

1. 

4 » 0 . 





12 

0. 

.01 

.03 .05 





13 

2100, 

2100. 

2100. 2100. 





14 

..0#\ . 

1 . 

4, .0 ... 





15 

0.' 

.01 

.03 ,05 





16 

.41 

.41 

.41 .41 





17 

,2944 

0. 






18 


K910962-13 FIG. 2 



LISTING OF INPUT DATA PUNCH CARDS: Sample Input 
TRANSIENT Ho/O n COMPUTER PROGRAM (Continuation of Non-Diffusion-Controlied Run) 


** transient h?/o 

>2 CASE 13 ** 

( P = 3 O 0 

11 10 

20 1 

3 0 


1,05 

492 , 

360 , 

300 . 

360 , 

300 . 

,1472 

0 ♦ 

.0047 

0.12 

.0066 

.0358 

0 . 

,00001 

.001 

,002 

,012 

.015 

,020 


0 .' 

1 . 

4 . 0 • 


0 , 

.01 

.03 .05 


.001 

.001 

.001 .001 


0 . 

1 , 

4 , . 0 * 


0 , 

.01 

.03 , 05 


2100 , 

2100 . 

2109 . 2100 . 


0 , 

1 . 

4 , 0 » 


0 , 

.01 

.03 ,05 


.41 

.41 

.41 ,41 


,00000 

. 10000 - 

*)4 * 100 . 00-02 . 

..£ 0000-02 


,12000-01 ,15000-01 ,20000-01 
M ,43575+03 .43575+05 ,45520+05 .47044+03 

c* ,43539+03 ,41906+03 ,39722+03 

,36026+03 ,38809+05 ,41217+05 .44521+03 
.42447+03 ,40009+03 

. .,45379-00. . 4 84 7.3+ 02_ ,5110.2.+ 0 JL51&5+ 0.3 


.10525+03 

,52932+92 


.30000+03 

.30000+03 

.30000+03 ,30000+03 

,29998+03 

.29908+03 


.14710-00 

,13630-00 

.12337-00 ,11852+00 

,12618-00 

.13376-00 


,14703-00 

,12992-00. 

, 1 15.76+.0.G. 103 06+ Q 0 

,74098-01 

.73103-01 


.ooooo 

.00000 

.ooooo .ooooo 

,00000 

.00000 


.ooooo 

.ooooo 

.00000 .00000 

.00000 

.00000 


.00000 

.42134+03 

,43569+03 .,44666+03 

,42554+03 

,41223+05 

.39352+03 

.ooooo 

.00000 

.00000 .00000 

.ooooo 

,00000 

♦OOOOO .ooooo 

,00000 

.000-00 

.00000 .00000 


Card No. 


G=1.05 

+ TF~36G 

* O/F-l * 25- 

30 MSH) TT 

1 

p 

0. 

.1472 

.1472 

0. 

CL 

3 

0, 

0. 

360, 

0.5 

4 

360, 

0. 

,00011 

0. 

5 

.004 

,00b 

.008 

.010 

7a 


7b 

9 

10 

11 

12 

13 

14 

15 

1 6 
17 


40000-02 

.60000-02 

.80000-02 

.10000-01 

18 a 

18 b 

48000+03 

.47417+03 

.46072+03 

.44732+03 

19 a 

19 b 

45768+03 

,45fe6l+03 

.44953+03 

.44040+03 

20a 

20b 

17444+03 

,17235+03 

.15794+03 

,13892+03 

21a 

21b 

29999+03 

.29999+03 

.29999+03 

,29999+03 

22a 

22b 

11599+00 

.11677+00 

.11892+00 

.12155+00 

23 a 

23 b 

88582-01 

,80914-01 

.76944-01 

,75017-01 

24a 

24b 

ooooo 

,00000 

.00000 

.00000 

25 a 

25 b 

OOOOO 

,00000 

.00000 

.00000 

26 a 
26 b 

45616+03 

,45292+03 

.44431+03 

.43466+03 

27 a 

27 b 

ooooo 

.00000 

,00000 

.00000 

28 a 

ooooo 

,00000 

.00000 

.00000 

28 b 

ooooo 

.00000 

.ooooo 

.ooooo 

28 c 
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LISTING OF INPUT DATA PUNCH CARDS (Cont.) 


ro 

-<i 


*00000 

.40655+03 

.42064+03 

,45805+03 

,51044+03 

.40781+03 

,42363+03 

,46704+03 

.49607+03 

,40777+03 

,42744+03 

,47696+03 

,39994+03 

, 40561+03 

,42712+03 

,47733+03 

,39631+03 

,40375+03 

,42641+03 

.47243+03 

♦39331+03 

.40215+03 

.42516+03 

.46143+03 

.39090+03 

,40078+03 

.42334+03 

.30580+03 

.30899+03 

.39962+03 

.42094+03 

,30184+03 

,38557+03 

,39557+03 

.41250+03 

,37578+03 

.37919+03 

.38658+03 

,39618+03 


,00000 
.40911+63 
,42361+03 
, 46525+03 
.50194+03 
,40879+03 
• .42720+05 
, 47506+03 
.45520+03 
.40916+03 
.43165+03 
.48559+03 
.40004+05 
.40723+05 
,43155+03 
.48528+03 
.39665+03 
.40556+03 
,43087+03 
.47722+03 
,39387+05 
.40409+03 
,42944+03 
,46072+03 
,39163+05 
.40280+03 
.42720+03 
.38600+03" 
.38985+03 
.40166+03 
,42420+03 
.38202+03 
.38643+03 
.39736+03 
.41480+03 
• 376Q1+03 
.37989+03 
.38776+03 
.39722+05 


,00000 

.40992+03 

.42700+03 

.47301+03 

.48435+03 

*41001+03 

.43124+03 

, 48353+03 

,40415+03 

,41078+03 

.43637+03 

,49432+03 

, 40033+ 03 

,4.0911+03 

.43647+03 

.49221+03 

.39716+03 

,40761+03 

,43578+03 

,47948+03 

.39457+03 

,40626+03 

,43404+03 

,38881+03 

, 39248+03 

,40502+03 

,43122+03 

,38614+03 

.39062+03 

,40388+03 

,42745+03 

,38228+03 

,38740+05 

,39926+03 

,41691+03 

,37629+03 

,38065+03 

.38898+03 


*00000 

.41098+03 

.43085+03 

.48122+03 

*43575+03 

.41150+03 

*43577+03 

.49217+03 

.40424+03 

.41271+03 

.44164+03 

.50241+03 

.40072+03 

.41125+03 

.44193+03 

.49649+03 

.39782+03 

.40993+03 

.44114+03 

,47417+03 

.39542+03 

.40868+03 

.43894+03 

.33901+03 

.39347+03 

.40747+03 

.43531+03 

.33636+03 

.39192+03 

.40629+03 

.43057+03 

.33262+03 

.38847+03 

.40128+03 

.41906+03 

.37663+03 

.33148+03 

.39022+03 


.00000 
,41230+03 
,43518+03 
• ,48968+03 
.40609+03 
,41326+03 
.44083+03 
,50054+03 
,40457+03 
,41492+03 
,44748+03 
,50809+03 
,40133+03 
,41371+03 
,44795+03 
,49648+03 
,39864+03 
,41254+03 
,44696+03 
.39219+03 
,39642+03 
,41136+03 
.44405+03 
,38918+03 
,39460+03 
.41015+03 
,43932+03 
,38669+03 
.39316+03 
,40888+03 
,43326+03 
,38304+03 
.38966+03 
,40341+03 
,37525+03 
,37702+03 
,38238+03 
,39147+03 


,00000 

,41391+03 

.44004+03 

,49798+03 

.40620+03 

,41533+03 

.44648+03 

.50756+03 

.40503+03 

.41747+03 

.45395+03 

.50869+03 

.40209+03 

,41649+03 

.45453+03 

,48203+03 

.39963+03 

,41547+03 

.45318+03 

.39233+03 

,39758+03 

,41434+03 

.44925+03 

.38944+03 

.39589+03 

,41308+03 

.44301+03 

.38711+03 

.39454+03 

,41166+03 

,43539+03 

,38354+03 

.39096+03 

,40562+03 

,37540+03 

.37747+03 

,38334+03 

.39271+03 


.40805+03 

,41581+03 

.44545+03 

.50539+03 

.40655+03 

,41773+03 

.45271+03 

.51101+03 

.40574+03 

.42039+03 

.46103+03 

,50216+03 

.40305+03 

.41963+03 

.46168+03 

.39592+03 

.40080+03 

.41873+03 

.45973+03 

.39258+03 

.39891+03 

,41762+03 

,45423+03 

.38981+03 

.39734+03 

.41626+03 

.44595+03 

.38763+03 

.39607+03 

.41462+03 

.38154+03 

.38412+03 

.39237+03 

.40789+03 

.37547+03 

.37799+03 

,38436+03 

.39392+03 


.40819+03 
.41805+03 
.45145+03 
.51047+03 
.40705+03 
.42048+03 
.45958+03 
.50773+03 
• 40664+03 
.42370+03 
.46872+03 
.47044+03 
.40423+03 
.42316+03 
.46932+03 
.39607+03 
.40217+03 
.42238+03 
.46634+03 
.39287+03 
.40043+03 
..42123+03 
.45851+03 
.39030+03 
.39897+03 
.41968+03 
.44732+03 
.38826+03 
.39776+03 
.41772+03 
.38172+03 
.38480+03 
.39391+03 
.41019+03 
,37560+03 
,37856+03 
.38544+03 
.39511+03 


Card No. 

28d 

28e 

28f 

28g 

28h 

28i 

28j 

28k 

281 

28m 

28n 

28o 

28p 

28q. 

28r 

28s 

28t 

28u 

28v 

28w 

28x 

28y 

28z 

28 aa 

28bb 

28cc 

28aa 

28ee 

28ff 

28gg 

28hh 
28ii 
28 jo 
28kk 

2811 

28mm 

28nn 

2800 

28pp 
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LISTING OF INPUT DATA PUNCH CARDS (Cont.) 


Card No. 


o 


fO 


CD 


.00000 ,00000 . C 000 G .00000 

,00000 ,00000 ,00000 .00000 
,00000 ,00000 ,00000 .00000 
.00000 ,00000 ,00000 .00000 
.95392-02 ,94819-02 .93789-02 .92345-02 
,02431-02 ,80131-02 .77951-02 .76036-02 
.76468-02 , 81119-02 ,88809-02 .10082-01 
,34662-01 ,47984-01 ,65471-01 .11025+00 
,12340-01 .12177-01 ,11982-01 .11765.-01 
.10408-01 .10132-01 .98811-02 .96673-Q2 
,10280-01 ,11156-01 .12561-01 .14748-01 
.56584-01 .97675-01 ,15770-01 .15718-01 
♦i5o7o-01 .14845-01 ,14585-01 .14300“0l 
.12705-01 ,12423-01 ,12189-01 .12031-01 
.14187-01 ,15946-01 .18667-01 .22802-01 
,22816-01 .22757-01 .22704-01 .22618-01 
.21745-01 .21488-01 .21212-01 .20920-01 
,19490-01 • 19307-01 .19210-01 .19241-01 
,24189-01 .27318-01 .31938-01 .33630-01 
.28838-01 .28782-11 .28708-01 .£8616-01 
,27897-01 ,27721-01 ,27542-01 .£7367-01 
,26968-01 ,27123-01 .27441-01 .£7983-01 
,38464-01 .43771-01 .51024-01 .72497-01 
.33692-01 .33659-01 .33620-01 .33575“01 
.33362-01 ,33347-01 .33356-01 .33398-01 
.34698-01 ,35386-01 .36340-01 .37647-01 
,55153-01 .71597-01 .37485-01 .37442-01 
,37518-01 .37552-01 ,37595-01 .37649“01 
.38200-01 .38405-01 .38661-01 .38982-01 
.42391-01 .43724-01 ,45438-01 .47644-01 
.40118-01. .40081-01 .40102-01 .40138-01 
.40563-01 .40704-01 ,40868-01 .41059-01 
,42561-01 .43023-01 ,43565-01 ,44206**01 
.50011-01 .52070-01 ,54624-01 .57798-01 
.42501-01 .42575-01 .42679-01 ,42813-01 
.43971-01 . 443jl-0l .44692-01 .45118-01 
,48099-01 .48922-01 .49852-01 .50909-01 
.59297-01 ,61971-01 ,65158-01 .73348-01 
,41931-01 ,42156-01 .42431-01 .42758-01 
,45192-01 .45849-01 .46566-01 ,47347"01 
,52347-01 .536*1-01 ,54966-01 .56489-01 
.66790-01 .72823-01 


,00000 .00000 .00000 .00000 
,00000 .00000 .00000 • .00000 
,00000 .00000 .00000 .00000 
,00000 .00000 .95972-02 .95618-02 

.90633-02 ,88778-02 .86808-02 .84691-02 
.74497-02 .73468-02 .73186-02 .74024-02 
.11921-01 .14712-01 .18932-01 .25276-01 
,12618-01 .12573-01 .12536-01 .12461-01 
.11527-01 .11268-01 .10989-01 .10698-01 
.95121-02 .94474-02 .95155-02 ,97691-02 
.18091-01 .23146-01 .30693-01 .41636-01 
,15669-01 .15581-01 .15451-01 .15282-01 
,13993-0.1 .13670-01 .13340-01 .13015-01 
.11984-01 .12097-01 .12433-01 .13084-01 
,29006-0.1 ,38104-01 .50854-01 .87845-01 
,22503-01 ,22358-01 .22183-01 .21978-01 
.20617-01 .20310-01 .20011-01 .19732-01 
.19455-01 .19925-01 .20759-01 .22108-01 
.48038-01 .76389-01 .28929-01 .28873-01 
.28505-01 .28375-01 .28228-01 .28067-01 
,27205-01 ,27067-01 .26969-01 .26927-01 
,28833-01 .30111-01 .31979-01 .34662-01 
,33806-01 .33756-01 .33741-01 .33719-01 
,33528-01 .33479-01 .33433-01 .33393-01 
,33485-01 .33634-01 .33866-01 .34208-01 
,39428-01 ,41841-01 .45089-01 .49428-01 
,37446-01 .37455-01 .37470-01 .37491-01 
,37716-01 .37801-01 .37906-01 .38037-01 
,39384-01 ,39892-01 .40534-01 .41350-01 
,50484-01 .54133-01 .58797-01 .71775-01 
,40189-01 .40257-01 .40341-01 .40442-01 
.41279-01 ,41533-01 .41826-01 .42166-01 
.44968-01 .45884-01 .46993-01 .48347-01 
,61738-01 .72349-01 .42484-01 .42458-01 
,42978-01 .43176-01 .43406-01 .43670-01 
,45592-01 .46120-01 .46709-01 .47365-01 
,52118-01 ,53514-01 .55139-01 .57045-01 
,41568-01 .41558-01 .41632-01 .41757-01 
.43137-01 .43569-01 .44054-01 .44595-01 
.48194-01 ,49113-01 .50107-01 .51182-01 
,58140-01 .59964-01 .61989-01 .64251-01 


29a 

29b 

29c 

29d 

29e 

29f 
29g 
29h 
29i 
2 93 
29k 
291 
29m 

29n 

29o 

29P 

29q 

29r 

29s 

29t ' 

29u 

29v 

29w 

29x 

29y 

29z 

29aa 

29bb 

29cc 

29dd 

29ee 

29ff 

29gg 

29hh 

29ii 

2933 

29kk 

2911 

29mm 

29nn 

2900 

29PP 


0 
NO 
ON 
IO 

1 

CO 


21 

p 

CO 

o 



LISTING OF INPUT DATA PUNCH CARDS (Cont.) : 

Card No. 


.00000 

•OOOoO 

.00000 

.00000 

,00000 

.00000 

.00000 

.00000 

30a 

,00000 

,00060 

. C -0000 

.00000 

,00000 

.00000 

.00000 

.00000 

30b 

,00000 

.00000 

.coooo 

.00000 

,00000 

,00000 

.00000 

.00000 

30c 

,00000 

,00000 

♦ 0000 

.00000 

,00000 

.00000 

.13658-00 

.13657-00 

30a 

.13656-00 

,13658-00 

.13658-00 

.13659-00 

,13660 -On 

. 136 bl-no 

.13662-00 

.13664-00 

30e 

, 13666-00 

-.13668-00 

,13670-00 

• 13673-00 

,13677-00 

,13681-00 

.13686-00 

.13692-00 

30 f 

,13700-00 

,13709-00 

,13720-00 

.13734-00 

,13753-00 

. 13778-00 

.13812-00 

.13860-00 

30 g 

.13926-00 

.14017-00 

.14133-00 

.14424-00 

,13233-00 

,13233-00 

.13232-00 

.13230-00 

30 h 

,13229-00 

.13226-00 

.13224-00 

.13221-00 

,13217-00 

.13213-00 

.13209-00 

.13204-00 

301 

.13196-00 

.13193-00 

♦13187-00 

.13180-00 

.13174-00 

.13168-00 

.13162-00 

.13157-00 

30a 

,13153-00 

.13150-00 

.13151-00 

.13157-00 

,13170-00 

,13193-00 

.13232-00 

,13292-00 

30k 

,13375-00 

, 13610-00 

.12923-00 

.12922-0 0 

,12920-00 

,12918-00 

.12914-00 

.12910-00 

301 

,12905-00 

,12899-00 

*.12892-00 

.12884-00 

,12875-00 

,12865-00 

.12853-00 

.12841-00 

30 m 

.12826-00 

.12815-00 

.12800-00 

.12785-00 

,12769-00 

,12754-00 

.12739-00 

,12725-00 

30n 

,12713-00 

.12705-00 

.12703-00 

,12709-00 

,12727-00 

.12764-00 

.12822-00 

.13008-00 

30o 

,12554-00 

,12553-00 

.12550-00 

.12546-00 

,12541-00 

.12534-00 

.12526-00 

.12516-00 

30 p 

,12504-00 

.12491+00 

.12477+00 

,12460+00 

,12442+00 

.12422+00 

.12400+00 

.12377+00 

30 q. 

,12353+00 

.12327+00 

*12300+00 

.12272+00 

,12243+00 

.12215+00 

. 12188+00 

.12164+00 

30r 

.121^3+00 

.12129+00 

.12123+00 

,12130+00 

,12154+00 

.12260+00 

.12472+00 

.12470+00 

30 s 

,12468+00 

.12464+00 

,12458+00 

.12450+00 

,12441+00 

.12431+00 

.12418+00 

.12404+00 

30 t 

,12388+00 

.12370+00 

.12350+00 

.12329+00 

,12306+00 

.12281+00 

.12254+00 

.12226+00 

30u 

,12197+00 

,12166+00 

,12135+00 

.12104+00 

.12074+00 

.12046+00 

.12020+00 

.11999+00 

30 v 

.11985+00 

.11981+00 

♦11988+00 

.12047+00 

,12541-00 

.12540-00 

,12538-00 

.12534-00 

30w 

,12526-00 

,12521-00 

.12512-00 

.12502"00 

,12490+00 

.12477+00 

.12461+00 

,12445+00 

30 x 

,12426+00 

, 124^6+00 

.12383+00 

,12360+00 

,12335+00 

.12308+00 

.12281+00 

.12252+00 

30 y 

,12223+00 

.12193+00 

.12165+00 

.12137+00 

,12112+00 

.12091+00 

.12075+00 

.12065+00 

30 z 

.12065+00 

.12097+00 

.12686-00 

,12685-00 

,12682-00 

.12679-00 

.12674-00 

.12667-00 

30aa 

,12660-00 

.12650-00 

.12639-00 

*12627-00 

,12613-00 

.12598-00 

.12581-00 

.12563-00 

30bb 

.12543-00 

,12522-00 

.12499+00 

.12475+00 

,12451+00 

.12425+00 

.12399+00 

.12373+00 

30cc 

,12347+00 

.12322+00 

,12300+00 

.12280+00 

,12264+00 

,12253+00 

.12249+00 

.12266+00 

3oaa 

.12863-00 

.12862-00 

• l286.CirrOQ 

,12857-00 

,12852-00 

,12846-00 

.12839-00 

.12831-00 

30ee 

,12821-00 

.12810-00 

,12798-00 

.12785-00 

,12770-00 

,12753-00 

.12736-00 

.12717-00 

30 ff 

,12697-00 

.12676-00 

,12655-00 

.12632-00 

,12609-00 

,12587-00 

.12564-00 

.12543-00 

30 gg 

,12523-00 

• 1250-5-00 

.12490+00 

.12480+00 

,12474+00 

,12483+00 

• 13190-00 

.13189-00 

30hh 

,13187-00 

.13185-00 

.13181-00 

.13176-00 

,13171-00 

,13164-00 

• 13156-00 

,13147-00 

3011 

.13137-00 

.13126-00 

.13114-00 

.13101-00 

,13087-00 

,13073-00 

.13057-00 

.13040-00 

30 j a 

,13023-00 

. 13005-00 

,12987-00 

,12969-00 

,12952-00 

,12935-00 

.12919-00 

.12905-00 

30kk 

.12893-00 

.12884-00 

,12878-00 

.12881-00 

,13780-00 

,13779-00 

.13778-00 

.13776-00 

3011 

.13774-00 

.13770-00 

.15766-00 

.13762-00 

,13756-00 

.13750-00 

.13743-00 

.13736-00 

30mm 
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LISTING OF INPUT DATA PUNCH CARDS «Eont.) 


,13727-00 

,13641-00 

,13565-00 

1 n i 1 

,29000-01 

,2944 


,13718-00 

-.13629-00 

.13563-00 

1 1 1 1 1 
,36046+05 

0 . 


.13709-00 

♦13617-00 

1 

.10000-02 


.13699-00 

.13605-00 


,13688-00 

,13594-00 


, 13677-00 
. 13585-no 


.13665-00 

.13576-00 


.13653-00 

.13569-00 


Card No. 

30nn 

30oo 

30pp 

31 

32 

33 


00 

o 
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LISTING OF INPUT DATA PUNCH CARDS: Sample Input 
TRANSIENT H2/O 2 COMPUTER PROGRAM (All Diffusion-Controlled Run) 

TRANSIENT H2/02 CASE 12 ** <P=3C0 * G=2*10 * TF=360 * 0/F=l * 14-18 MSH) TF 

24 20 20 1 4 i 

2*10 ' 492 . 360 # 300 # 0 # * 1,472 .1472 0 . 

360 , 300 * .1472 0 * 0 , 0 . 360 . 2.0 

, 004" 7 0.12 ,'0066 ,0358 360 . 0 . ,00011 0 . 

0 * #0002 .0004 .0006 .0008 ,0010 .0012 .0015 

.0020 .0025 ,0030 .0040 ,0050 .0060 .0070 .0080 

#0090 ,0100 .DUO ,0120 ,0140 .0160 .0180 .0200 

#001 ,001 #001 .002 .002 .002 .005 .005 

#005 ,005 .010 .010 #020 .020 .020 ,020 

,'850 .050 .050 ,10 

0 # 1 . 4 . 0 . 

0 # .01 *05 #05 

•0016 #0016 *0016 #0016 
0 # 1 , 4 , 0 , 


*01 
♦ 0016 

1. 


.03 
* 0016 

^ * 


#05 

#0016 

0 # 


0, 

.01 

.03 

• 05 

1300, 

1300. 

1300, 

1309. 

0. 

1. 

4. 

0, 

0. 

.01 

*03 

,05 

• 41 

.41 

.. *41 . 

.ju41. 


.00000 .lQQOO-02 .2.0000-02 .40000“02 ,60000“02 

.15000-01 .20000-)! 

.49425+03 .49425+05 ,49579+Co .50970+03 ,52166+03 
.50783+03 , 49804+05 

.14720-00 .12909-00 .11385+00 *83390-01 .69297-01 
.26540-01 .16162-01 

,14720-00 .13956-00 *13330-00 .12282+OU .11501+00 
.10247+00 .10149+00 

.00000 .13505-01 .25079-01 .44412-01 ,58977-01 

.97996-01 .11017+00 

.36000+03 ,37705+03. ,39226+03 *42021+03 .44399+03 
.48514+03 .48626+03 
.31000-01 .36067+03 
,2944 0. 


.80000-02 .10000-01 .12000-01 
.51481+03 .51263+03 .51010+03 
,55351-01 .44606-01 .36167-01 
.110i3+00 .10683+00 .10460+00 
,70682-01 .80296-01 .88311-01 
,45987+03 .47091+03 .47834+03 
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Output Description 

Output from the B^/Og transient program consists of both printout and punch cards. 
The printing is basically the same for all of the three types of runs discussed under 
the input description; however, the form of the punch card output varies. 

Printout 


There are two different types of printed output from this program: (l) summary 

printing (when input indicator DEBUG=F) , and (2) detailed printing (when DEBUG=T) . 

The detailed form yields a considerable amount of printed output since everything • 
produced by the summary output is printed out as well as a full station-by- station 
and time-by-time account of reactor calculations. This can yield over 2000 lines 
of printing per pass. The detailed printing should be used only when the summary 
form doesn't give enough detailed information (as for example, at the very beginning 
of an initial non-diffusion-controlled run for fast reacting cases), where the 
summary printing for only every twentieth time step may not allow enough information 
for plotting. An abridged sample of the detailed output from a transient run 

is shown in Figs. 5a through 5t. The output shown in the top of Fig. 5a and in Figs. 
5e - 5t is what would have been produced if the summary printing alone had been 
requested. A description of the printed output produced for both the detailed and 
summary forms follows, along with the name of the routine from which it is printed. 

Detailed Form 


F874: 

A listing of the program input including the data case title (see top of Fig. 5& ). 
The extra punch card input for the continuation and all diffusion-controlled runs is 
not printed. 

The following message which concludes each run (see Fig. 5t): 

t 

■ x - x - x - x - x x x - OPERATIONS COMPLETE - x-x - x - x - x - x - x - x - 
TOTAL EXECUTION TIME = XX. XX MINUTES 


Subroutine MCALC: 

The following message which is printed when the non-diffusion-controlled phase 
of a data case is completed (see Fig. 5t); 

xxxxxxxxxxxxxxxx x - x - x - x - xx - x - x - x - x - x - x-x - xx - xx - x - 
GO - TO -ALL -DIFFUSION-CONTROLLED RUN 
TIME = XX.XX SECONDS 
X X X X X X - X - X X X XX X X X X X - X x x x - x x - x - x - x X - X - X - X - X X X X - 


Subroutine DEBUG3: 

1. A complete list of the following values for each axial station and every time 
increment of each pass (see Figs. 5b - 5&) : 

32 
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o 

a. interstitial concentrations of hydrogen in lb /ft 

b. interstitial concentrations of oxygen in lb/ft3 

c. interstitial concentrations of water in its solid phase in lb/ft^ 

d. interstitial concentrations of water in its liquid phase in Ib/ft^ 

e. interstitial concentrations of water in its vapor phase in lb/ft3 

f. interstitial concentrations of helium in lb/ft3 

g. interstitial temperatures in deg R 
*h. particle surface temperatures 

i. enthalpies in BTU/lb 

j . interstitial pressures in psia 

K. corresponding axial positions in ft. 

2. The cumulative time (sec) and wall temperature (deg R) at the end of every time 
step. (See Rigs. 5b - 5d) . 

Subroutine PRPLOT: 

Printer plots of the following values -vs. axial position (in feet), corresponding 
to the last time step of each pass. Each plot summarizes the last five passes or 
whatever passes remain. (See Figs. 5e - 5k). 

a. interstitial temperatures (deg R) 

*b. particle surface temperatures (deg R) 

c. chamber pressures (psia) 

d. hydrogen mole-fractions 

e. oxygen mole-fractions 

f . water mole- fractions 

g . helium mole- fractions . 

Subroutine OUTPUT: 

Summary printing of the seven variables mentioned under subroutine PRPLOT vs. 
axial position (ft) at the end of every pass. This printing is produced after every 
tenth pass and/or at the end of the last pass. (See Figs. 5-£ - 5^) . 

Subroutine TMPENT: 

Summary printing of wall temperature (deg R) and interface location (axial position) 
in feet of water ehange-of-phase) vs. time (values printed correspond to the end of 
each pass). (See Fig. 5s). Also, if H2LEAD=F, pressure in psia and mass flow rate 
in lb/ft 2 -sec vs', time are printed. This printing is done after every tenth’ pass 
and/or at the end of the last pass. 

*Irregularities in the values for particle surface temperature, caused by water 
change- of- phase, may be exaggerated in some cases due to use of too few axial positions. 
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Summary Form 


Everything printed under the detailed printing described above is also produced 
for the summary form except those things printed from subroutine DEBUG3. 

Punch Cards 


Two different types of punch card output can be produced by this program. For 
the all diffusion-controlled runs, no punch cards will be produced at all. For 
both of the two types of non-diffusion-controlled runs (initial and continuation) 
the two possible types of punch card output are described below. 

Continuation Cards 


These cards are described in detail under the input discussion section (cards. 
18-32 under "continuation run"). They will be produced when HJNCH=1 or 3 and when 
a run completes 3JPASS number of passes without reaching the all diffusion- controlled 
phase . 

Diffhsi_on-Controll^ed_Cards_ 

These cards are described in detail under the input discussion of all diffusion- 
controlled runs (cards 18-25). They will be produced whenever the non-diffusion- 
controlled phase of a data case is completed. 
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LISTING OF OUTPUT FORSAMPLE DATA CASE_ 

***** PROGRAM INPUT ***** 


TRANSIENT ANALYSIS OF H2/02 CATALYTIC REACTOR 
** TRANSIfNT H2/02 CASE 13 ** (P=300 * G=1.05 * TF=360 * 0/F=l * 25-30 MSH) 


NOFZ NPASS INPERP DEBUG MOONO PUNCH 

11 10 ‘20 T 1 1 

\ 

G Ts MH2' M02 MH20 MHE D0H2 ' D002 DOH20 KP ALPHA 

,10500+01 .49200+03 .20200+01 ,32000+02 ,ld020+02 .40030+01 .16760-02 .75720-03 .63550-03 .40000-04 ,16000+14 

AUPHAK GAMMA 

.31740+08 .50000+04 

AW CW MW DC TA HA HA1 HA2 

.47000-02 .12000+00 .66000-02 .35800-01 .36000+03 .00000 .11000-03 .00000 


INITIAL TEMPERATURES, PRESSURES, AND CONCENTRATIONS 
T • P CIH2 CI02 . CIH20 ClHE TPS 

.36000+03 ,30000+03 .14720-00 ,00000 .00000 .00000 ,36000+03 


INLET TEMPERATURE, PRESSURE, ENTHALPY, AND CONCENTRATIONS 
T P H CIH2 CI02 CIHE 

.36000+03 ,30000+03 .00000 .14720-00 .14720-00 .00000 


lo INPUT AXIAL STATIONS 

^,00000 .10000-04 .10000-02 .20000-02 ,40000-02 .60000-02 .80000-02 .10000-01 .12000-01 .15000-01 .20000-01 


A, AP, DELTA VS 2 TABLES 


,U0000 

,10000-02 

.10000-01 

.10000-02 

.30000-01 

.10000-02 

.00000 

.50000-01 

.10000-02 

.10000+01 ,40000+01 

,00000 




00000 

21000+04 

,10000-01 

,21000+04 

.30000-01 

.21000+04 

.00000 

.50000-01 

.21000+04 

.10000+01 .40000+01 

.00000 




00000 

41000-00 

.10000-01 

.41000-00 

.30000-01 

.41000-00 

.00000 

,50000-01 

.41000-00 

.10000+01 .40000+01 

.00000 




1472U-00 

15858-00 

.14722-00 

,14898-00 

INTERSTITIAL CONCENTRATIONS OF 
.15034-00 .15229-00 .15357-00 

H2 

,15444-00 

.15507-00 

.15553-00 

.15604-00 

14720-00 

«l4b89-00 

,11922+00 

INTERSTITIAL CONCENTRATIONS OF 
.97876-01 .67544-01 .48046-01 

02 

.34870-01 

.25661-01 

,19071-01 

.12376-01 


.01781-02 P 

CJ1 

INTERSTITIAL CONCENTRATIONS OF H20(S) 
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. uouoo 
.oooou 

.00000 

.00000 

.00000 .00000 .00000 

.00000 

,00000 

.00000 

.uooou 

.OOOOU 

.ooooo 

.00000 

INTERSTITIAL CONCENTRATIONS OF H20(L> 
.00000 ,00000 .00000 .00000 

.00000 

.00000 

, uouoo 

.UOUOO 

.00000 

,00000 

INTERSTITIAL CONCENTRATIONS OF H20CV) 
.OOOOO ,00000 .00000 .OOOOO 

.00000 

.00000 

,00000 

,00000 

.OOOOO 

.00000 

INTERSTITIAL CONCENTRATIONS OF HE 
,00000 .00000 ,00000 

,00000 

.00000 

.00000 

,36000+03 

,36000+03 

.36000+03 

,36000+03 

INTERSTITIAL TEMPERATURES 
.36000+03 .36000+03 .36000+03 

,36000+03 

.36000+03 

.36000+03 

,36000+03 

,36000+03 

.36000+03 

.36000+03 

PARTICLE SURFACE TEMPERATURES 
.36000+03 ,36000+03 .36000+03 

.36000+03 

.36000+03 

.36000+03 

,00000 

,00000 

.ooooo 

,00000 

ENTHALPIES ' 

.00000 .00000 ,00000 

,00000 

.00000 

.00000 

,30000+03 

,29996+03 

.30000+03 

.30000+03 

INTERSTITIAL PRESSURES 
,30000+03 ,30000+03 ,29999+03 

.29999+03 

.29999+03 

.29999+03 

,00000 

,20000-01 

.10000-04 

.10000-02 

AXIAL POSITIONS 

.20000-02 .40000-02 .60000-02 

,60000-02 

.10000-01 

.12000-01 


TIME s 

,20000-02 

WALL TEMP = .36000+03 




,14720-00 

,15543-00 

.14720-00 

.14756-00 

INTERSTITIAL CONCENTRATIONS OF H2 
,14814-00 .14964-00 .15098-00 

.15209-00 

.15290-00 

.15370-00 

14720-00 

15674-01 

.14716-00 

.13862-00 

INTERSTITIAL CONCENTRATIONS OF 02 
.12832-00 .10386+00 .82945-01 

.65949-01 ' 

.52380-01 

,41606-01 

00000 

00000 

.00000 

,00000 

INTERSTITIAL CONCENTRATIONS OF H20CS) 
,00000 .00000 .00000 ,00000 

.00000 

.00000 

OOUOO 

00000 

.00000 

.OOOOO 

INTERSTITIAL CONCENTRATIONS OF H20(L> 
.00000 ,00000 .00000 ,00000 

.00000 

.00000 

ooooo 

00000 

.00000 

,00000 

INTERSTITIAL CONCENTRATIONS OF H2CHV) 
,00000 .00000 .00000 .00000 

.00000 

.00000 

UOUOO 

uouoo 

.00000 

.00000 

INTERSTITIAL CONCENTRATIONS OF HE 
,00000 ,00000 ,00000 

.00000 

.00000 

.00000 


.ooooo 

.00000 

.00000 

.00000 

.36000+03 

.36000+03 

.00000 

,29990+03 

.15000-01 

.15453-00 
.29094-01 
.00000 - 
, 00000 ' 
.00000 
.00000 
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INTERSTITIAL temperatures 


o6uoy+o3 
> 36015+03 

,3b00l+n3 

,36038+03 

.36059+03 

*36069+03 ,36064+03 

.36055+03 

,36046+03 

.36038+03 

56159+03 

36012+03 

.3bl59+03 

.36133+03 

PARTICLE SURFACE TEMPERATURES 
,36113+03 ,36083+03 ,36062+03 

.36048+03 

.36037+03 

.36029+03 

UOOOU 

ooooo 

.92068-02 

.65744-00 

.10206+01 

enthalpies 

.12052+01 .11143+01 

,92575-00 

.71480-00 

.51278-00 

30000+03 

29998+03 

.30000+03 

,30000+03 

.30000+03 

interstitial pressures 

,30000+03 .29999+03 

.29999+03 

.29999+03 

,29999+03 

UOUOO 

20uQQ-01 

. 10000-04 

.10000-02 

.20000-02 

AXIAL POSITIONS 
.40000-02 .60000-02 

.80000-02 

.10000-01 

.12000-01 


TIME = 

,40000-02 

•WALL TEMP = 

,36000+03 





< ♦ 




* 


■ . 36027+03 
. 36021+03 
. 23538-00 
. 89999+03 
. 15000-01 
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cn 


,14720-00 
• 05301-01 

.14574-00 

.48469-01 

INTERSTITIAL CONCENTRATIONS OF H2 
.37510-01 .26211-01 .29236-01 

.32801-01 

.37604-01 

.42859-01 

.51447-01 

.14720-00 

,t>8054-03 

.14549-00 

.25569-01 

INTERSTITIAL CONCENTRATIONS OF 02 
.14965-01 .56506-02 .35308-02 

.24959-02 

.18540-02 

.13951-02 

.97574-03 

,00000 

,00000 

.00000 

.00000 

INTERSTITIAL CONCENTRATIONS OF H20(S) 
,00000 .00000 .00000 .00000 

.00000 

.00000 

.00000 

,00000 

,00000 

.00000 

.00000 

INTERSTITIAL CONCENTRATIONS OF H20(L) 
,00000 .00000 .00000 .00000 

.00000 

.00000 

.00000 

,00000 

,o7235-01 

.32600-03 

,29557-01 

INTERSTITIAL CONCENTRATIONS OF H20CV) 
.29075-01 .29144-01 ,33216-01 .39165-01 

.46204-01 

.53591-01 

.65236-01 

,U0Q00 

,00000 

,00000 

,00000 

INTERSTITIAL CONCENTRATIONS OF HE 
,00000 .00000 ,00000 

.00000 

.00000 

.00000 

.00000 

,36000+03 

,46241+03 

.36353+03 

.74477+03 

INTERSTITIAL TEMPERATURES 
.92759+03 .11766+04 .11279+04 

.99467+03 

.85695+03 

.74387+03 

.60733+03 

,11596+04 

45104+03 

.11596+04 

.14381+04 

PARTICLE SURFACE TEMPERATURES 
.12711+04 .13648+04 .11413+04 

,95012+03 

.79626+03 

.68695+03 

.57105+03 

00000 

20486+03 

.60869+01 

,66409+03 

ENTHALPIES 

.10224+04 .14719+04 .13838+04 

.11442+04 

.89827+03 

.69743+03 

.45599+03 

30000+03 

29991+03 

.30000+03 

.30000+03 

INTERSTITIAL PRESSURES 
,29999+03 .29998+03 .29997+03 

.29995+03 

,29994+03 

.29993+03 

.29992+03' 

00000 

20000-01 

.10000-04 

.10000-02 

AXIAL POSITIONS 

,20000-02 .40000-02 .60000-02 

.80000-02 

.10000-01 

.12000-01 

.15000-01 


TIME = 

.75000-01 

WALL TEMP = .37495+03 
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time profiles: H2 mole fraction vs axial position 

•97000-00 — + — — — + — — + — — — -+ — ■ +■ 


1 


1 


1 


1 


,95000-00 + 


1 


1 


1 


4 



,93000-00 + 


MOLE FRACTION 
(DIMENSIONLESS) 


2 


.91000-00 + 


3 

4 

2 


4 


.89000-00 + 2 

: 4 

: 2 

; 4 


. 87000-00 + — + — 

-.291 -.073 .145 , 364 ' ,582 


2 

4 2 

4 




-+ — + — 

.800 1.018 


4 


1.236 1.455 


4 


2 


! 

+ 

: 

i 

: 

i 


1,673 1.891 


— + 
2.109 


AXIAL POSITION (FT) TIMES (IN SEC) REPRESENTED BY CURVES 1 THRU 4 

.0290 .0490 ,0690 .0750 


x-aais scaling: x values*io** 2 


MS 'Did i UI-2960I6M 



• 59500—01 + 


time profiles: 02 mole fraction vs axial position 


.47500-01 + 


.35500-01 + 


MOLE FRACTION 
(OlMEhlSIONLESS) 


,23500-01 + 


.11500-01 + 


50000-03 + 

-.291 



AXIAL POSITION (FT) 


TIMES (IN SEC) REPRESENTED BY CURVES 1 THRU 4 
.0290 ,0490 .0690 .0750 


X-AXIS SCALING! X VALUES*10** 2 


K&10962— 13 ' ■ FIG.' Si 



,14000-00 +- 


TIME PROFILES: H20 MOLE FRACTION VS AXIAL POSITION 
+ + 


■x 

to 

$ 

CT5 

r 


.11000+00 + 


.60000-01 


MOLE FRACTION 
(DIMENSIONLESS) 


,50000-01 + 


. 20000-01 + 


'. 10000-01 

-.291 


.073 


.145 


.364 


.582 .800 1.018 

AXIAL POSITION (FT) 


1.236 


1.455 


-+ 

1.673 


1.891 


— + 
2.109 


TIMES (IN SEC) REPRESENTED BY CURVES 1 THRU 4 ^ 

,0290 .0490 .0690 .0750 5 


X— ax is scaling: x values*io** 2 



.10000-10 + 


TIME pfcOFU.ES: HE MOLE FRACTION VS AXIAL POSITION 
-+ ■ •— +- — - — : + +- 


.00000 + 


-. 10000-10 + 


MOLE FRACTION 
(DIMENSIONLESS) 


-. 20000-10 + 


-.50000-10 + 


40000-10 + 

-.291 



AXIAL POSITION (FT) 


TIMES (IN SEC) REPRESENTED BY CURVES 1 THRU 4 
.0290 .0490 .0690 .0750 


X-AXIS SCALING: X VALUES*10** 2 


, \ $ 1 - 296016 ) 



TIME PROFILES OF INTERSTITIAL TEMPERATURE VS AXIAL POSITION 
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R) 
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4 
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TIME PROFILES OF PARTICLE TEMPERATURE VS AXIAL POSITION 

AXIhL * ■ TEMPERATURE (DEGREES R) 

POSITION * 

(FT) * 
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TIME PROFILES OF CHAMBER PRESSURE VS AXIAL POSITION 
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TIME PROFILES 6F HZ MOLE FRACTION VS AXIAL POSITION 
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TIME PROFILES OF 02 MOLE FRACTION VS AXIAL POSITION 
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TIME PROFILES OF H20 MOLE FRACTION VS AXIAL POSITION 
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TIME PROFILES OF HE MOLE FRACTION VS AXIAL POSITION 
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WALL TEMPERATURE* INTERFACE VARIATIONS VS TIME 
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************************************** 

* SO TO ALL 0 1 FFUS I ON-CONTROLLED RUN * 

* TIME = .7500-01 SEC * 

********* **************** ************* 


****** OPERATIONS COMPLETE ****** 
TOTAL EXECUTION TIME: 4.0 MINUTES 
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Description of Routines 

The following is a list and brief description of all the routines in the H2/O2 
transient computer program;- 

F874: This routine handles the reading and writing of the standard 

deck input. Subroutine MCALC is called from it. 

MCALC: This routine handles the main body of reactor calculations. 

Most of the other subroutines are called from it. 

DERTV: This subroutine calculates the values of the derivative 

expressions. 

PCALCT: This subroutine does the non- diffusion^ controlled catalyst 

particle analysis. 

PCALCS: This subroutine does the diffusion-controlled catalyst particle 

analysis. 

GREER: This subroutine finds the values of the Green's functions used 

in calculations in subroutine’ PCALCT. 

ROOT: This subroutine finds the roots of the characteristic equation 

x*cot(x)+C=0. using a modified method of successive substitu- 
tions. These roots are used to find the Green's functions in 
subroutine GREEN. 

This subroutine inserts the interface axial positions (positions where 
water changes phase) into the original Z array. ' 

This subroutine deletes the inserted interface axial positions at 
the end of each time step. 

OVSTEP: This subroutine makes appropriate adjustments to temperature, 

pressure, concentrations, etc., when axial positioning has over- 
stepped the interface boundaries. 

UMAX: This routine is used to find the maximum value of elements in an 

array . 

MIN: This routine is used to find the minimum value of elements in an 

array. 


INSERT: 

DELETE: 
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RDPNCH: 


OUTHJT: 


SETUP: 


PRPLOT : 


TMPENT: 


DEBUG 3 : 


UUBAR: 


TABLES: 


This subroutine handles the punching and reading of the extra 
input needed by the continuation and all diffusion- controlled runs. 

This routine prints the temperature, pressure, and mole- fraction 
profiles -in summary form. 

This subroutine is used to setup the arrays needed for the printer 
plots in subroutine PRPLOT. 

This subroutine is used to develop and print the printer plots of 
temperatures, mole- fractions, etc. vs. axial position at various 
times. 

This subroutine handles the printing of wall and interface position 
values vs. time. For non-hydrogen-lead cases, it also handles 
printing of the mass flow rate and pressure vs. time. 

This subroutine handles the detailed printing which is produced 
only when input logical indicator DEBUG=T. 

This subroutine will do linear quadratic or cubic interpolation 
between tabulated points. Interpolation may be either univariate 
or bivariate. 

This subroutine contains tables of heat of reaction, species 
viscosities, species specific heats, vapor pressure, heat of con- 
densation, and species thermal conductivities. 
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LISTING OP COMPUTER PROGRAM 

LOGICAL DEBUG, H2LEAD~ 

INTEGER TITLE (12)* P.UNCH . . • 

REAL KAPPA » KP , KCH2 » KC02 , KCH20 # MH2 , M02 , MH20, MHE , MW , MU , MB AR » MBAR5S 
PARAMETER NAX=7Q, NT.I=4Q. 

PARAMETER NPP=30, NRTS=200 

COMMON /CB1/ .TPPJ.NP.Pf NAX),fC.Q2JNEP_f.NAX )..iCH2U«^ 

X T(NAX) fP(NAX) »PT(NAX) ,H(NAX) ,CIH2(NAX> ,CI02(NAX) »CIH20(NAX) # 

X . CIH2QSCNAX) »CIH2QL(NAX) »CIH20V (NAX).f CIHEJNAX) ,DELT (NTI) * 

X PPT(NAX) 

COMMON /C.B2/ Z(NAX) f ZSSL(NTl)_,ZSLUNTI.).*ZLLV(N.TI) ,ZLV.V(NTI) 

COMMON /CB3/ ALPHA t GAMMA , KAPPA , HC » KP , A , DELTT » KCH2 , KC02 , KCH20 , TCH2 # 
X • TC0.2 1 TCH.20f P.CH2 f .PCQ2f ECJH20i-MH2 MHEiHR a HA.rHAlf.MA2j 

X A W , CW , MW , DC , T A f D0H2 , D0Q2 , 00H2O , ALPHAK » AP , DELTA , G » MU » HS , HL , 

X HV » TW , RHOM * WH2 , WQ2 , WH20S , WH20L # WH20 V , WHE tCFBAR , MB AR , C I CPI » 

X C ICP2 ,C ICP3 , RHO , RHOV , DELHC t TlTPS t DTDZ » DHDZ * TS , TL , DPDZ » DC1DZ , 

X DC.2DZ r DC3S0Z , 0C3LDZ » 0C3.VDZf DC4DZ» RHQP , CP f DELTP » DPTPS * SDT I ME » 

X PF , TSS , AC , VC , PREST > GSS » MBARSS 

COMMON /CB4/ NQFZf.NP.ASS » NAXI ALi.lNPERPf.MODPCj.M.ODNQ*Pi^iCtU I.ISTEP* . 

X JDP f I SHUT » I HOT 

COMMON /TABLES/ AVSZ ( 40 } , APVSZ ( 40) , UELVSZ < 40 ) » HRVST (36 ) , TCi ( 22 ) , 

X TC2(20) »TC3V(2Q) »TC4(18) ,MUlVST(34) ,MU2VST(34) ,MU3VST(32 ) r 

X MU4VST ( 34 ) » CFI VST ( 34 ) # CF2VST ( 34 ) , CF3 VST ( 34 ) , CF3SVT (20 ) * 

X CF4VST (20) r VPVST (60) , TVSVP(48) ,DHCVST(24) 

DATA ALPHA f ALPHAK > GAMMA t KPjJMH2 j MQ2 *Ji1H20 iMHEtPCHS# PCQ2*RCH2Q / 

X .160 E+14, ,3174 £+8» .50 E+4, ,40 E-4, 2.02, 32. Or 18,02. 

X • 4.0Q3r 188,16. 730. 59» 3200.2 / 

DATA TCH2»TC02,TCH20,DOH2fD002,DOH20,RHOP»CP / 

X 59.94, 277. 92» 1164,6, .1676 E-2, ,7572 £-3, ,6355 E-3. 

X 124.9, .177 / 

KAPPA = Kpy.(RHOP*CP) 

0 CALL TWINS (TSTART) 

READ (5ilfEND=9V) TITLE, H2LEAD, DEBUG 

1 FORMAT (13A6,2L1) 

READ (5,3) NOFZ » NPASS flNPERP# MODNO » PUNCH iITSTEP 
3 FORMAT (2014) 

READ (5f2). .Gf.T5.rTJl ).»P.( 1 )jHJ 1 ) rjCIH211) rC 102.111 fCIHE.U). 

2 FORMAT (8E1Q.5) 

READ (5,2). TIfPJ,CH2I,C.02IfCH20I,CHEI#PTI.fS0.TIME 
HEAD (5,2) AW r CW , MW , DC » T A , HA » H A 1 , HA2 
N AXIAL = .NOFZ ■ 

IF (PUNCH. EG, 4) NAXI AL=N0FZ+4 
IF (H2LEAD) GO T.O 1,01 
GSS = G 

READ (5,102) PFfTSSfMBARSSf ACfVCfPREST 
102 FORMAT (6E10.5) 

101 READ (5f2) (ZJI), 1=1, NOFZ) 

IF ( ITSTEP,EQ»0) GO TO 7 
READ (5r2) (DELKID I=1»NPASS) . 

7 WRITE (6,10) TITLE 

10 FORMAT (*D 53X. ****** PROGRAM INPUT ****** /// 

X 43X * TRANSIENT ANALYSIS OF H2/02 CATALYTIC REACTOR* / 

X 3QX 14A6 ) 

WRITE (6,17) NOFZ, NP ASS, INPERP, DEBUG, MODNO, PUNCH 
17 FORMAT l*Q' .// 37X *NQFZ» 6X . *NPAS.SJ . 5X. • INPERP* 3X ♦ DEBUG* 7X 
X ’MODNO* 6X ’PUNCH* / 37X, 13, 9X> 12, 8X, 12, 9X, L2, 9X, 

X 12, IQXf ID 

WRITE (6,12) G , TS , MH2 , M02 , MH20 , MHE , D0H2 , DQQ2 , DOH20 , KP , ALPHA 
12 FORMAT ( *0* 12X *G* 10X *TS* 8X *MH2* 8X *M02* 8X *MH20* 7X * MHE * 
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x 8X *D0H2’ 7X *D002’ 6X ’DGH20’ 8X ’KP* 7X * ALPHA’ / 

X 6X, UEil.5 ) 

WRITE (6* 13) ALPHAK, GAMMA 

13 FORMAT (’O’ 9X ’ALPHAK* 6X ’GAMMA* Yj5X, 11 E.Ul».5J. _ 

WRITE ( 6 ,5 j AW,CW»MW,DC,TA,HA,HA1,HA2 

5 F0RMA' r (’O’ 27X * AW* 9X «CW* 9X *MW* 9X »DC» 9X *TA» 9X *HA' 8X 
X *HA1» ax *HA2* / 22X , 6E11.5/ ) 

WRITE (6,19) TI f PI»CH2IfC02IfCH20IfCHEIrPTI 
19 FORMAT (*0* 40X ’INITIAL TEMPERATURES » PRESSURES, AND CONCENTRATIO 

XNS* / 33X *T* 1QX *P» 8X *CIH2* 8X *CI02* 6* . * CJLH2Q.1 . 7.X _’CIHE* 

X ‘ 8X ’TPS* / 28X , 7E11.5 /) 

WRITE (6,14) T{1),P(1)»H(1),CIH2U),CI02<1),CIHEU> 

14 FORMAT (*0* 37X ♦ INLET TEMPERATURE, PRESSURE, ENTHALPY, AND CONCE 

XNTRaTIONS* / 39X *T* 10X *P* 10X *H* 9X *CIH2* 7X »CI02* 7X *CIHE* 
X / 33X, 6611,5 / ) 

WRITE (6,16) UUJ ,I=1,N0FZ) 

16 FORMAT (’O’ 56X ’INPUT AXIAL STATIONS’ / (12EH.5) ) 

IF (ITSTEP,EQ,0> 60 TO 8 
WRITE (6,18) (DELT ( I ) , 1=1 ,NPASS ) 

16 FORMAT (’O’ SIX ’TIME INCREMENTS FOR EACH PASS’ / (12EU.5) ) 

8 IF (H2LEAD) GO -TO 104 

WRITE (6,103) PF r TSSfMBARSS»ACfVC 

103 FORMAT (’O’ 54X ’VALUES USED TO CALCULATE DPDT* / 39X »PF* 12X 

X ’TSS* 10X ’MBARSS’ 10X *AC» 12X ’VC* / 35X,E10.b,4X» 

X ElO«5,4x,EiO »5,4X,E10 »5,4X,E10«5 /) 

104 READ (5,4) (AV5Z(I) ,1=1,4) 

NA = AVSZ ( 3 ) +4 ♦ 

!ib = NA+1 

NC = 2,*AVSZ(3)+4. 

BEAU (5,4) (AVSZ ( I ) , 1=5, NA) 

READ (5,4) (AVSZ(I) »I=NLS,NC) 

WRITE (6,24) 

24 FORMAT (’0* // 53X *A, AP, DELTA VS Z TABLES’) 

WRITE (6,21) (AVSZ(I) ,1=1,4) 

WRITE (6,22) (AVSZ(I) ,1=5, NA) 

WRITE (6,22) (AVSZ(I) ,I=NB,NC) 

WRITE (6,23) 

READ (5,4) (APVSZ(I) ,1=1,4) 

NA = APVSZ(3) +4. 

Hb = NA+1 

NC = 2,*APVSZ(3)+4, 

READ (5,4) (APVSZ(I) »I=b,NA) 

READ (5,4) {APVSZ(l) »I=Nb,NC) 

WRITE (6,21) (APVSZ(I) ,1=1,4) 

WHITE (b,22) ( APVSz ( 1 ) , 1=5 ,NA ) 

WRITE (o,22) (APVSZ(I),I=NU,NC) 

WRITE (6,23) 

READ (5,4) (UELVSZ(I) ,1=1,4) 

NA = DELVSZ ( 3 ) +4 , 

NB = NA+1 

NC = 2**DELVSZ( 3) +4. 

BEAU (5,4) (DELVSZ (I), 1=5, NA) 

READ (5,4) (DELVSZ ( I ) , I=NB , NC ) 

WRITE (6,21) (OELVSZ(I) ,1=1,4) 

WRITE (6,22) ( DELVSZ ( I ) , 1=5, NA ) 

WRITE (6,22) (UELVSZ(I) ,I=NB,NC) 

4 FORMAT (10E8.4) 

20 FORMAT (’O’ // S3X ’A, AP, OLLTa VS Z TABLES’ ) 

21 FORMAT (38X, 4L13«5) 

z2 format ax, ioti3.5) 

23 FORMAT (’O’) 

I HOT = 0 
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IF (PTI.GT.TI) IH0T=1 
CIH20Q) = 

CIH20SU) s 0, 

.CIW£QLiH..=L.Qa 

CIH20VU) s 0. 

IF. .{ PUNCH, EQ. 2* GR» PUNCH* EQ»3) 60 TO 90. 

C SET INITIAL CONDITIONS THROUGHOUT REACTOR BED 
DO. 50 Is2»NAXIAL 
CIH2U) = CH2I 
.C1Q2U) js. CQ22 
CIH20 ( I ) s CH20I 
C.IH20S(X) = CH20I 
CIH20UI) = CH20I 
CIH20V ( I ) = CH20I 
CIHE(I) = CHEI 

nil 5 t.x . 

PT< I ) = PTI 
bQ P ( I ) = PI 

C START REACTOR CALCULATIONS 

90 CALL MCALC (H2LEAD»DEBUG) 

91 CALL TMINS (TEND) 

. TOTELT s iend-tstart- 

IF (TOTELT »LT .0 0 ) T0TELT=3 440 ,-tSTART+TENU 
WRITE (6»100) TOTELT 

100 FORMAT (*0* /// 49X ******* OPERATIONS COMPLETE ******* / 
X 49X f TOTAL EXECUTION TIME! ♦ »F5.1 » » MINUTES 1 ) 

60 TO 6 
99 STOP 
END 
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c 

c 

c 


c 

c 

c 

c 

c 
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c 

c 

c 

c 

c 

C 

C 

L 

c 

c 

c 

c 

c 
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' SUBROUTINE' Me Ale ’ ( H2LEAD VDEBU G) 

ROUTINE TO’ HANDLE MAIN BODY OF REACTOR CALCULATIONS 

LOGICAL DEBUG »H2LEADrSSTATE . 

INTEGER PAS SNO* PUNCH r REG ION _ 

R~E AL MH2 » M02 1 MH20 , MHE p MUH2 , MU02 , MUHaO r MUHE , MU , MB AR » MFH2 r MFQ2 1 
X MFH20 f MFHE , KAPPA * KP , KCH2 , KC02 , KCH20 , MW , MU1VST t MU2VST » MU3VST , 

X MU4VST »KCF ,MBARSS 

PARAMETER NAX=70r NTI=4p 
PARAMETER NPP=30, NRTS=200 

PARAMETER NXI.0=KAX«10 . . __ 

COMMON /CBl/ TPP (NPP ,NAX ) »C02(NPP,NAX) ,CH2(NPP,NAX) ,CH20 (NPP»NAX) , 

X . T(NAX) rPCNAX) rPT(NAX) ,H(NAX) »CIH2(NAX) »CI02 (NAX) ,CIH2Q(NAX) , 

X CIH20S (NAX) , CIH20L (NAX ) ,CIH20V(NAX) ,CIHE(NAX) »OELT { NT I ) , 

X . PPT(NAX) 

COMMON /CB2/ Z (NAX) , 2SSL (NTI ) »ZSLL (NTI ) »ZLLV (NTI ) ,ZLW (NTI ) 

.. COMMON /CB3/ ALPHA t GAMMA , KAPPA.r HC f KP » A » QELTT t K.CH2 * KC02 1 KCH2Q-1 ICH2.*. 

X TC02»TCH20»PCH2fPC02,PCH20tMH2»M02,MH20»MHE,HR»HA,HAirHA2f 

X AW , CW » MW , DC » TA , D0H2 » DO 02 , D0H20 , ALPHAK , AP , DELTA » G * MU »HS , HL , 

X HV > TW , RHOM » WH2 , W02 , WH20S» WH20L » WH20V t WHE t CFBAR , MB AR » C ICP1 » 

X ClCP2rCICP3fRH0f RHOVfDELHCVTITPS>DTDZtDHDZ»TS»TL»DPDZ*DClDZt 

X DC2DZ,DC3SDZ,DC3LDZ,0C3VDZ,DC4DZ,RH0P,CP»DELTP»DPTPS,SDTIME, 

X PF t TSS* AC r.VC t PREST » GSS r MDARSS 

COMMON /CB4/ NOFZ*NPASS,NAXIAL , INPERP t MODPC , MODNO , PUNCH, ITSTEP , 

X JDP, ISHUT » IHOT 

COMMON /TABLES/ AVSZ (40 ) , APVSZ ( 40 ) , DELVSZ (40 ) ,HRVST ( 36 ) ,TC1<22) , 

X TC2 ( 20 ) , TC3V ( 20 ) , f C4 ( 18 ) f MUl'VST ( 34 ) , MII2VST ( 34 ) , MU3VST ( 32 ) » 

X MU4VST (34) ,CF1VST<34) ,CF2VST(34 ) ,CF3VST{34) ,CF3SVT(20) , 

X . CF4VSTC20) ,VPVST(faO) pTVSVP(4Q) »DHCVST.(24) 

DIMENSION MFH2 (NAX ) ,MF02(NAX) ,MFH20 (NAX } ,MFH£ (NAX) ,HW(NAX> , 

X - MINT (NAX) ,TM( NTI ) ,TWT (NTI ) ,GT (NTI ) ,PRESTT (NTI ) »PP (NX10 ) « 

X TT(NXIO) fTP(NXiO) ,FH2(NXI0) ,F02(NX10) ,FH20(NX10) ,FHE(NX10> 

DATA R, PI, OMEGA, KM /10.73, 3,1415927, 35800., 5 / 

DIF (UO,P»T) = 14,7*D0/P* ( T/492 . ) **1 ,823 

KCF (G,RHO,AMfUI,AP) = ,'6'1*G/RH0*( AM/RHO/DI ) **-.67*(G/AP/AM)**-.41 

*** i**************************************** ************************* 

PROGRAM FLAGS AND INDICATORS 

NOFZ I MO » OF INPUT AXIAL STATIONS 

NAXIAL: NOFZ+4 (ALLOWS FOR INTERFACE AXIAL STATIONS) 

■ * 

INPERP: NO, OF TIME INCREMENTS PER PASS THRU MESH 

MODPC : MODULAR PASS COUNT INDICATOR USED TO DETERMINE WHEN 

SUMMARY OUTPUT SHOULD BE CALLED 

’ MODNO : MODULAR NUMBER OF INDEPENDENT VARIABLES USED FOR 

’ PRINTER PLOT OUTPUT 

IOCT : PROGRAM COUNT OF THE NO, OF TIMES SUMMARY I/O HAS 

' • BEEN CALLED 

IT J PROGRAM COUNT OF CURRENT TIME STEP INCREMENT 

iz : program count of current axial station 
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PUNCH AUXILIARY I N PUT / OUTPUT INDICATOR j 

=1 NO AUXILIARY I/O 

=2_R£AQ. AUXIUABY._1NEUT.,QNLX 

=5 PUNCH AUXILIARY OUTPUT ONLY 
, =4 READ AUXILIARY .INPUT* AND PUNCH AUXILIARY GUIPUT_ 

IT.STEP; TIME INCREMENT .INDICATOR . . 

=0 PROGRAM DETERMINES DELT INCREMENTS 
. =J OEL.T INCREMENTS. ARE.. INPUT. 

IHOT J HOT START INDICATOR — — 

=0 NORMAL START 
=1 HOT START 

ISHUT ; REACTOR SHUTDOWN INDICATOR — * 

=0 NORMAL REACTOR OPERATION 
=X REACTOR SHUTDOWN HAS OCCURRED 

JDP : CONST AMT DELT INDICATOR (USED ONLY WHEN ITST£P=0) 

C =0 PROGRAM DETERMINES DELT AS USUAL 

C =-l DELT SET CONSTANT AT .001 SEC 

C 

C MINT3 ; DIFFUSION CONTROL INDICATOR 

C =0 REACTOR CALCULATIONS NOT YET DIFFUSION CONTROLLED 

C =1 REACTOR CALCULATIONS ALL DIFFUSION CONTROLLED 

C 

C **#.***************************************************************** 

C INITIALIZE FLAGS AND COUNTERS 
ID CT = 0 
ITP = 0 
MINTS = 0 
JDP = 0 
ISHUT = 0 
IT = 1 
IZ = I 
K - 1 
Ni = 1 
PASSNO = 1 

n = i 

MOD PC =10 

IF INPASS.LT.IO) MOuPC=NPASS 
DC TEMP = (29160. -6. *P Cl) )/57, 

SSTATE = .FALSE. 

ZEND = Z(NOFZ) 

C INITIALLY ALL INTERFACES COINCIDE WITH END OF REACTOR BED 
DO 7 1=1, NTI 
ZSSL(I) = ZEND 
ZSLL(I) = ZEND 
ZLLV(I) = ZEND 
7 ZLVVCI) = ZEND 

0.0 TO (102*102*102*111)* PUNCH 
102 -ALPHA = 0. 

IF (IHOT.EQ.l) 00 TO 33 
IF (TCU-DCTEMP) 25*25,24 

24 SSTATE = .TRUE, 

60 TO 111 

25 DO 32 I=1,NAXIAL 

32 MIN T( I ) = 1 
GO TO 35 

33 DO 34 1=1 *NAXIAL 

34 MIN'I(I) = 2 
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60 TO 35 

111 DO 112 I=1»NAXIAL 

112 MINT ( X ) = 3 

C -DETERMINE APPROPRIA.TE DE].JT 

35 "IF" {iTsTEP.EQ.1) GO TO 19 

TPSMAX = PT (2) 

CALL DELTAT 
GO TO .17 

19 DELTT = DELT(l) 

17 PTUi = PT<2>„ _ 

C SET INITIAL PARTICLE SURFACE VALUES FOR TEMPERATURE, CONCENTRATION 
TPS = PT (1 ) 

TW s TPS 
CSH2 = CIH2 (2) 

CS02 = CI02(2) 

CSH2Q = CIH20(2) 

IF IH2LEAD) GO TO 121 
G = GSS*SQRT ( (PF-PREST)/ (PF-P(l) ) ) 

CI02(1) = CI02 { 1 ) *PREST/P { 1 ) 

CIH2 ( 1 ) s CIH2U)*PREST/P(1) 

CIHE(1) s CIHEU)*PREST/P(1) 

121 TIME = ,0. 

IF (PUNCH. NE. 4)- GO TO 113 
CALL PREVUS 

113 IF (PUNCH, EG. 2, OR. PUNCH. EG, 3) CALL RDPNCH{TIME,HW,MINT) 

CMSUM = CIH2 ( 1 ) /MH2+CI02 ( 1 )/M02+ClH20 { 1 ) /MH20+CIHE (1 ) /MHE 
MF02 ( 1 ) = CI02 { 1 ) /M02/CM5UM 

REGION = .1 

CALL UNBAR (TVS VP * 1 ,P (1 ) , 0 » ,TL , KK ) 

IF (TU)-TS) lfl»0 
B IF (T(l)-TL) 9,9, 1U 

9 REGION s 3 
GO 10 1 

10 REGION = 5 
t 

C BEGIN STEPPING AXIALLY THROUGH THE REACTOR 
C 

1 12 = 12+1 

IRHOM = 0 

C CHECK IF CALCULATIONS DIFFUSION CONTROLLED AT CURRENT AXIAL STATION 
M2 = MINT (12) 

GO 10 (37,39,47) , M2 

39 IF i P ASSNO , EG, 1, AND. IT, EG, 1) GO TO 3c 
IF (PPT(IZ>-PT(iZ) ) 37,37,36 

36 SSTATE = .FALSE, 

GO TO 36 

37 IF (PT(IZ)-DCTEtMP) 31,30,30 

30 MINT ( IZ) = 3 

GO TO 47 

C NOT DIFFUSION CONTROLLED 

31 sstate = .false, 

MINT(IZ) = l 
GO TO 3ri 

C DIFFUSION CONTROLLED 

47 SSTATE = .TRUE. 

IF (PUNCH, NE. 4) REGION=5 
36 UELT2P = DELJZ 

DELTZ = Z(I2)“4(U"1) 

TAV = (PT (12-1 ) +T( 12-1 ) ) /2 . 

CALL UNBAR ( AVSZ , 1,2 ( 12-1 ) , 0 , , A .KK) 

CALL UNBAR ( APVSZ ,1 tl ( 12-1 ) ,0 . , AP»Kh 1 
CALL UNBAR (DELVSZ r 1,2 ( 12-1) ,0. , DELTA, KK) 
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CALL UNBAR (HRVSTr IfTdZ-l) pO. pHRpKK) 

CALL UNBAR (DHCVSTp 1 »T aZ-1 > ». 0.i fDELHC P KK ) 

CALL UNBAR <MUlVSTplfTAVf0.pMUH2pKK) 

CALL UNBAR (MU3VSTp1pTAVp 0. f MUH20»KK) 

CALL . UNBAR {MU4VST TAV » Q.*j.MUH£j.KK.) . 

CALL UNBAR (CFlVSTp 1 p T( IZ-1 ) 1 0 . #CFH2pKK ) 

CALL UNBAR LCF2VSTi.lL r T { IZ-1) t ft. tCFOZ t KK 1 
CALL UNBAR (CFAVSTplrT(IZ-l) pO. pCFHEfKK) 

CALL UMBAR...XT.Cli If .T.AV.fO . 

CALL UNBAR (TC2p1pTAVpO. pTH02pKK) 

CALL UNBAR (TC3Vp Ip TAV p 0. r THH20pKK) 

CALL UNBAR (TC4 p 1, TAVpO. pTHHEp KK) 

00 TO (2i2r3rAr5)f REGION 

2 CALL UNBAR CCF3SVT,1pT( IZ-1) p 0, pCFH20Sp KK) 

CFH2QL s 1«.Q 

CFH20V = 0, 

RHOM =0, 

GO TO 6 

3 CFH205 = 0. 

CFH20L =1.0 
CFH20V = Q. 

RHOM = 0. 

GO TO 6 

4 CALL UNOAR (CF3VST , 1 ,T ( IZ-1 ) , 0 , , CFH20V fKK ) 

CFH20S = 0» 

CFH20L =1.0 
RhOM = ft. 

GO TO 6 

5 CALL UNBAR (CF3V5T» 1 , T( IZ-1) , 0. »CFH20V»KK) 

IF (IZ.EQ.2) GO TO m 

IF (H(IZ-l) ,L.T.H{IZ-2) ) GO TO 83 
SA IF (T(IZ-l)-TS) 81 f 81 f 82 
81 CALL UNOAR (CF3SVTp 1 1 T( IZ-1) 1 0 . »CF’H20Vi KK) 

GO TO 83 

62 IF CT(IZ-l) .LT.TL) CFH20V=1.C 
83 CFH205 = 0. 

CFH20L =0. 

RHOM = ALPHA*CiH2(IZ-l)*Cl02(IZ-i)*EXP(-0MEGA./T(IZ-l) ) 

6 RHO = CIH2( IZ-l]+Cl02{ IZrl 1+CIH2QSI IZ~1)+CIH2QL( IZ-D+CIH2QV ( IZ-1 3 

X +CIHEUZ-1) 

RHOS = RH0-CIH20L( IZ-1 )-ClH20V( IZ-1) 

RHOL = RHQ-CIH20S ( IZ-1 )-CIH20V ( IZ-1 > 

RHOV = KH0-CIH20S { I Z-l) -CIH20L { IZ-1 ) 

SUM = C IH2( 12-1) /MH2+CI02( IZ-1) /M02+ ( CIH2QS < IZ-1) +CIH20LC IZ-1 )+ 

X CIH20V { IZ-1 ))/MH20+CIHE( IZ-1 J/MHE 

MU = (MUH2*CIH2 ( IZ-1 ) /MM2+MU02+CI02 ( IZ-1 ) /M02+MUH20* (CIH20S ( IZ-1 )+ 
X C1H20L ( IZ-1 ) +CIH20V { IZ-1 ) ) /MH20+MUHE+CIHE ( IZ-1 )/MHE) / SUM 

MBAR = RHOV/ { CIH2 ( IZ-1 ) /MH2+CI02 ( IZ-l ) /M02+C IH20V ( IZ-1 ) /MH20+ 

X CIHEUZ-D/MHE) 

CFBAR = <CFH2*CIH2( IZ-1 ) +CF02*CI02 (IZ-1 ) +CFH205*CIH2QS { IZ-1) + 

X CFH20L*CIH20LUZ-1)+CFH20V*CIH20VU2-1)+CFHE*CIHE{IZ-1)} / 

X RHO 

DIH2 = DIF (D0H2 1 P ( IZ-1 ) p T ( IZ-1 ) ) 

DI02 = DIF{D002pP(IZ- 1) rT(lZ-l) ) 

0IH2O = DIF (DOH20 pP (IZ-1 ) pT (IZ-1 ) ) 

KCH2 = KCF(GpRHOpMUpDIH2»AP) 

KC02 = KCF(GpRH0»MUpDI02pAP) 

KCH20 = KCF(G»RH0»MU»DIh20»AP) 

SUMC = C IH2 ( IZ-1 ) /MH2+CI02 ( IZ-1) /MQ2+CIH20 (IZ-1 ) /MH20+CIHE ( IZ-1) / 

X MHE 

THERMC = (THH2*CIH2 ( IZ-i) /MH2+TH02*CI02( 12-1 ) /M02+THH20* 
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X CIH2Q.UZ-1 ) /MH204-THHE*CIHE { IZ-1 ) /MHE> /SUMC 

HC = ,61*CFBAR*G*(G/(AP*MU))**-,41*(CFBAR*MU/THERMC)**-,667 
HWCIZ) = HC 

H W < U a HW (2 )_ „ 

58 " IF (SSTATE) GO TO 60 

C PARTICLE ANALYSIS TRANSIENT 

CALL PC ALCT ( CSH2 r CS02 » CSH20 1 TPS , DCTEMP t IT t ITP » IZ , PASSNO * NC ♦ DEBUG ) 
ITP = IT 

IF (NC,EG,0) GO TO 5S 

. T?SJ = iPTlIZJ -PPJ ( IZJL> /pELJPj?,lDELJp+DELTT. )L±£EI LIZ) 

IF (TPST.LT, DCTEMP) TPST=DCTEMP 
TPS = TPST 

CS02 s CI02(IZ)-HC/(HR*KC02)*(T(IZ)~TPS) 

CSH2 = CIH2(IZ)-2.*KC02/KCH2+MH2/M02*(CI02(IZ)-CS02) 

CSH20 = 0, 

GO TO 59 

C PARTICLE ANALYSIS — STEADY SfAfE 

60 CALL PCALCS (TPS r CS02 1 CSH2 » CSH2Q t C I0£PR » IZ # M I NT ) 

59 TITPS = T ( IZ-1) -TPS 
CICP1 = C.IH2 ( IZ-1 ) -CSH2 
CXCP2 = CI02 ( IZ-1) -CS02 
C.ICP3 CIH20V ( IZ-1 ) -CSH20 
IF (CSH20) 85»85?86 

85 CICP3 = 0. 

86 WH2 = CIH2(IZ-1)/RH0 
W02 s C 102 ( IZ-1 >/RHO 
WH20S = CIH20S(IZ-1)/RHO 
WH20L = CIH20LUZ-U/RH0 
WH20V = CIH20V ( IZ-1 ) /RHO 
IF (REGION, NE. 5) GO TO 87 
WH2US = 0. 

WH20L = 0. 

87 WHE s CIHEUZ-U/RHO 

C make derivative calculations 

900 CALL DERI V (ClOZPRrDELTZ r IZ , REGION ) 

C 102 ( IZ ) = CI02(IZ-1)+DC2DZ*UELTZ 

IF (CI02(IZ).LT.O..AND.RHOM.GT.Q,) GO TO 901 
IF (iRHOM.EO.i) CIO2(IZ)=0, 

00 TO 903 

901 RHOM = 0. 

I WHOM s 1 
GO TO 900 

903 P(IZ) = P ( IZ— 1 ) +DPDZ*DELTZ 
C TEST FOR NEGATIVE PRESSURE 
IF (P(IZ)) 905fbl»61 

61 TilZ) = T(IZ-1)+DT0Z*PELTZ 

IF (T(IZ) »LT ,T(1) ) T(IZ)=T(1) 

PPT(IZ) = PT (IZ) 

PT ( IZ) = TPS 
PTU) a PT (2) 

H ( IZ) = H(IZ-1)+DHDZ*DELTZ 
IF (H(IZ).LT.H(lZ-l)) REGI0N=5 
IF (H(IZ) .LT.O.) H(1Z)=0. 

CIH2(IZ) = CIH2(IZ-1)+DC1DZ*DELTZ 
CIH20S(IZ) = ClH20S(IZ-l)+DC5SDZ*DELTZ 
CIHaoL(IZ) = CIH20L (IZ-1 ) +DC3LDZ*DELTZ 
CIH20V ( IZ) = CIH20V ( IZ-1 ) +DC3VDZ*0ELTZ 
IF (REGION, NE. 5) GO TO 88 
CIH20S(IZ) = 0, 

CIH20L(IZ) a 0, 

88 CIHt(IZ) = CIME( 2Z-1 ) +DC4DZ*DELTZ 
IF (CIH2 ( IZ) ,LT ,0 • ) CIH2(IZ)=0, 
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IF (CI02( IZ) .LT.O. ) CI02UZ)=0, 

IF -1CIH2QV L1Z1 «.LT., Q » ) CIH20y.(.IZ)sD. . 

IF (CIHE(IZ) .LT.O.) CIHEUZ)=0, 

~ CALL UNBAR (TVSVP»1#P(IZ) »0.»TL»KK) 

. .CMSUM S ClH2aZl/MH2+Cl02ClZi/MQ2+CIH2Q(IZ>y«H20+CIHE:(lZ)/MHE 
MF02(IZ) = CI02(IZ)/M02/CMSUM 

IF. UMFP2,UZ)..kT . «Q1*MFQ2 ULLL.AND. (PmZ)_.gT...(T (1) +30. J )J. MIHTXIZ). 
X < c 3 

flfl. S1.X=2 jNAX-IAL 
IF (MINT ( I)-3) 140»91,140 
91 CONTINUE 

BEGIN TESTING FOR INTERFACE LOCATIONS 

GQ XO .Ufli.0.f.aQrI2ai.I3.0,i40)f REGION .. 

C SOLID/SOLID-LIQUID INTERFACE 

100 IF (T(IZ)-TS) 5l,i01»101 

101 CALL OVSTEP (ZEND»IZ,PASSNO, REGION) 

GO TO 140 

C SOLID-LIQUID/LIOUID INTERFACE 
110.. IF. (CIH2QS11Z) ) lOir.lQliSl 
C LIQUIU/LIQUID-VAPOR INTERFACE 
120 IF (T(IZ)-TL) Sip lOlrlOl 
C LIQUID-VAPOR/VAPOR INTERFACE 
130 IF (CIH20L ( IZ) ) lOlrlOlpSl 
C End OF REACTOR BED 
140 IF (ZCXZ)-ZEND) 51r70#lQl 
Si IF (IZ.LT.NOFZ) GO TO 1 
70 NM1 s NOFZ-1 

DTWOT = (PI*DC*HSUM(NMI)-AW*(HA*(TW-TAHHA1’MTW-TA)**1.£5+HA?+ 

X (TW**4-TA**4) ) )/(MW*C«i> 

T w = TW+DTWUT*UELTT 
C CALCULATE MOLE FRACTIONS 
54 DO S3 I=1 pNAXIAL 

C.MSUM = CIH2(I)/MH2+CI02(I)/M02+CIH20(I)/MH20+CIHE(I)/MHE 
MFH2 ( I ) = CIH2U)/MH2/CMSUM 
MF02(I) '= CI02{I)/M02/CMSUM 
MFH*iO(I) = CIH20 ( I ) /MH20/CMSUM 
53 MFHE(I) = CI.HEU}/MHE/CMSUM 

IF (H2LEAD.OR.TIME.GT.S0TIME) Go TO Ufa 

QPDT = R*T ( IZ > /MBAR+AC/VC* ( G-PREST*GSS*SQRT ( HBAR*TSS/ ( T ( 12 ) *.MBARSS 
X ) )/P( 1 ) ) 

PPP r R*AC*GSS/ ( VC#P(1 >*SQRT (MBARSS/TSS) ) *SuRT (T( IZ) /MBAR) 

DPDT = OPDT* ( 1,-EXP C-PPP*DELTT) J/(PPP*DELTT) 

PKEST = PREST+DPDT*QELTT 
6 = GSS*SQRT( (PF-PREST)/(PF-P(1) ) ) 

C102(l) = CI02 ( 1) *PREST/P ( 1 } 

CIH£(1) = CIH2 ( 1} +PRE.ST/P { 1 ) 

CIHL(1) = CIHE(1)*PREST/P(1) 

C END OF REACTOR HAS BEEN REACHED FOR CURRENT TIME STEP 
116 TIME = TIME+DELTT 
C CHECK FOR DEBUG PRINT 

IF (DEBUG) CALL DEBUG3 (TIMEp IT ) 

IF (PUNCH. EG. 4) GO TO 104 
DO 107 I=2pNAXIAL 
IF (MINT ( I )-3) l04p 107t 104 
107 CONTINUE 

-CALL RPNCH4 (TIME) 

MINT3 = 1 
NPASS = PASSNO 
GO TO 108- 
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104 IF (ISHUT.EQ.l) GO TO 42 
C TEST IF SHUTDOWN TIME HAS BEEN REACHED 

IF (TIME-SDTIME) 42,42,40 
C SHUTDOWN _ _ 

40 I SHUT = 1 

READ (5,41) CIH2(1> ,CI02<1> 

41 FORMAT (2E10.5) 

42 IT = IT+1 
IZ = 1 

REGION s x 

IF (T(l)-TS) 78,78,76 
76 IF (T(D-TL) 74,74,75 

74 REGION = 3 
GO TO 7B 

75 REGION = 5 

76 JF (PUNCH,NE,4) GO TO 89 
CALL DELETE (PASSNO) 

NOFZ = NAXIAL-4 

C DETERMINE APPROPRIATE DELTT FOR NeXT TIME STEP 
69 DELTP = DELTT 

TPSMAX a FMAX(PT,NAXIAL,NAX,1) 

CALL DELTA! 

IF (ITSTEP.NE.I.OR.JDP.EG.I) GO TO 26 
IF (TPSMAX-38Q. ) 26,23,23 
23 JDP = 1 

DELTT = .001 
DpTPS = TPSMAX 

2d IF IIT.LE.INPEKP) GO TO 1 

C COMPLETION OF A PrtSS CHECK OUTPUT REQUIREMENTS ■ 

108 IT - 1 

M 2 s PASSNO 

MODPC = NO. OF COLUMNS OF PRINT ON SUMMARY OUTPUT PAGES — MUST BE 

RESET AS PASS NUMBER COUNT (PASSNO) EXCEEDS 10 MAX. NO. OF OUTPUT 

COLUMNS PER PAGE IS 10 IOCT IS COUNT OF HOW MANY TIMES SUMMARY 

OUTPUT HAS BEEN PRINTED MAX. OF 3 TIMES (IOCT = 0,1,2) 

IF (IOCT.EQ.Q) GO TO 55 
GO TO (56,57), IOCT 
bo MODPC = NPASS-10 
GO TO 55 

57 MODPC = NPAS3-20 

55 TWT (PASSNO) = TW 

TM( PASSNO) = TIME 

IF IH2LEAD.0R.TIME.GT .SOT I ME) GO TO 28 
GT (PASSNO) = G 
PRESTT( PASSNO) = PREST 

C PASS COMPLETE SET Up TEMP, PRESSURE, MOLE FRACTION CURVES FOR 

C PRINTER PLOTS 

28 CALL SETUP (TM,HFH2,MF02,MFH20,MFHE, PASSNO, K, KM) 

C STORE TEMP, PRESSURE, MOLE FRACTION PROFILES FOR LATER PRINTOUT 
DO 62 1=1, NOFZ 
MKT = M+MODPC*(I-l) 

TTIMKT) = Til) 

TP (MKT) = PT('r) 

PP (MKT) s PCI) 

FH2 (MKT ) = MFH2 ( I ) 

F02CMKT) = MF02 ( I ) 

FH20 ( MKT ) = MFH80 ( I > 

62 FHEIMKT) = MFHE(I) 

M = M+l 

IF (ITSTEP.EQ.O) GO TO 27 
DELTT = DELT (PASSNO) 

C CALL OUTPUT EVERY 10TH PASS AND AETER LAST PASS 
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27 IF (MOD(PASSNO,10).EG.Q,OR.PA5SNO.EG.NPASS) GO TO 63 
RASSNO = PA_SS.NO.tl 
C RETURN TO INLET 

fifi .JDQ.JL 

63 NNAX = NAXIAL 
NAXIAL = NOFZ - 

CALL OUTPUT (TM,TT»TP»PP»FH2»F02,FH20,FHE,N1»N2,IOCT) 

CALL TM.PRNT (TM,.TWT.,N1.,N2J .. ... 

IF (.N0T.H2LEAD, AND. TIME, LT.SDTIME) CALL TMPRP (TM » GT » PRESTT » N1 » N2 ) 

NAXIAL. s. NNAX 

M = 1 

N1 = N2+1 

PASSNO = PASSNOtl 

IF (MINT3.EQ.Q), GO TO 73 

WRITE (6,79) TIME 

79 FORMAT 4.6X 3.8 (.**») / 47X ** GQ TQ ALL D IFFUS ION-jCQN TROLLED RU 

XN ** / 47X ♦* TIME =*,E9.4,* SEC*, ?X*** / 47X, 38(**») 

X ) 

IF ( I SHUT »EQ, 1 ) GO TO 72 

C DUMMY READ STATEMENT POSITIONS PAST FINAL INPUT CARD IF 5HUTS0WN TIME 
C HA5 NOT BEEN REACHED 

71 READ (5,41) DUMMY 1,.UUMMY2 

C RETURN TO MAIN PROGRAM IF CASE HAS REACHED ALL DIFFUSION CONTROLLED 

C STATUS (MINT3=1) EVEN IF AMIDST A PASS 

72 RETURN 

73 IF (PASSNO. GT.NPASS) GO TO 77 
C RETURN TO INLET 

GO TO 1 

77 IF (PUNCH. EG.l. OR. PUNCH. EG. 3) CALL RPNCH3 (TIME, HW, MINT) 

IF (ISHUT.EU.O) GO TO 71 
RETURN 

ERROR MESSAGES FOLLOW 

905 WRITE (6,902) P ( XZ ) , UPUZ , IT , PASSNO , IZ , REGION 

902 FORMAT (*0* 5X , 120(***) // b6X ’ERROR ERROR* / 53X ’PRESSURE 

X HAS GONE NEGATIVE* / 48X 'P =*,£11.5, 5X *DPDZ =*,E11.S / 

X 52X ’TIME INCREMENT* ,13, * OF PASS*, 13 / 

X SIX 'AXIAL STATION*, 14, * ... REGION', 12 // 5X»120('#«) ) 

CALu EXIT 


INTERNAL FUNCTION AND SUBROUTINE SUB-PROGRAMS FOLLOW 


FUNCTION RSUM (NMD 
R5UM z 0. 

DO 101 1=1, NMl 

101 RSUM z RSUM+,5*(HW (1+1)*(T( 1+1 )-TW)+WW ( I ) * (T ( I )-Tw ) )*(2(I+1)-Z(D) 

RETURN 

1 


SUBROUTINE PKEVUS 

ROUTINE TO SET UP AND INTERPOLATE TEMP, PRESSURE, CONCENTRATION VS Z 

DIMENSION TABLE1 (64) ,TABLE2(64) ,TABLE3(84) ,TABLE4(a4) ,ZT(40) , 

X PTT(40) ,C02T(40) ,CH2T(40) ,CH20T(40) ,TABLE5(84) ,TTI(40) 

DATA (TABLED!) ,1=1,4) / 0., 1., 0., 0. / 

DATA ( TABLE2 ( I ) , 1=1 , 4 ) / 0., 1., 0., 0, / 

DATA (TABLL3( I) ,1=1,4) / 0,r 1., 0., 0, / 
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DATA (TABLE4(I).,I=1,4) / o.# i., 6., 6. / 

DATA (TABLE5U) , 1=1,4) / 0,* i.r 0_ e r Q . yf 

CAUL RPNCH2 (TIME, NAXT, ZT,PTT»C02T,CH2T,CH20T,TTI) 

TABUE1< 3) a NAXT 

TABLE2(3) = NAXT 

TABLE3(3) a NAXT 

T A0LE4 { 3 ) a NAXT 

TABLE5(3> = NAXT 

DO 1 1=1, NAXT 
TABLE 1U +4) = ZXLI.T 
TABLES (X+4) = ZT(I) 

TABLE3(I+4) = ZT(I) 

TABLE4 ( 1+4) = ZT { I ) 

X TABLES (144) a ZT(.I) 

DO 2 1=1 , NAXT 

TABLE! t T+4+NAX.T.) = . P.T,T.(.I ) 

TABLES ( I 444NAXT) = C02T(I) 

TABLE3( I+4+NAXT) = CH2TII) 

TABLE4(l444NAXT) = CH20T(Ij 

2 TABLE5(I44+NAXT) = TTI(I) 

DO 3 I=1,NAXIAL 

CALL UNBAR (TABLE!, if 2(1.) ,Q..,PTCI) ,KK) 

CALL UNBAR (TABLE2 pi ,Z{X),0.,CIQ2(I)»KK) 

CALL .UNBAR (TABLE3 1 X »Z ( I ) » 0 . » CIH2 (I ) ,KK) 

CALL UNBAR (TABLE4 , 1 , Z ( I ) , 0 . , CIH20 ( I ) » KK) 

call unbar (tables, i,zu) »o.,tci) ,kk> 

3 CIHLOV(I) = 0. 

return 


subroutine oeltat 

routine to DETERMINE APPROPRIATE uELTA T INCREMENT BA5ED ON • TEMPERATURE 

and pressure 



IF (ITSTEP.EQ.l) GO TO 9 
IF (57« *TPSMAX+12, *P ( 1 ) -17280 . ) 

10,10,1 

1 

IF (S7.*TPSMAX+12.*P(1)-1999G.) 

11,11,2 

2 

IF ( 100. *TPSMAX+21.*P(1) -38000. 

) 12,12,3 

3 

IF (TPSMAX-380 . ) 13,13,14 


14 

OELTT = .001 


13 

JDP = 1 

DPTPS = TPSMAX 
GO TO 9 
DELTT = ,002 


12 

GO TO 9 
DELTT = .01 


11 

GO TO 9 
DELTT = .05 


10 

GO TO 9 
DELTT = .10 


9 

RETURN 



END 
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SUBROUTINE BERTV ~\Cl OgPRTEELTZ , I Z . REG I ON ) " 

C ROUTINE TO CALCULATE DERIVATIVE EXPRESSIONS IN EACH REGION 
C 

INTEGER REGION 

REAL KCH2 p KC02 p KCH20 p MHg t M02 * Mh|20 1 MHE » MU p MB AR 

PARAMETER NAX=7o f NT 1=40 
PARAMETER NPP=30 

COMMON /CB1/ TPP(NPPpNAX) »C02 (NPPpNAX) pCH2(NPPpNAX) pCH20 (NPPpNAX) * 
X T(NAX) ,P(NAX) fPT(NAX) ,H(N aX) ,CIH 2(NAX) pCI02(NAX) pClH20(NAX) p 

X CIH20S (NAX) pCIH20L(NAX) pCIH20V(NAX) pCIHE(NAX) #DELT (NTI ) t 

X PPTJNAXl 

COMMON /CB2/ Z (NAX) pZSSL (NTI > pZSLL (NTI ) pZLLV (NTI ) ,ZLVV (NTI ) 

COMMON /CB3/ ALPHA p GAMMA pKAPPA t HC pKPpApDELTT pKCH2pKC02p KCH20pTCH2p 
X TC02 p TCH20 p PCH2 p PC 02 p PCH20 »MH2 1 M02 , MH20 p MHE . HR p HA , HA1 p HA2 p 

X AW p CW p MW p DC p T A p D0H2 p DO 02 p D0H2O p ALPHAK p AP t DELT A » G t MU » HS » HL » 

X HV p TW pRHOMp WH2 p W02f WH20S» WH20Lp WH20V p WHEp CFBAR t MBARpClCPl p 

X C I.CP2 p C I.CP3.P RHO p RHOV p DELHC p IMPS p.DTDZ i D.HD2 p.TS i.TL i B£DZ*Q£li2Z*. 

X . DC2DZ pOC3SDZp 0C3LQZpDC3VD2pDC4DZpRH0PpCPp DELTP * DPTPS pSDTlMEp 
X . PF pTSSp AC p VC pPRESTpGSSpMBARSS 

COMMON /CB4/ NOFZ t NPASS p NAXI AL p INPERP p MODPC t MODNO * PUNCH p I TSTEP p 
X JDPpISHUTpIHOT - 

UATA DELHF/143.4/ 

IF (IZ.EQtNOpZ) GO TO 10 
QZPK = Z(IZ+1)-Z(IZ-L) 

.10 FACT1 = (l.-EXP(-KCH2*AP*RHQV*0ELTZ/G> )/(KCH2*AP*RH0V*QELTZ/G> 
FACT2 s (l.-EXP(-KC 02 *AP*RH 0 V*DELTZ/G) )/(KC02*AP*RH0V*DELTZ/G) 
FACT3 = (l.-EXP<-KCH20*AP*RH0V*DELTZ/G> )/{KCH2O*AP*RH0V*DELTZ/G) 
FACT2P = (l.-EXP(-KC02*AP*RH0V*DZPR/G) )/<KC02*AP*RH0V*DZPR/G) 

QwlOZ = {.-(RH0M*DELTA+KCH?,*AP*CICP1)/G)*FACT1 
Dw2l)Z = <-(RH0M+OELTA*.5*MO2/MH2+KCO2*AP*CICP2)/G)*FACT2 
UH3DZ = ( <RH0M*DELTA*MH2Q/MH2-KCH20*AP*CICP3)/G)*FACT3 
DW20ZP = DW2DZ+FACT2P/FACT2 . 

UW4UZ = 0. 

DRDZ = (DELTA-1 . ) / ( 144 . *DELTA**3) * ( 1 , 75+75 . *MU* ( 1 . -DELTA ) / ( A*G ) * 

X (G*G/(B4.4*A*KH0}}) 

GO TO (1p2p3p4p5)p REGION 
C ICE REGION 

1 HRD = 0, 

DELTAH = DELHC+DELHF 
DW3DZ = DW3DZ/FACT3 
QWDZ3 = 0. 

UW3SDZ = DW3Q2 
DW3LDZ = 0, 

DW3VDZ = 0, 

W3 - 0* 

HFACT = 0» 

RHOi = RHO/ ( 1 , -WH20S ) 

DteJUZ = DW30Z 
RHO<; = RHOV/KHO 
GO TO 6 

C ICE-LIQUID REGION 

2 hKO = 0. 

HSL = H(IZ-i)+UELHF*WH20S 

DELTAH = DELHC+DELHF* (HSL-H(IZ-i) )/{ hSL-HSJ 
Dw3L)Z = DW3DZ/FACT3 
Dpv'023 = 0, 

DW3SDZ = DW3DZ*( (HSL-H(IZ-l)}/(HSL-iiS))-OHDZ/(HSL-HS) 
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SUBROUTINE PCALCT <CSH2,CS02#CSH20#TPS#DCTEMPiIT»XTP»IZ»PASSN0# 

X,. NCQN# DEBUG) .... .. 

ROUTINE TO. HO PARTICLE ANALYSIS. QFTEMP# CONCENTRATION DISTRIBU TIONS . 

AT .PRE-STEADY STATE (TRANSIENT) CONDITIONS 

LOGICAL DEBUG"'" 

INTEGER PASSNOr PUNCH 

REAL LAMDA # LAMDAO , KAPPA # KP , KC02 #KCH2 # KCH20 # M02 # MH2 # MH20 

PARAMETER NAX=7G# NTI=4Q. 

PARAMETER NPP=30» NRTS=200 

COMMON /CBi/ .TP.P 1NP.P rNAXl.tX.QE.lNPJ.tNAX.l r.(^INPP.r!iAXJ.»J^i2QlNPE<MAX)^. 

X T(NAX) »P<NAX) »PT (NAX) #H(NAX) #CIH2(NAX) »CI02(NAX) #CIH20(NAX> » 

X CIH20S (NAX) t CIH20 L(NAX) iCIH2QV (NAX) *CIHE (NAX) #DELT (NTI) » 

X- PPT(NAX) 

COMMON /CB3/ ALPHA t GAMMA # KAPP A , HC » KP » A rOELTT # KCH2 »KCQ2 » KCH20 # TCH2 # 

X TC02 , TCH20 , PCH2 # PC 02 » PCH20 , MH2 # M02 , MH20 # MHE # HR # HA # HA1 »H A2 # 

X .. AWfCW»MWia.C»J.AfDQH2»Di)02»DOH20^ALEHAK*AP»DELIA*S4.MU*HSjddL J 

X HV , TW , RHOM , WH2 , W02 , WH20S , WH20L » WH20V , WHE » CF8AR , MBAR , C I CPI # 

X C I CP2 r C I CP3 » RHO , RHOV » DELHC » Tl TPS » DTDZ * DHDZ t TS # TL # DPDZ . DC1DZ # 

X DC2DZ » DC3SDZ # DC3LDZ » DC 3VDZ # DC4DZ# RHOP , CP # DELTP # DPTPS , SDT I ME » 

X PF » TSS » AC # VC » PREST t GSS t MBARSS 

COMMON /CB4/ NOFZ»NPASS»NAXIAL»lNPERP»MODPC r MODNO»PUNCH*ITSTEPf 
X JDPiISHU.T #.IHOT 

DIMENSION GRT(NPP»NPP) #GRTI (NPP,NPP) » GRTT I ( NPP » NPP ) ,GR02(NPP»NPP ) , 

X GR02I (NPP, NPP) »GR02TI (NPP# NPP) »GRH2(NPP,MPP) ,GRH2I (NPP# NPP) »' 

X GRH2TI(NPP#NPP) # TP (NPP #2) #CP02(NPP#2) »CPH2(NPP»2) #CPH20(NPP# 

X 2) » RHET ( NPP » 2 ) # LAMDA (NPP ) f LAMDAO (NPP) » RATE (NPP) #RTT (NRTS) » 

X RT02 (NRTS) #RTH2 (NRTS) #F(NPP) 

DATA PI/3. 1415927/ 1 RAOP/ .656E-7/ 1 TOL/. 0001/# MAX/25/ 

PARTICLE ANALYSIS NOMLNCLATURE FOR 02 + 2H2 — 2H20 REACTION 


CS02 J PARTICLE SURFACE CONCENTRATION OF 02 

CSH2 : PARTICLE SURFACE CONCENTRATION OF H2 

CSH2Q J. PARTICLE SURFACE CONCENTRATION OF H2Q 

TS : PARTICLE SURFACE TEMPERATURE 

CP02(IrX) J PARTICLE CONCENTRATION DISTRIBUTION OF 02 AT TIME 1 
CPH2 ( I » 1 > : PARTICLE CONCENTRATION DISTRIBUTION OF H2 AT TIME 1 


CPH2Q ( I # 1 ) 
TP ( I » l ) 
CP02 ( I #2) 


: PARTICLE CONCENTRATION DISTRIBUTION OF H20 AT TIME 1 
: PARTICLE TEMPERATURE DISTRIBUTION AT TIME I 
J PARTICLE CONCENTRATION DISTRIBUTION OF. 0.2.AT XI ME 2. 


CP02A(D 


(GUESSED) 

PARTICLE CONCENTRATION DISTRIBUTION OF 02 AT TIME 2 


CPH2 ( I r 2) 
CPH20(I#2) 
TP ( I #2) 


(CALCULATED) 

PARTICLE CONCENTRATION DISTRIBUTION OF H2 AT TIME 2 
PARTICLE CONCENTRATION DISTRIBUTION OF H20 AT TIME 2 
PARTICLE TEMPERATURE DISTRIBUTION AT TltffLJt. _ 


CIH2 J INTERSTITIAL -CONCENTRATION OF H2 
CI02 J INTERSTITIAL CONCENTRATION OF 02 

CIH20 : INTERSTITIAL CONCENTRATION OF H20 

T 5 INTERSTITIAL TEMPERATURE 

P : INTERSTITIAL PRESSURE 

002 : KNUDSEN 

DH2 ! DIFFUSION 
DH20 : COEFFICIENTS 
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RT02 (.N ) J ARRAYS CONTAINING. 

RTH2(N) : CHARACTERISTIC ROOTS OF 

RTH20 ( N ) 5 THE E Q UATION X*COT (X) + (PHl*A-l) a _0_ 

RTt(N) : HEAt TRANSFER (RTT) » AND MASS TRANSFER ( RT02 » RTH2 # R TH20 ) 

LAMDA(I) J NORMALIZED RADIAL POSITIONS 
LAMDAOU) WITHIN THE CATALYST P ARTI CLE 

GRT : GREEN’S F UNCTIONS USED IN PARTICLE 

GRTI : TEMPERATURE DISTRIBUTION CALCULATIONS 

GR02 J GREEN’S FUNCTIONS USED IN PARTICLE CALCULATION 

GR02I : OF OXYGEN DISTRIBUTION 

GRH2 : GREEN’S FUNCTIONS USED IN PARTICLE CALCULATION 
6RH2I : OF HYDROGEN DISTRIBUTION 

GRH20 J GREEN’S FUNCTIONS USEO IN PARTICLE CALCULATION, _ ___ 

GRH2QI ! OF "H20 DISTRIBUTION 

RADP ! PORE RADIUS OF CATALYST PELLET 

TOL : TOLERANCE used FOR CONVERGENCE CHECK ON 02 PROFILE 

MAX I MAXIMUM NO. OF ITERATIONS ALLOWED 

NPP : NO. OF. POINTS/PROFILE WITHIN CATALYST PELLET 

NRTT : MAXIMUM- NUMBER OF TERMS USED IN SERIES EXPANSION FOR 
NRT02 J OBTAINING GREEN’S FUNCTIONS ... 

NRTH2 : NRTT TEMPERATURE# NRT02 OXYGEN# 

NRTH20! NRTH2 HYDROGEN# NRTH20 WATER 

DJ\F(RADP#.T#WT) = 237,*RADP*SQRT(T/WT) 

DF(PC»PCH2»WT #WTH2#T »P»TC#TCH2) = ,442E-6*CBRT ( PC*PCH2) +SQRT ( 1 ,/W' 

X, +1./WTH2)*T**1,823 / (P*(TC*TCH2)**.495) 

R = A 

PT = DELTT 
CAPPA = KAPPA 
NPl = NPP“1 
ITER = 0 
KGUNT = 0 
NCON = 0 
IRATE = 0 
NEWOT = 0 

IF (PASSNO. EG, 1. AND. IT. EQ.l) GO TO 21 
IF ( I TSTEP * EQ , 1 . AND . PASSNO • EG . 1 ) GO TO 20 
IF { I TSTEP , EQ , 1 , AND . JDP . EQ . 0 ) GO TO 25 
IF (JDP, EG. 0) GO TO 28 
TPSX = DPTPS+50. 

IF (PT(IZ) .GT.TPSX) GO TO 23 

CHECK TIME INCREMENT WITH PREVIOUS DELTT MUST CALCULATE NEW GREEN’S . 

FUNCTIONS IF DIFFERENT 
28 IF (DELTP-DELTT) 21 #20 #21 
25 IF (DELT (PASSNO) -DELT (PA5SN0-1 ) ) 2l#20#21 
NEWDT IS INDICATOR FOR NEW TIME INCREMENT ... 

NEWDT = 0 INCREMENT UNCHANGED 

NEWDT = 1 — NEW INCREMENT IN USE 
23 OPTPS = PT(IZ) 

21 NEWDT = 1 

20 IF (PUNCH. EG. 2. UR. PUNCH, EQ. 3) GO TO 93 
IF (PASSNO. EG. 1. AND. IT, EG, 1) GO TO 43 

c set particle concentration, temp# rhet profiles from surface values 
C CALCULATED AT. PREVIOUS TIME INCREMENT 
93 DO 42 I=1»NPP 

CP02(I#1) = C02 ( I # I Z) 

CPH2 ( I # 1) = CH2 ( I # IZ ) 

CPH20 ( I # 1 ) = CH20(I#IZ) 
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TP(I,1> s TPP(I,IZ) 

RHET (I»l) = ALPHA.K*CP02ll. e J,l ** . 8*£XP (-flAMMA/JE LU1U. . - 

42 RHET { I ? 2) = RHETdd) 

. _ GOTO 44 

43 DO 1 IrlfNPP 

. CPQ2U>1) = 0S02 

CPH2(I»1) = CSH2 

CPH2Q(I,1). .= CSH20. . . . 

RHET(X,1) = ALPHAK*CP02(I»l)**.8*EXP(-GAMMA/TP(Iri) ) 

AU£UXjl2Jl =..BHEX(.I.*Ju). 

1 TP(I,1> s TPS 

44 NRTT 5 A*SQRT (23*03/ (KAPPA^DELTT) J/PI+.5 
IP (NRTT.GT.NRTS) NRTTrNRTS 

PHIT s HC/KP 

SOLVE FOR ROOTS OF EQUATION X*COT(X)+C =0. ... C=HC/KP*A-1. 

. X S PHIT*Arl« . . _ . - 

C HEAT TRANSFER 

RLOWER = PI/2. 

RUPPER = PI 
DO 2 N=1,NRTT 

CALL ROOT (RLOWER»RUPP£R,C f RTT(N> »N] 

RLOWER = RLQWER+Pl 

2 RUPPER = RTT (N)+Pl 

DK02 = DKF(RADP,T(IZ)»M02) 

DKH2 = DKF(RADP,T(1Z) »MH2) 

IF CP(IZ) .LT..0G1) 00 TO 13 
IF (PUNCH. EQ. 2. OR. PUNCH. EG. 3) GO To 71 
IF (PASSNO.EQ.i. aND.IT .EQ.1} GO TO 26 
71 TEMP = PT(IZ) 

60 TO 24 

26 TEMP = T(IZ) 

24 U02 r DF(PC02fPCH2»i'i02fHH2»TEMP,P(IZ) fTC02rTCH2) 

UH2 = DF(PCH2#PCH2*MH2»MH2»TLMPrPUZ)rTCH2»TCH2) 

D02 = D02*U.-E,XP(-DK02/D02)) 

UH2 = DH2* ( 1 « “EXP ( -UKH2/DH2 ) ) 

GO TO 50 
16 Do 2 = DKQ2 
DH2 s DKH2 

50 PHI02 = KC02/DQ2 
PHIH2 = KCH2/DH2 

C SOLVE FOR ROOTS OF THE EQUATION X*COT(X)+C = 0. C=KCH2/0H2*A-1 

C C=KC02/D02*A-l.t 

C C=KCH20/DH20*A-1. 

6 NRT02 = ASSORT (23.03/ (D02*D£LTT ) ) /P 1+ ,5 

NRTH2 = A*SQRT(23.03/(DH2*DELTT) )/Hl+,5 
IF (NRT02.GT.NRTS) HRT02=NRT$ 

IF (MRTH2.6T.NRTS) NRTH2=NRT5 
KLOWER = PI/2, 

RUPPER = PI 

C MASS transfer 

CA = PHI02+A-1. 

DO 3 N=1,NRTQ2 

CALL ROOT (RL0tfERrRUPPER»CA»RT02(N) »N) 

RLOrtEK = RLOWER+Pl 

3 RUPPER = RT02 (N ) +PI 
CB = PHIH2*A-1. 

RLOWER = PI/2. 

RUPPER = PI 

DO 4 Nzi r NRTH2 

CALL ROOT ( RLOWER » RUPPER »CB»RTH2(N) #N) 

RLOWER = RLOWER+Pl 

4 RUPPER = RTH2(N)+PI 
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c SET UP RADIAL STATIONS WITHIN PARTICLE ON FIRST ENTRY TO ROUTIN! 

DEL a 1. /FLOAT (NPP) 

LAMDAll) a 0. 

LAMDAO (1} _a 0, 

DO a f=2»NPl 

LAMDA(I) = LAMDA (I-D+DEL 

8 LAMDAO (I) a LAMDA (I) 

LAMDA(NPP) = 1,0 
LAMDAO (NPP) a 1,0 

CALCULATE INITIAL GUESS FOR P ART ICLE CONCENTRATION UISTRIBU » ium ui- u* 

19 NT = NRT02 

DO 150 L=1,NPP 
DO IOC I=2»NPP 

IF ( N.EWDT ..Q . QR..U P. . EQ ,Xt.OR . ITER ..GTU.Q1. SO ...TO 9 . 

IF (I.LT.L) GO TO 7 

CALL GREEN (LAMDA(L) ,LAMDA0(I) ,UT02,D02,PHI02,GR02CL»I) » 

X GK02ICUI) ,GR02TI(L,I) »DT,R,NT,L) 

GO TO 9 

C GREEKS FUNCTIONS SYMMETRIC W.R.T. LAMDA, LAMDAO G (L,L0 ) =G (LO ,L) 

7 GRQ2 (L» X ) a GR02(.IiL) 

GR02I (L, I ) a GR02I ( I »L ) 

GR02TI (L,I) a GR02TI(I,L) 

9 IF (L.EQ.l) GO TO 30 

F (I ) = 2.*LAMDAQ < I ) /LAMDA (L) * (-A**4/C02**2* (RHET ( 1 ,2) -RHET ( I » 1 ) )/ 

X DELTT*GR02TI (L, I ) -A**2/D02*RHET ( I , 1 ) *GR02I (L, X ) + 

X (CP02(I,1)-C102(IZ))*GR02(L,I)) 

60 TO 100 

30 F ( I ) a 2,*{-A#*VD02**2*(RHET(I,2)-RHET{I#l))/DELTT*GR02TI(LfI)-A* 

X *2/D02*RHET (1,1) *GRQ2I (L,I)+ (CP02 { 1 , 1 ) -C 102 < IZ ) ) *GR02 C L » I ) > 

100 CONTINUE 
RSUM - 0, 

DO 133 1=1 f.NPl 

133 RSUM = RSUM+.5*(F(I )+F(I + l) )*(LaMDAU+1)-LAMDA(I) ) 

CP02 (L ,2) = CI02{IZ)+RSUM 
IF (CP02(L,2).LT,0.) CP02(L,2) a 0. 

C02 (L , IZ) = • CP02 (L , 2 ) 

150 CONTINUE 

C CALCULATE PARTICLE TEMPERATURE DISTRIBUTION 
C 

15 NT a NRTT 

UO 650 L=i»NPP 
DO 600 1=2, NPP 

IF INEWDT.EQ.O. OR. ITP.EO.IT. OR, ITER. GT.O) GO TO 10 
IF (I.LT.L) GO TO 11 

CALu GREEN (LAMUA(L) .LAMDAO ( I ) »RTT »CAPPA»PHlT ,GRT (L, I) , 

X GKTKL*I)»GRTTI<L,I) ,DT»R,NT»L) 

GO TO 10 

C GREEN’S FUNCTIONS SYMMETRIC W.R.T. LAMDA AND LAMDAO G(L,L0)=G(L0,L) 

11 GRT(L.I) = GKTU.L) 

GRTKL.I) = GRTKI.L) 

GKTTKL.I) = GRTTKI.L) 

C CALCULATE INTEGRAND FUNCTION 
10 IF (L.EO.l) GO TO 31 

F ( I ) = 2 . +LAMDA0 ( I ) /LAMDA (L) * { ( -MR*A**4/ (KAPPA*KP ) ) * C RHET ( 1 , 2 > “ 

X RHET (1,1)) /DELTT*GRTT I ( L , I > -HR*A**8/KP*RHET (1,1) *GRT I ( L , I ) + 

X ( TP ( 1 ,1 ) -T ( I / ) ) *GRT ( L , I ) ) 

GO TO 600 

51 F ( I ) = 2,*((-HK*A**4/CKAPPA*KP) ) ♦ (RHET ( I , 2) -RHET ( I , 1) J/DELTT* 

X GRTTI (L , I ) -HR*A**2/KP*RHET (1,1) *GRTI (L , I ) + (TP ( I , 1) -T ( IZ)) * 
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X GRT (L»I) ) 

6QQ CONTINUE 

c calculate integral expression for given lamda, ALL LAMDAO'S. 

RSUM S Q, 

DO 55 1=1 » NP1 

55 RSUM s RSUM+.5*tF{I}+FU+l))*tLAMOA-im)-LAMDA(I)) 

TP (L ,2) = TCIZJ+RSUM 

TPP (LilZj „ = TP.CL,2) . - 

RATE(L) = ALPHAK*CP02(L»2)**,8«EXP(-GAMMA/TP(L»2) ) 

65CL ..CONTINUE - — - 

CHECK FOR CONVERGENCE 
DO 49 I=1»NPP 

IF <ABS<RHETU»2)-RATE(I> )-.02*RATEC I ) ) 49*49,32 
32 IF {RATEtI)-.02*RATE{NPP)) 49,49,27 . 

49 CONTINUE 

CONVERGENCE — CALCULATE CONCENTRATIONS FOR H2t H20 
CALCULATE PARTICLE CONCENTRATION DISTRIBUTION FOR H2 

NT = NRTH2 
DO 550 L=i»NPP 
DO 500 1=2, NPP 

IF (NEWOT.EQ.O.OR.ITP.EO.IT) GO TO 13 
IF (I.LT.L) GO TO 12 

CALL GREEN (LAMbA (L) f LAMDAO { I ) »RTH2,0H2,PHIH2,GRH2 (L» I ) » 

X GRH2I (L# I) »GKH2TIIL» I ) »DT»fi»NT»L) ■ • 

GO TO 13 

L GREEN’S FUNCTIONS SYMMETRIC W.H.T. LAMDA AND LAMDAO - — G (L*L0 ) =G (LO «L) 

12 GKH2 ( L , I ) = 6RH<d ( I , L ) 

GRH2I (L» I ) = GKH2I ( X ,L) 

GKHLTI(L,I) = GRH2T I ( I ,L ) 

13 IF (L.EQ.l) GO TO 33 

F ( I ) = 2.*LAMDAO(I)/LAMDA(L)*(-A**4/DH2**2*2.*MH2/M02* 

X . (RH£T(I,2)-RHET{I,1) >/DELTT*GRH2TI (L»I )-A**2/DH2*2.*MH2/ 

X M02*RHET(I»1)*GRH2I(L»I)+(CPM2(I,1)-CIH2(I2> )*GRH2(L»I)) 

GO TO 500 

55 FID = 2»*{-A**4/DH£**2*2«*M!i2/M02*(RHET ( I »2)-RHET (1,1) ) /DELTT 
X +GRH2TI ( L , I) -A**2/DH2*2 . *MH2/M02*RH£T [ 1 , 1 ) *GRH2 1 C L » I ) 

X + (CPH2 { I » 1 ) -CIH2 ( 12) ) +GRH2 (L» I ) ) 

500 CONTINUE 
RSUH = 0, 

DO 75 1=1 ,NP1 

75 RSUM = RSUM+ ,5* (F ( I ) +F (1+1 ) ) + (LaMDA ( 1+1 ) —LAMDA ( I ) ) 

CPH2(L,2) = CIH2CUJ+RSUM 

CH2(L»I2) = CPH2 (L, 2) 

350 CONTINUE 

CPH20(NPP,2) = 0. 

C ■ SAVE SURFACE VALUES OF TEMP, CONCENTRATION 
CS02 = CP02 ( NPP , 2 ) 

CSH2 = CPH2 (NPP,2) 

CSH20 = CPH20(NPP»2> 

TPS r TP(NPP,2) 

ITER = ITER+1 
IF (TPS-UCTEMP) 77 , 77 , 78 
78 NCOiM = 1 
GO TO 99 

77 IF (CS02) 880,76,76 

76 IF (CSH2) 881 r 99 ,99 
99 RETURN 

C 
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C NO CONVERGENCE ..... SAVE CALC. RHET PROFILE AND REPEAT TEMP AND CP02 CALCS 
C 

27 IF (KOUNT.EQ.l) GO TO 38 

IF ( JTER.GT.O) GO TO 35 

IF ( RATE ( NPP ) . GT . RHEf (NPP , 2 } > I RATE=1 
GO TO 36 

35 IF ( (RATE (NPP) • GT • RHET ( NPP , 2 ) . AND. IRATE. EG. 0) .OR. (RATE(NPP) ,LT . 

X RHET (NPP .2) .AND, IRATE ,EQ.i ) ) GO TO 38 

36 DO 37 1=1 » NPP 

37 BHEJJA,2J s RATEJI ) 

GO TO 41 

38 00 39 1=1, NPP 

39 RHET (1,2) = , 5*(RATE(I)+RHET (1.2) ) 

kount = 1 

41 ITER = ITER+1 

IF ( I TER , GE . MAX ) GO TO 86 
GO TO 19 

86 CS02 = CP02 (NPP, 2) 

CSH2 = CPH2(NPP,2) 

C5H20 = CPH20(NPP,2) 

TPS = TP (NPP* 2) 

IF (T.PS-DCIEMP) 161*161,162 

162 NCON = 1 
RETURN 

161 IF (CS02) 880*90*90 

90 IF (CSH2) 881,92,92 

92 RETURN 

ERROR STATEMENTS FOLLOW ... 

880 WRITE (6,890) CS02,CSH2,CSH20,TpS, IT,PASSNO 

890 FORMAT CO* 5X, 120{'*') // 58X * ERROR ERROR* / 38X 'PARTICLE 

XSURFACE CONCENTRATION OF 02 HAS GONE NEGATIVE' / 24X *CS02 =*, 

X Ell. 5, 5X,_ *CSH2 =* ,E11 .5, 5X, *CSH20 =*,£11.5, 5X, *TPS = »., 

X Ell. 5 / S2X 'TIME INCREMENT ', 13, * OF PASS', 13 // 120C*') ) 

CALL EXIT 

8dl WKl IE (6,891) CbH2, CS02,CSH20, TPS, IT rPASSNO 

891 FORMAT ('O' 5X, 120('*') // 58X 'ERROR ERROR* / 38X 'PARTICLE 

XSURFACE CONCENTRATION OF H2 HAS GONE NEGATIVE' / 24X 'CSH2 = ', 

X Ell ,5* 5X» ' C502 ='»E11.5» 5X, 'CSH20 =',EU.b, SX, 'TPS =*, 

X Ell . 5 / 52X 'TIME INCREMENT *, 13 , ' OF PASS', 13 // 120C'**) ) 

CALL EXIT 
END 
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XXX XXXXXX X XXX 


SUBROUTINE PCALCS {TPS»CS02,CSH2,CSH20,CI02PR# IZ»MINT> 

ROUTINE TO DO PARTICLE SURFACE TEMPERATURE CALCULATIONS AT STEADY STATE 
INTEGER PUNCH 

. REAL _KCQ.2L/J&JH&rKCH20.f MQ2 f.MH2i.MH20._„. 

PARAMETER NAX=70, NTI=40 
PARAMETER NPP=3Q, NRTS=2QQ 

COMMON /CB1/ TPP(NPPfNAX) >C02(NPP»NAX) »CH2<NPP»NAX) ,CH20 (NPP#NAX) # 
T<NAX) ,P(NAX) »PT (NAX) ,H(NAX) ,CIH2(NAX) ,CI02(NAX) f CIH20(NAX)» 
CIHHOS(NAX) »CIH20L (NAX) rCIH20V(NAX) rCIHE(NAX) »DELT (NTI ) t 
PPT.(NAX) 

COMMON /CB3/ ALPHA » GAMMA , K ARPA , HC » KP t A r DELTT t KCH2 , KC02 r KCH20 » TCH2 1 
TC02 1 TCH20 1 PCH2 1 PC02 f PCH20 » MH2 1 M02 » MH20 1 MHE t HR » HA » HA1 » H A2 , 

A W t C W , MW t DC t T A ♦ D0H2 , DQ02 , 00H20 » ALPHAK > AP » DELTA t G t MU » HS r HL , 

HV i TW # RHOM » WH2 , W02 , WH20S f WH20L i WH20V , WHE t CFBAR « MB AR • C I CPI , 
ClCP2»ClCP3»RHO»RHOV»DELHC»TlTPSfDTDZ,DHDZ#TSfTL»DPDZ»DClOZ» 
DC2DZrPC3SDZ»0C3LQZ»DC3VD2»DC.4DZ».RHQPrCPiDELTP»DPTPSiSDTlME, 

PF r TSS # AC r VC » PREST t GSS , MB ARSS 

COMMON /CB4/ NOF Z f NP ASS i NAX I AL » I NPERP » MODPC t MODNO , PUNCH » I TSTEP * 

UDP» ISHUT i IHQT 

COMMON /TABLES/ AVSZ<40) ,APVSZ(40) »DELVSZ(40) #HRVSTC36) ,TCi(22> , 
TC2(20) »TC3V(20) ,TC4US) ,MU1VST{34) »MU2VST(34) »MU3VST<32) * 
MU4VST (34) iCFIVST (34) /CF2VST (34) »CF3VST (34) »CF3SVT (203 * 

CF4VST (20 ) ? VPVST ( bO ) >TVSVP<48) rDHCVST<24) 

CALL UNBAR (HRVSTi l t PT ( IZ) » 0. »HR»KK) 

GO TO (l»l'lr2)' PUNCH 
IF UZ.NE.2) GO TO 3 
C102CIZ) = CI02(IZ-i) 

go ro l 

CI02(IZ) = CI02PR 
IF (CI02(XZ).LT.O«) CI02(IZ)=0. 

DTPSDT s 3./(A*RH0P*CP)*(HC*<T(IZ)-PT(lZ) ) -HR*KC02*CI02 ( 12-1 > ) 

L'XTRM = U.-EXP(-3.+HC*DELTT/(A*RH0H*CP)))/(3.*HC/{A*RH0P*CP) ) 

TPS = PTtI2)+DTPSDT*(l.-EXP(-3.+HC*DtLTT/CA*RH0P*CP) ) )/ 

X (3,*HC/(A*RH0P*CP) ) 

CS02 = 0. 

CSH2 = CIH2(IZ-1)-(KC02/KCH2*Mh2/(M02*,5)*(CI02(IZ~1)-CS02) ) 

CSH20 = ClH20VUZ-i)-(KC02/KCH2o*MH20/(M02*,B)*(CS02-CI02(I2-l) ) ) 

RETURN 

END 



=o.cr> 


0W3LD2 = DW3DZ*( (H(IZr , l)-HS)/(HSL'~HS) J+DHDZ/ (HSL-HS) 

J1W3.VDZ = 0.. 

W3 = 0, 

’ RHOi sVh0/(1.-WH20S-WH20L) 

DWJDZ = DW3SDZ+DW3LDZ 
RHOa = RHOV/RHO 
60 TO 6 

C LIQUID REGION 
.3 HRD = ,Q. 

DELTAH = DELHC 
DW3DZ = DW3DZ/.FACT3 
DWDZ3 = 0, 

DW3SDZ = 0, 

DW3LD2 = DW3DZ 
DW3VDZ = 0.. 

W3 = 0. 

HFACT = 0. 

RHOI = RHO/ ( 1 « -WH20L ) 

DWJDZ = DW3DZ 
RHOi = RHOV/RHO 
GO TO 6 

C LIGUIU-VAPOR REGION 

4 HRD = 0, 

HLV = H(IZ-l)+DtLHC*WH20L 

DELTAH = DELHC*(HLV-H(IZ-1) )/(HLV-HL) 

DW3DZ = DW30Z/FaCT3 
DW3SDZ = Q. 

DW3LDZ = DW30Z*({HLV-H(IZ-1))/(HLV-HL))-DHDZ/(HLV-HL) 

DW3VDZ = DW3DZ*((H(IZ-1)-HL) /{HLV-HLJ ) +DhDZ/ (HLV-HL) 

DWDZ3 x Ofc3VQZ 
Vi 3 - WH20V 

HFACT s -KC02*AP*CICP3*DELHC/{e*(HLV-H(IZ-l) ) ) 

RHOI = RH0/U.-WH20L) 

DWJDZ = UW3LDZ 
RhOZ x RHOV/RHO 
GO TO 6 

C VAPOR REGION 

5 HRD = HK*RH0M*UELTA/G 

DELTAH = 0, 

0WDZ3 = DW3DZ 
DW3GDZ = 0, 

DW3LDZ = 0. 

DW3V0Z = DW3DZ 
W3 = WH20V 
HR ACT = Q. 

RHOI = RHO 
DWJUZ = 0. 

RHOa = 1, 

6 FACTOR = a.-EXP(-(hC*Ap/{G*CF8AR)+HFACT)*DELTZ))/((HC*AP/(G* 

X CFBAR)+HFACT)+DELTZ) 

DhDZ .= {~HRD-HC*AP*TITPS/G-KCH20*AP*CICP3*DELTAM/G-4,*HC 
X *{T(IZ-1)-Tfl)/(G*DC) )*F ACTOR 

DTDZ = DHDZ/CFBAR 
GO TO (7»8»7 » 8r7 ) t REGION 

7 TDZ = DTDZ/T ( IZ-1 ) 

GO TO 9 
TDZ = 0. 

DM8KDZ = “MBAR*(DW1DZ/MH2+DW2DZ/M02+DWDZ3/MH20)/(WH2/MH2+W02/M02+ 
X W3/MH20+WHE/MHE) 

URHODZ = RH01*{DWJDZ+RH02* IDPDZ/P ( IZ~i ) +DMBRDZ/MBAR -TDZ)) 

RCOR s -RH01 /RHo*(DwJDZ+RH02*(0PDZ/P(IZ-1)+DNBRDZ/MBAR-TDZ) ) 


77 



DRODZP = DRHODZ*(1.-EXP(-RCOR*DZPR))/{RCOR*DZPR> 

DC2DZP = RHO*DW2DZP+W02*DRODZP 

CI02PR = CI02<IZ-1)+OC2DZP*DZPR 

DRHODZ_ s DRH0DZ*{1.-EXP(-RC0R*DELTZ1)/.(R£0R*DELJZ) 

DClDZ = R'H0 *DWIdZ+WH2*DRH0DZ 

DC2DZ = RH0*DW2DZ+WQ2*DRH0DZ 

DC3SDZ = RHO*DW3SDZ+WH20S*DRHOOZ 

DC3LDZ r RH0*DW3LDZ+WH20L*DRH0DZ 

DC3VDZ = RHO*0W3VbZ+WH2OV*DRHODZ 

PQ40Z = WHE*DRHODZ 

RETURN 

end 
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SUBROUTINE GREEN ( LAMDA , LAMDA0»RT» D # PHI , GR » GR I » GRTI , OT » A r NR , LL ) 

ROUTINE TO CALCULATE GREEN’S FUNCTIONS NEEDED FOR TRANSIENT 
PARTICLE ANALYSIS. 

. .REAL LAMDAjlLAMDAG 
DIMENSION RT(1 J 
SUMi a 0. 

SUM2 = 0, 

SUM3 = 0, 

SUMLl a 0. 

S.U.M.L2 = 0, 

SUML3 = 0. 

DO 2 Ial#NR 

R.TERM = (RTU)**2+(PHI*A-1.)**2)/(RT(I)**2+PHI*A*{PHI*A-1.) > 
rad = RT<I)*LAMQA 
RADO = RT(I)*LAMDAQ 
S1NL = SIN (.RAD) 

SINLO a SIN(RADO) 

EXTRM1 a EXP ( - ( RT ( I ) /A ) **2*D*DT ) 

EXTRM2 a EXTRM1/RT { I ) **2 

EXTRM3 = U.-EXTRM1)/RT(I)**4 

IF (LL.EQ.l) GO TO 1 

SUMI = SUMURTERM*SIHL*SINLQ*EXTRM1 

SUM2 = SUM2+RTERM*SINL*SINLG*EXTRM2 

SUM3 = SUM3+RTERM*SINL*SINL0*EXTRM3 

GO TO 2 

1 SUMLl = SUML1+LAMDA0*RTERM*RT(I)*SINL0*EXTRM1 
SUML2 a SUML2+LAMDAO*RTERM*1./RT(I)*SINLO*EXTRM1 

SUML3 = SUML3+LAMDA0 *RTERM*1 «/RT ( I ) **3*SINLQ* ( 1 .-EXTRM1 ) 

2 CONTINUE 

IF CLL.EQ.l) GO TO 4 
GR = SUMI 

IF (LAMDAO.GT. LAMDA) GO TO 3 

S = 0.5*LAMDA*LAMDA0*<1,/LAMDA-(PHI*A-1.)/<PHI*A) ) 

GRI a S - SUM2 

GKTi = D*DT/A**2*S - SUM3 

RETURN 

3 S = 0.5*LAMOA*LAMDAO*(1./LAMDAO-(PHI*A-1.)/(PHI*A)) 

GRI a S - SUM2 

GKTI = D*DT/A**2*S - SUM 3 
RETURN 

4 GR = SUMLl 

S = LAMDA0/2.-LAMDA0**2/2.*( (PHI+A-i . ) / (PHI*A) ) 

GRI = S - SUML2 

GkTI = D*DT/A**2*S - SUML3 

RETURN 

END 
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SUBROUTINE ROOT (RLOWER,RUPPER,C,RTO,N) 

ROUTINE TO FIND ROOTS OF THE CHARACTERISTIC EQUATION X*COTCX)+C = 0, 
modified method of SUCCESSIVE SUBSTITUTIONS 

1=0 

RTLO = RLOWER 
RTH1 = RUPPER 

1 IF CABS (RTHI-RTLO) •LT«2»E-5) GO TO 9 
I = 1+1 

IF { I *GT 0 25) GO TO 7 
RTO s IRTL0+RTHU/2. 

RT = -C+TANCRTO) 

IF CRT-RUPPER) 2,2r4 

2 IF (RT-RLOWER) 5,3 ,3 

3 IF (AB5((RT-RTQ)/RT)-,G0D05) 9,9,6 

«f RTLO = RTO 

GO .TO 1 

5 RTH1 s RTO 

GO TO 1 

b IF CRT.LT.RTU) RTHI=RT0 

IF CRT. GT. RTO) RTLO=RTO 

go ro l 

7 WRITE (6,8) N 

a format <*o» // 2ix, *no convergence on root number*, 13,* after 25 

XITERaTIONS USE VALUE CALCULATED ON 25TH ITERATION') 

9 RETURN 

END 
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CIQ2(K) s CCI02( I) 
CIH20.CKJ = CCH20(I) 
K s k-1 

-5fl ,£M.UMJ£ 

RETURN 

END 



SUBROUTINE INSERT (ZIF,I2) 


C 

C 

X. 


50 


51 


ROUTINE TO INSERT INTERFACE AXIAL POSITIONS INTO ORIGINAL Z ARRAY 
PARAMETER NAX=70 r NT I =40 

~ COMMON /CB1/ TPP ( NPPTNAX )’ rc 02 { NpP » NAX > »CH2 (NPP»NAX) »CH20(NPP»NAX) t 
X , T(NAX) »P(NAX) »PT (NAX) rH(NAX) rCIH£{NAX) »CI02(NAX) ,CIH2Q(NAX) » 

X. CIH20S{NAX) »CIH20L(NAX) fClH20V(NAX) »CIHE(NAX) »DELT (NTI > » 

X . PP.T(NAX) 

COMMON /CB2/ 2(NAX)rZSSL(NTl)#ZSLL(NTI) f ZLLV(NTl),ZLW(NTI} 

XOMMQN /Cfi3/. ALPHA * GAMMA » KAPPAiHC t K.P * AtOELT T-t KCH2 #KCD2.*iCCH20 « T CH2* 
X TC 02 1 TCH20 , PCH2 , PC02 1 PCH20 1 MH2 1 M02 1 MH20 » MHE p HR > H A r HA1 r H A2 , 

X . AWjCy»MW#DC»TA»DOH2,D002,uOH20 ? ALPHAK,AP,DELTA»G»MU#HS»HL > 

X ■ HV f TWrRHOM»WH2»W02rWH20S»WH20L»WH20V»WHErCFBARfMBAR»ClCPl, 

X . C ICP2 # C ICP3 f RHO » RHOV # DELHC t T ITPS # DTD2 # DHDZ » TS » TL # DPOZ , DC1DZ p 

X . DC2DZ > DC3SDZ»DC3LDZ,DC3VDZfDC4DZ»RH0P,CP»DELTP#DPTPS»SDTIME» 

. X PF i.TSS » AC * V£.» PREST » GSS * MBARSS .... 

COMMON /C84/ NOFZ pNP ASS » NAXIAL , INPERP p MODPC pMODNO , PUNCH p ITSTEP , 

X JDPp ISHUT p IHOT 

DIMENSION ZZ(NAX) pPTZ(NAX) pTIZ(NAX) ,CI02Z(NAX) pCIH2Z<NAX> 

DO 50 1=1 p NOFZ 
PT2( I) = PT( I ) 

. . UZ1I} .K XU) 

CI02Z(I) = CI021I) 

CIH2Z(I) = CIH2(I) 

ZZil) = Z(I) 

DO 51 I=IZpN0FZ 
PT ( 1+1 ) = PTZ ( i ) 

m+u = t.iz{.u 

CI02 ( 1 + 1 ) = CIQ2ZU) 

CIH2(I+I) = CIH2Z(I) 

Z(I+1) = ZZtl) 

ZUZ) = ZIF 
NOFZ = NQFZ+1 


52 RETURN 
END 
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SUBROUTINE DELETE (PASSNO) 

C 

C ROUTINE TO QELETE INTERFACE AXIAL STATIONS FROM INITIAL INPUT 2 ARRAY 
C 

INTEGER PASSNO . 

EARJW&IERJIMSiat Jill 54a*. NPPs3n 

COMMON /C81/ TPP(NPP,NAX) * C02 (NPP»NAX) ,CH2(NPP,NAX) ,CH20(NPP»NAX) » 

X T(NAX) fPlNAX) ,PT(NAX) ,H(N aX) ,CIH2(NAX) ,CI02 (NAX) ,CIH20(NAX) » 

X CIH20S (NAX) »CIH20L(NAX) »ClH20V(NAX) »CIHE(NAX) »DELT(NTI) * 

X . PPT(NAX) 

COMMON /CB2/ Z (NAX) * ZSSL (NTl ) »ZSLL (NTI ) tZLLV (NTI ) »ZLVV (NTI ) 

COMMON /CB4/ NOFZ t NP ASS , NAXI AL , INPERPjMODPC * MODNO , PUNCH tIT£TEP 

X JDP» I SHUT 1 1 HOT ' • 

O.IMENS ION 22 (NAX) »TT ( NAX ) f CC I H2 ( NAX ) , CC 102 ( NAX > , CCH20 ( NAX ) 

C SAVE CURRENT 2 ARRAY AND CORRESPONDING TEMPS* CONCENTRATIONS 
00 5 1=1 t NOFZ 
22(1) = Z(I) 

IT. LI) .= Til) 

PPT(I) s PT{ I ) - 
CCIH2(I) = CIH2(I) 

CCI02(I) = CI02 { I ) 

5 CCH20(I) = CIH?.0(I) 

C LOCATE INTERFACE POSITIONS 
ROF.Zl = NOfZ-1 
00 10 I=l»NOFZi 
IF (Z(I)-2SSL(PASSN0) ) 10*1*10 
10 CONTINUE 

RETURN 

1 IF1 = I 

DO 20 I=IF1»N0FZ1 
IF (Z (I )-ZSLL (PASSNO) ) 20*2*20 
20 CONTINUE 

IF2 = NAXIAL+1 
IF3 = NAXIAL+1 
IF4 r NAXIAL+1 
GO TO 6 

2 IF2 = I 

DO 30 I=IF2*N0F21 
IF (2(1) -ZLLV ( PASSNO ) ) 30*3,30 
30 CONTINUE 

IF3 = NAXIAL+1 
IF4 = NAXIAL+1 
GO TO 6 

3 IF3 = I 

DO 40 I=IF3,N0F4l 
IF (Z (I )-ZLVV (PASSNO) ) 40,4,40 
40 CONTINUE 

IF4 = NAXIAL+1 
GO TO 6 

4 IF 4 = I 

C ADJUST 2 ARRAY TO INITIAL VALUES 

6 K = NAXIAL-4 

DO 50 I=N0FZ,l*-i 

IF (I.EQ,IF4,0R.l,Efc,IF3.0R.I.Eu.IF2.0R.I.EU.XFl) GO To 50 
Z(K) = ZZ(I) 

T (K) = TT ( I ) 

PT(K) = PPT(I) 

CIH2(K) = CCIH21I) 
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SUBROUTINE OVSTEP l ZEND *IZ*PASSNO* REGION) 

C .... 

C THIS ROUTINE MAKES APPROPRIATE ADJUSTMENTS TO TEMPERATURE* PRESSURE* 
C. CONCENTRATION. EICm WHEN AXIAL. POSITIONING. HAS.- QVERSTEPPED. INTERFACE 
C BOUNDARIES FROM ONE REGION TO ANOTHER 

integer 1*Xs1n67reg ion 

PARAMETER NAX=7Q* NTI=4Q 
PARAMETER NPP=30 

COMMON /C(31/ TPP ( NpP t N AX > * CO?, (NPP * NAX ) i.CHS t NPP * NAX ) * CH20 INPP iNAX } * 
X T (NAX) *P(NAX) *PT(NAX) *H(NAX) *CIH2(NAX) *CI02(NAX> *CIH20(NAX> * 

.X...- . .CIH2QSlNAXI».CIH20LtNAX}.»OlW2QVINAXl.*.CIHECNAXl*DELTtNTll^- 
X PPT(NAX) 

COMMON /CB2/ Z(NAX) *ZSSL(NTl) »ZSLL(NTI> »ZLLV(NTI) *ZLVV(NTI} 

COMMON /CB3/ ALPHA , GAMMA * KAPPA * HC * KP * A * DELTT * KCH2 * KC02 * KCH20 * TCH2 * 
X . TC02 » T CH20 » PCH2 * PC02 * PCH20 * MH2* M02 * MH20 » MHE* HR * HA *H Al * HA2 * 

X A W » CW * MW * DC * T A » D0H2 * D0O2 » UOH20 * ALPHAK * AP * DELTA , G * MU * HS • HL * 

X . .. HV* T W i RHOM* WH2 * WQ2 * WH20S *.WH201_» WH2Q V * WHE * CEBAR* MBAR«.CICPI*.. 

X CICP2*CICP3* RHO * RHO V * DELHC * T I TPS » DTDZ * DHDZ , TS » TL , DPDZ * DC lOZ » 

X DC2DZ*DC3SDZ»DC3LDZ»DC3VDZ#DC4DZ,RH0P*CP,DELTP*DPTPS»SDTIME* 

X PF»TSS*AC,VC*PREST*GSS*MbrtRSS 

COMMON /CB4/ NOFZ * NP ASS * NAXI AL » INPERP * MODPC * MODNO * PUNCH * ITSTEP » 

X JDP* I SHUT* I HOT 

GO TO Ui2i 3*4»5)» REGION 
c solid/solid-liquiu interface 
1 DELTZ = ITS-TUZ) )/DTDZ 
TUZ) = TS 

H ( 12.) = H(IZ)+DHDZ*UELTZ 
HS = H { IZ) 

CIH20SUZ) = CIH20S(IZ)+DC3SCZ*DELTZ . 

ZSSL(PASSNO) = Z { IZ) +DELTZ 
ZIF = ZSSL(PASSNO) 

PT(IZ) s (ZIP-Z(IZ-I) )*IPT(1Z)-PT(IZ-1) )/(Z(IZ)-Z(IZ-l) )+PT(IZ-l) 
CALL INSERT (ZIf,I2) 

REGION = 2 
GO TO 10 

C SOL 1D-L I GUI D/LIQUID INTERFACE 
2' uLLTZ = -CIH20S HZ) /DC3SDZ 
H(IZ) = H(IZ)+DhDZ*dELTZ 
CIHZOS(IZ) = 0. 

CIH20LUZ) = CIH20L ( IZ) +DC3LDZ* DEL 1 ) Z 
ZSLL(PASSNO) = Z(IZ)+DELTZ 
ZIF = ZSLL(PASSNO) 

PT t IZ) = (ZIF«Z(IZ-1))*(PTUZ)-HTUZ-U )/ (Z ( U) -Z UZ-1) )4pT(IZ-l) 
CALL INSERT (ZiFfIZ) 

REGION = 3 
GO TO 10 

C liquio/liouid-vapor INTERFACE 
3 DELTZ = (TL-T ( IZ) ) /UTDZ 
T(I2) = TL 

H(IZ) = H<IZ>+QHD2*0ELT2 
HL b H{ IZ) 

CIHZOL ( IZ ) = CIH20LCIZ)+DC3LDZ*DELTZ 
ZLLV(PASSNO) = Z(IZ)+OELTZ 
ZIF = ZLLV(PASSnO) 

PT ( IZ) = UIF-ZaZ-l))*(PT(IZ)-PT(IZ-l))/(ZUZ)-Z£I2-l) )+PT(IZ-l) 
CALL INSERT (ZIF*IZ) 

REGION = 4 
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60 TO 10 

c uquid-vapor/vapor jnterface 

4 MLTZ a -CIH20L(iZ)/DC3LDZ 

JH1IZL.= H(XZ)+ DHD2 *PELTZ 

C|H20L(iZ) a 0. 

CIH20VUZ) s C IH20V ( IZ) +DC3VDZ*DELTZ 
ZLVV(PASSNO) s Z{IZ)+DELTZ 

Z_IF a ZLVV CPASSNO) „ 

pfaz) = (ziF-z(xz-i) )*{pt(iz)-pt(iz°i) )/(zaz)-Z(iz-i))+PT(iz-i) 

.CALU_XN§£RX-tlI£ilZl 

REGION = 5 
6.0 TO 10 

C END OF REACTOR BED 

5 DELTZ a ZEND-Z(IZ) 

H(IZ) = H(IZ)+DHDZ*DELTZ 
CIH20V(IZ) = CIH20V(IZ)+DC3VDZ*DELTZ 
Z(IZ) = ZEND 

ZIF = ZEND 

PT(XZ) = (ZIF-Z(IZ-l) )*(PT (IZ)-PT (IZ-1) )/(Z( IZ)-Z(1Z*»1) l+P.TUZ-l) 
C ADJUST TEMP# PRESSURE# CONCENTRATIONS 
10 CIH2(IZ) = CIH2(IZ)+DC1DZ*DELTZ 
CXOji(U) = C.I.02. ( |.Z) +DC2DZ*REUTZ 
CIH20CIZ) = CIH20S(IZ)+CIH20L(IZ)+CIH20V(IZ) 

CIHE(IZ) = CIHECIZ)+DC4DZ*DELTZ 
IF (REGION, EQ ,2*0R »REG ION, EQ »4) 60 TO 12 
T(IZ) = T ( IZ} +DTDZ+DELTZ 
12 P(1Z) = P(IZ)+QPDZ*DELTZ 

RETURN ' 

END 
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n r n , o n rv 


FUNCTION FMAX (X,N,MAX»M) 


ROUTINE TO FIND MAXIMUM' VALUE OF ALL ELEMENTS IN ARRAY X 

DIMENSION X {MAX f M) 

FMAX.. = .XL1j.1). . 

DO 10 J=1»M 
DO 10 1=1 rN 

IF (X ( I » J) -FMAX) 10 f 10 f 5 
FMAX = X(I»J) 

CONTINUE 

RETURN.. 

END 


FUNCTION FMlN (X,N f MAXrM) 

ROUTINE TO FIND. MINIMUM VALUE OF ALL ELEMENTS IN ARRAY X 

DIMENSION X(MAXrM) 

FrtlN = X<1.1) 

UO 10 J=1»M 
DO 10 1=1 »N 

IF (FMJN-X ( I »U) } 10»1Q»5 . 

5 FMIN = X(I»J) 

10 continue 
RET URN 
END 
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XXXXXX XXX 


SUBROUTINE RDPNCH < T IME * HW * MINT ) 

C ... 

C ROUTINE TO READ AND PUNCH AUXILIARY INPUT 
C 

INTEGER PUNCH 

- PARAMETER.. NAX-.7JU. NXIsftQ 

PARAMETER NPP=30* NRTS=200 

COMMON /CB1/ TPP { NPP * NAX ) » C02 (NPP » NAX ) * CH2 C NPP » NAX ) , CH20 { NPP i NAX } , 
T (NAX) »P(NAX) ,PT (NAX) ,H(NAX) *CIH2(NAX) ,CI02(NAX) *CIH20(NAX) , 
CIH20S (NAX) ,CIH20L(NAX) »CIH20V(NAX) »CIHE(NAX) ,DELT (NTI) » 
PPT(NAX) 

COMMON /CB2/ Z(NAXi.».ZSSUNTl)jZSLLlNm iZLLV.lNTll ,ZLU.VXNTX) 

COMMON /CB3/ ALPHA * GAMMA * KAPPA * HC * KP , A , DELTT * KCH2 , KC02 , KCH20 * TCH2 * 
TCG2,TCH20*PCH2*PC02»PCH2a»MH2fM02,KH20,MHE,HR»HA,HAl*HA2, 

AW , C W * MW * DC'* TA * D0H2 * D0O2 , DOH20 * ALPHAK * AP » DELTA * G , MU , HS * HL » 

HV * TW * RHOM * WH2 , W02* WH20S , tf H20L * WH20V * WHE # CFBAR * MBAR . C iCPi » 
CICP2*CICP3,RH0,RH0V,DELHC#TITPS*DTDZ*DHDZ,TS,TL*DPDZ,DC1DZ* 
DC2DZ * D.C.3SD.Z.I DC3LDZ , DC.3VDZ ».DC4DZ» RHOP , CP , DELTP .DR.TRS tEDIlMEL. 
PF * TSS * AC * VC * PREST * GSS * MBARSS 
COMMON /CBV MOFZ » NPASS , NAX I AL * INPERP , MODPC , MODNO * PUNCH , I TSTEP * 

X JDP * ISHUT * IHOT 

DIMENSION HWU) ,C02T(1) ,CH2T(1) ,CH20T(1) »PTT(1) ,TTl(i) ,ZT(1) * 

X MINT (I) 

REAU (5*1) (ZU) ,1=1* NAXIAL) 

READ (5*1) (PTCl) *I=i*NAXIAL) 

READ (5*1) (T ( I ) * 1=2 *NAXIAL) 

REAU (5*1) <H(I) *I=2*NAXIAL) 

REAU (5*1) (P(I) *I=2*NAXIAL) 

READ (5*i-> (CXH2 ( I ) » I =2* NAXIAL) 

READ (5*1) (CI02(I).*I=2*NAXIAL) 

REAU (5,1) (CIH20V( I) ,1=2, NAXIAL) 

READ (5*1) (CIHE(I) *I=2*NAXIAL) 

READ (5,1) (PPT ( I ) * 1=1 , N AXIAL) 

READ (5,1) ( (TPP(I*J) *I=l,NPP),J=l f fjAXIAL) 

READ (5*1) ( (C02(I*U> *1=1, NPP) *J=1*NAXIAL) 

READ .(5,1) ( (CH2 ( I *J) * I=1»NRP) *J=1 * NAXIAL) 

READ (5*2) (MINT ( I ) » 1=1 *NAXI AL) 

READ (5*1) TIME, TW, DELTT 
1 FORMAT (8E10.5) 

RETURN 

ENTRY RPNCH2 (T IME »NAXT *ZT *PTT *C02T *CH2T *CH20T * TTI ) 

READ (5*2) NAXT 

READ (5*1) ( ZT ( X ) * 1=1 » NAXT) 

READ (5*1) (PTT(I) *1=1, NAXT) 

READ (5*1) (C02T(I) *I=1*NAXT) 

HEAD (5*1) (CH2T(I) *1=1, NAXT) 

•READ (5,1) (CH20TU) *1=1, NAXT) 

READ (5,1) (TTKI), 1=1, NAXT) 

READ (5*1) TIME *TW 
RETURN 

ENTRY RPNCH3(TIME*HW,MINT) 

10 WRITE (7,1) <ZU)*I=1*NAXIAL) 

WRITE (7,1) (PT ( I ) * 1=1 ,NAXIAL) 

WRITE (7,1) (Til) *I=2*NAXIAL) 

WRITE (7*1) (H ( I ) * I=2»NAXIAL ) 

WRITE (7*1) (PlI) *I=2»NAXIAL) 

WRITE (7*1) (CIH2(1) *I=2*NAXIAL) 

WRITE (7,1) (CI02( I) * I =2, NAX I AL) 
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WRITE <?,1> (CIH20VU), 1=2, NAXIAL) 

WRITE (7,1) (CIHE ( 1 ) , I =2, N AXIAL) 

WRITE (7,1) (PPT< I), 1=1, NAXIAL) 

WRITE (7,1) ( (TPP ( I ,D) , 1=1 ,NPP) ,J=1, NAXIAL) 
WRTTE (7,1) ( tC02(l,J) ,I=1,NPP) ,J=1, NAXIAL) 
WRITE (7,1) ( (CH2(I»U) ,I=1»NPP) ,J=1, NAXIAL) 
WRITE (7,2) (MINT(I) ,1=1, NAXIAL) 

WRITE (7,1) TIMt»TW,DELTT 
2 FORMAT (4012) 

RE_TUR_N _ 

ENTRY RPNCH4 (TIME) 

WRITE (7,2) NAXIAL 

WRITE (7,1) (Z(I), 1=1, NAXIAL) . 

WRITE (7,1) (PT(I)» 1=1, NAXIAL) 

WRITE (7,1) (CI02(I) ,1=1, NAXIAL) 

WRITE (7,1) , ( Cl H2 (I ),.I=1, NAXJ A.U. 

WRITE (7,1) (CIH20(I), 1=1, NAXIAL) 

WRITE (7,1) (T(I) ,1=1, NAXIAL) 

WRITE (7,1) T3ME,TW 

RETURN 

END 
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. r, n .n 


SUBROUTINE OUTPUT "<TM,TT# TP i PPV f'H27fo 2 r FH20 ,FHE,N1 ,'n 2» f0CT) 

ROUTINE TO PRINT TEMP, PRESSURE,' MOLE FRACTION PROFILES‘Tn^UMMaW?ORM 

’ PARAMETER" NAX=70T NTI=4o' 

INTEGER FORMT (9) / 54H ( UX , 1H* , 5X , 3HSEC , { 9X ,.3HSEC ) /U 

XX, 1H*) / 

COMMON /CB2/ Z^AX},2SSL{NTl),^SLL(NTI),ZLLV(NTl),ZLVy{NT.U 

COMMON /CB4/ NOFZ ,NPASS ,N AXIAL , INPERP ,MoDPC , MOONO , PUNCH , ITSTEP , 

X . UQP* I SHUT j I HOT _ _ __ 

DIMENSION NUMBER (9) 

DIMENSION TM(1) ,TTU) ,TP(D ,PP( 1} ,FH2(1) ,F02(1) ,FH20(1) »FH£(1) 

DATA NUMBER / ’ 12345678 

X 9' / 

C NO. OF COLUMNS OF OUTPUT = (N2-N1+1) 

F.QRMT15) = NUMBER (N.2-N1) 

WRITE (o,10Q) 

100 FORMAT {*!» / 36X 'TIME PROFILES OF INTERSTITIAL TEMPERATURE VS AX 
XIAL POSITION* // 3X ’AXIAL *’ 42X ’TEMPERATURE {DEGREES R)’ / 

X ’ POSITION *’ / 3X '(FT) *♦ ) 

WRITE (6,101) (TM(I) »I=N1,N2) 

101 FORMAT ( UX »* * FG.4, 9(4X,F.8.4) .) . _ 

WRITE (fa, FORMT) 

Kl = 1 

K2 s MODPC 

DO 10 I=1,NAXIAL 

WRITE (6,102) 2(1) , CTT(K) »K=K1,K2) 

Kl = K2tl 

10 K2 = (I+D+MODPC 

102 FORMAT ( £9,4, ’ EU.5, 9(E12.5> ) 

WRITE (6,109) 

109 FORMAT t ’ 1 • / 38X ’TIME PROFILES OF PARTICLE TEMPERATURE VS AXIAL 
XPOSITION* // 3X ’AXIAL ** 43X ’ TEMPERATURE (DEGREES R)* / 

X » POSITION *' / 3X ’(FT) *’ ) 

WRITE (6,101) (TMU),I=N1»N2) 

WRITE (6, FORMT) 

Kl - 1 

K2 = MODPC 

DO 70 I=1,NAXIAL 

WRITE (6,102) 2(1), (TP(K) ,K=K1,K2) 

Kl = K2+1 

70 K2 = (I+l)*MODPC 
WRITE 16,103) 

103 FORMAT (*1» / 40X ’TIME PROFILES OF CHAMBER PRESSURE VS AXIAL POSI 
XTION’ // 3X ’AXIAL ** 46X 'PRESSURE (PSIA) ' / 

X * POSITION *• / 3X ’(FT) *• ) 

WRITE (6,101) ITM(I) ,I=N1,N2) 

WRITE (6, FORMT) 

Kl = 1 

K2 = MODPC 

DO 20 I=1,NAXIAL 

WRITE (6,102) 2(1), (PP(K) ,K=KlrK2) 

Kl - K2+1 

20 K2 = (I+l)*MODPC 
WRITE (6,104) 

104 FORMAT (»!’ / 40X 'TIME PROFILES OF H2 MOLE FRACTION VS AXIAL POSI 
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XTlON’ // 3X * AXIAL ** 38X »H2 MOLE FRACTION (DIMENSIONLESS) » / 

X ’ POSITION *’ / 3X . *_(F_T). *♦ X 

WRITE (6,101) (TM(I),I=N1,N2) 

WHITE, JLfajJr.QRMIl 

K1 = 1 

K2 s MODPC 

DO 30 I=1,NAXIAL 

WRITE (6*1.08) 2.(11* <FH2(K)jK=K:uK21 . ... 

K1 s K2+1 

30. .., .A£..a-II±ngMflQgC. 

WRITE (6,107) 

107 FORMAT (’1* / 4QX ’TIME PROFILES OF Q2 MOLE FRACTION VS. AXIAL POS1 
XTION' // 3X 'AXIAL *♦ 38X *02 .MOLE FRACTION (DIMENSIONLESS) ♦ / 

X . * POSITION *' / 3X ’(FT) *’ ) 

WRITE (6,101) (TM(I) ,I=N1,N2) 

WRITE (6*F.0RMT)’ . 

K1 s 1 

K2 s MODPC 

DO 40 I=1»NAXIAL 

WRITE (6,102) 2(1), (F02 (K) ,K=K1,K2) 

K1 = K2+1 

40 . K2 = (ItU.2FM0.DPL. 

WRITE (6,105) 

106 FORMAT (’!’ / 40X ’TIME PROFILES OF H20 MOLE FRACTION VS AXIAL POS 
X I T I ON * // 3X 'AXIAL *• 38X ’H20 MOLE FRACTION (DIMENSIONLESS)’ / 
X ♦ POSITION *’ / 3X ’(FT) *' ) 

WRITE (6,101) (TM(I),I=N1,N2) 

WRITE (6,FORMT) 

K1 = 1 

K2 = MODPC 

DO 60 I=1»NAXIAL 

WHITE (6,102) 2(1), (FH20 (K) »K=K1 »K2) 

K1 = K2+1 

50 K2 = (I+.D+MQDPC 
WHITE (o,106) 

106 FORMAT (»1» / 40X ’TIME PROFILES OF HE MOLE FRACTION J/S AXIAL POSI 
XTION* // 3X * AXIAL *» 38X 'HE MOLE FRACTION (DIMENSIONLESS)’ / 

X ’ POSITION *’ / 2X ’(FT) *♦ ) 

WRITE (o,101) (TM(I) »I=N1,N2) 

WRITE (6,F0KMT) 

K1 = 1 

K2 = MODPC 

DO 60 I=1»NAXIAL 

WRITE (6,102) 2(1), (FHE(K) ,K=K1,K2) 

K1 - K2+1 

60 K2 = ( X+l) *MODPv 
IOC( = IOCT+1 

RETURN 

END 
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non 


SUBROUTINE SETUP ( TM »MFH2 r MF02 » MFH20 > MFHE » PASSNO ? K i KM ) 
ROUTINE TO SET UP ARRAYS FOR PRINTER PLOTS 


C 


L 


b 


io 


INTEGER PASSNO 

REAL,. MFHE.MF.OEi MFH2Q tMFilE 

PARAMETER NAX=7q» NTI=40 
PARAMETER NPP=3Q, NRTSS200 

COMMON /CB1/ TPP ( NPP * NAX ) * C02 ( NPP » NAX ) * CH2 ( NPP t NAX ) ♦ CH20 ( NPP t NAX ) , 
X . T (NAX) »P (NAX) #PT (NAX) »H(NAX) ,CIH2(NAX) ,CI02(NAX) #CIH20(NAX) # 
X CIH20S (WAX) »CIH20L(NAX) »CIH20V<NAX) .CIHE(NAX) »OELT ( NT I ) r 

X PPT(NAX) .... ... .... 

COMMON /CB2/ Z(NAX) >ZSSL(NTl) »2SLL(NTI) »ZLLV(NTI) rZLVV(NTl) 

COMMON /CB4/ N0F2»NPASS#NAXIAL* INPERP »MODPC f MODNO»PUNCH# ITSTEP# 

X JDP» ISHUT t IHOT 

DIMENSION LABEL (91) »LBL(15) ,LABLY(91) ,LABLX(24) ,TZ{5) 

DIMENSION Z1(NAX»5) ,Yl(NAX r 5) »Y2(NAX»5) ,Y3(NAX»5) ,Y4(NAX»'5) r 
X Y5(NAXi5) »Y6 (NAXiS)*Y7(NAXi5) 

DIMENSION TM(1) f MFH2(1) ,MF02(I) fMFH20(l) »MFHE(1) 

LABELS FOR PRINTER PLOTS ARE IN FOLLOWING DATA STATEMENTS ... 

DATA LA0EL/b*» *»»TIME PROFILES*. INTERSTITIAL TEMPERATURE VS 

X AXIAL POSITION »»6*» 1 1 * TIME PROFILES: CHAMBER PRESSURE VS 

XAXIAL PQ5ITION » »6*» •■•TIME PROFILES: H2 MOLE FRACTION V 

XS AXIAL POSITION «» 6** PROFILES! 02 MOLE FRACTION 

X VS AXIAL POSITION »,6*« «,«TIME PROFILES! H20 MOLE FRACT 

XION VS AXIAL POSITION »»6*» »»*TIME PROFILES! HE MOLE FRAC 

XT I ON VS AXIAL POSITION */ 

DATA LBL/6*' 1 r ’TIME PROFILES: PARTICLE TEMPERATURE VS AXIAL 


XP0S1TI0N ’/ 

DATA LABLX/* AXIAL POSITION (FT) *** TIME (SECONDS) »» 
X* TIMES (IN SEC) REPRESENTED BY CURVES 1 THRU '» 'AXIAL POSITION 
XS (IN FT) REPRESENTED BY CURVES 1-*/ 

DATA LABLY/' •M»,tp»^»E»f»R*r »A»»»T'r»Ul.»_'R»»»E»,* » 


X t 


,» i,» «,» ♦ («. *D*f 'E» » »G» » 'R* r ♦£♦ > »E* » 'S' » * «»»«'»»)'»< '» 


X » *»*M'r*0*,»L*»'E'»' *#*F«» , R*» , A*r , C»# , T , »»I , »*OS , N'»* 

X * (»,»D , » f I*i »M»»*E» PN'PS'iN'i'O' PN'.'L’PE-S'S'.'S' »*)♦»♦ 
X ♦ , 'P* r *k* , »E» * 'S' , »S* t *U* , 'R» , 'E* » * ♦,» »»• *.» '»♦ *#» 

X »»♦ »»»(*»*P*» , S , f »I'f »A'f»)S* '»» *»' '»• »/ 

IF (K.GT.l) GO TO 5 
MAX s NAX 

CALCULATE NUMBER OF POINTS PER PLOT 
NPTS = NAXIAL/MuDNO 

IF (MOD(NAXIALrMODNO) ,NE.O> NPTS=NpTS+l 


1 = 1 

DO 10 J=1 fNAXIAL rMODNO 
Zl(i»K) = Z(U) 

Yl(i.K) = T ( J) 

Y7 ( I f K) = PT(J) 

Y2U»K) = P(J) 

Y3( I »K) = MFH2 ( J) 
Y4U.K) = MF02(J) 

Y5 ( I »K ) = MFH20(J) 

Yb ( I >K ) = MFHE(J) 

I = 1+1 

Z1 (NPTS t K ) = Z(NAXIAL) 
Y1(NPTS»K) s T(NAXIAL) 
Y7(NPTS»K) = PT(NAXIAL) 
Y2(NPTS»K) s P(NAXIAL) 
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Y3(NPTS»K) = MFH2(NAXIAL) 

Y4(NPTS,K> = MF02 (NAXIAL) 

Y5(NPTS»K) = MFH20 (NAXIAL) 

Y6 ( NPTS »K)= MFH E ( NAXIAL ) _ 

K s k+I 

IF ( P ASSNO • NE # NPASS ) 60 TO 15 
KM 5 K-l 

15 IF (K.LE.KM) RETURN 

ITM = PASSNO-KM 

00 20 ?=1»KM 

20 TZ(I) = TMUTM+I) 

C CALL PRINTER PLOT ROUTINE WHEN ALL Yl, Y2,...»Y6 ARRAYS FULL 

CALL PRPLOT (Zl,Yl»NPTS,KM,MAX,TZ»LABEL(i) ,LABLX(1) ,LABLX(9) , 

X LABLY (1) , LABLY ( 17) ,16*0) 

CALL PRPLOT (Z1,Y7,NPTS,KM,MAX,TZ»LBL(1),LABLX(1),LABLX(9), 

•X LABLY(l) ,LABLY(17),15,0) 

CALL PRPLOT < Z1 » Y2, NPTS, KM, MAX, TZ, LABEL (17) *LABLX(1) ,LABLX'(9) , 
X LABLY (62) , LABLY (77) ,15,0) 

CALL PRPLOT (Zi ,Y3» NPTS, KM, MAX, TZ, LABEL (32) » LABtX ( 1 ) ,LABLX (9) , 
X LABLY{32), LABLY (47), 15,0) 

CALL PRPLOT (Z1,Y4,NPTS,KM,MAX,TZ*LABEL(47) ,LABLX(1) ,LABLX(9) , 
X LABLY (32), LABLY (47), 15, Q) 

CALL PRPLOT (Zl ,Y5> NPTS, KM, MAX ,TZ» LABEL (62) ,LABLX(1) ,LABLX{9) , 
X LABLY(32) ,LABLY(47) ,15,0) 

CALL PRPLOT (Zl,Y6,NPTS,KM,MAX,TZrLABEL(77) ,LABLX(1) ,LABLX(9) , 
X LABLY (32) * LABLY (47) , 15, 0) 

K = 1 

return 

end 
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noon 


SUBROUTINE PRPLQT (XrY,NP»NV»MAX»T2,LBL»LXl,LX2,LYl,LY2,NL>NC) 

ROUTINE TO DEVELOP PR INTER PLOTS OF Y VS X 
MAXIMUM OF FIVE SUPERIMPOSED .CURVES PER PLOT 



DIMENSION LBL(l) *LX1(1) f LX2(l) #LYl£l) #LY2U) 

DIMENSION LINE £126) fXGRAD(12) »YGRAD(12) » XSCAL £ 12 ), YSCAL (12J 
INTEGER NUMBR(IO) /*1» r *2» r *3» r *4» * *5* t *6» t *7* » »8» , »9N »0»/ 

INTEGER ALPHA £26} /*A» t »B* t *C * , «D* » *E* * «F» * »G' t 'H* , » I » » * J» * »K»* 

X ’LS«M»»*N»» , 0»»»P*»*Q’»*R»»*S»#*T*»*U»f»V«, 'W')'X*i'Y* f *Z»/ 

DATA KBLNK/A . .*/ ^ IRLUS/ '.t.*/^.MlNUSy « - ' J t, J£ QLN/ U •/. . 

WRITE (6» 1000) (LBL£I) »Isl»NL) 

1000 FORMAT £»1* // 20A6) 

NSCAL s 120 
LK = 12 
NTRS = 1 
. NNN z, I 
XYMRGN =0, 

XMAX = FMAX £ X r NP » MAX » NV ) 

XMlN = FMIN (X »NP » MAX t NV ) 

C SCALE AXES FROM MAX AND MIN VALUEh (SAME SCALING PROCEDURE FOR BOTH AXES) 
C LXSF! POWER OF 10 TO WHICH X VALUES ARE SCALED 

L LYSF? POWER OF 10 TO WHICH Y VALUES ARE SCALED 

920 IF (ABS(XMAX) *LT.1» ) GO TO 100 
DO 30 11=1 i 30 

IF (ABS(XMAX),LT.10,**II) GO TO 31 

30 CONTINUE 

31 LYSF = -(II-l) 

GO TO 22 

100 DO 40 II=1»10 

IF (ABS(XMAX) .01 ,10.**(-II) ) 60 TO 41 

40 CONTINUE 

41 LYSF = II 

22 MANT = IFIX(XMAX*lO.**(LYSF+2) ) 

DO 302 K=lf.lQ 

IF £MOD(MANT+K»5) •EQ«0) GO TO 3t;3 

302 CONTINUE 

303 AMM = FLOAT (MANT+K 3/100, 

XSPAN = amm-xmin*io,**lysf 

MANT = IFIX ( XSPAN*100 , ) 

• DO 304 K=lr NSCAL 

IF (MOD£MANT+K»NSCAL) ,EQ,0) GO TO 306 

304 CONTINUE 

30b XINCR s { FLOAT (MANT+K)/100.) /FLOAT (lK-1) 

HALFX = XINCR/2, 

IF £(XMIN-(AMM-FL0AT(LK-1)*XINCR)),GT. HALFX) XYMRGNsHAlFX 
C ADD .HALF A SCALING INCREMENT TO FORM MARGIN IF NEEDED 
YSCAULK) = AMM+XYMRGN 
DO 306 I=LKf2»-i 

306 YSCAL £1-1) = YSCAL ( I ) “XINCR 

IF ( ABS £ YSCAL( I) ) ,LT»10«) GO TO 338 
LYSF = LYSF-1 
XINCR = XINCR/10. 

DO 339 l=ltLK 

339 YSCAL ( I ) = YSCAL ( I ) /10 . 

C SCALING COMPLETE ON X-AXIS IF NTRS=0» SCALING COMPLETE ON Y-AXIS ALSO 

338 IF lNTRS.EQ.0) GO TO 42 
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c xscalj scaling 'gradations on x-axis""'” 

c yscal: scaling gradations on_ y-axis 

LXSF = LYSF 

XINCX ~__XINCR 

'■"DO" 43 UJ=f»LK 

XGRAD(JU) = YSCAL(JJ) 

43 XSCAL(JJ) = YSCAL (UU) 

C RESET SCALING PARAMETERS TO SCALE FOR Y-AXIS 
NSCAL = 50 
LK = 6 _ __ 

NTRS = 0 
MM s LK 
XYMRGN =0, 

XMAX = FM AX ( Y , NP , MAX , NV ) 

XMlN = FMI N ( Y , NP , MAX , NV ) 

RANGE = ABS ( XMAX-XMIN ) 

C IF Y-AXIS RANGE LESS THAN 1 PERCENT OF MAX VALUE, ASSUME CURVE IS CONSTANT 
IF (RANGE. LT.ABS(,01*XMAX)) XNAX=XMIN 
GO TO 920 
42 DO 44 U=i»LK 

44 YGRAD(J) = YSCAL (J) *10.** (-LYSF) 

C GENERATE TOP AND BOTTOM LINES OF GRAPH 
9 DO 13 K=16,126 

IF (MOD (K-lfa , 10) ) 14,15,14 

14 LINE(K) r MINUS 
GO TO 13 

15 LINE(K) = IPLUS 
13 CONTINUE 

NMN = 1 1-NNN*5 

WRITE (6,52) YGRAD(NMN) , (UNE(M) ,M=lb,l26) 

52 FORMAT (6X, E10.5, IX, 111A1) 

IF (NNN.EQ.2) GO TO 712 

c generate internal lines of graph 
C = XINCR*10.**(-LYSF)/10, 

A s YSCAL (6) *10 .** (-LYSF) +C/2 . 

U = A-C 

C UEGIN LOOP TO SET UP EACH INTERNAL LINE OF GRAPH — LINES ARE LABELED 
C 1 THRU 49 FROM TOP TO BOTTOM 

DO 401 JS=1»49 
DO 403 K=l,126 
403 LlNt(K) = KBLNK 

IF (MOD( JS,10) ) 18,19,18 

18 LINE (lb) = ICOLN 
LINE (126) = ICOLN 
60 TO 20 

19 LINEU6) = IPLUS 
LINE (126 ) = IPLUS 

C A, 8 ARE UPPER AND LOWER BRACKETING VALUES FOR Y-COORDINATE OF POINT 

C C IS SCALED HEIGHT OF PRINT CHARACTER (A, B ARE ADJUSTED BY THIS AMOUNT 

C ON EACH NEW PRINT LINE) 

20 A - A-C 
B - B-C 

DO 427 NN=i»NV 
DO 427 JR=1,NP 

C GET SUBSCRIPT OF UR Til POINT OF NN TH CURVE 
UD = JR+ (NN-1 ) *MAX 

C TEST IF Y-COORDlNATE LIES WITHIN RANGE OF A MINUS B OF LINE NO. JS 
IF (Y(JD).GT.A.OR.Y(JD),LE.B) GO TO 427 
C VALUE IS POSITIONING VALUE FOR X-COOROINATE OF POINT 
VALUE = X(JQ)*10.**LXSF 
DD = XINCX 

NPOS = (VALUE-XSCAL(i))/DD*10.U6.5 
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C TEST IF X-COOROINATE LIES WITHIN CHARACTER POSITIONS 17 THRU 125 OF LINE JS 

IF ( NPQS « LT • 17 *OR j NP0.S.,0T.»-12 5 30. „I.0_JNSL7. . . . 

C CURVES IDENTIFIED WITH DIGITS IF NC=Q 

Ll NE(NFQSl . = NUMBR1NN) 

URS = MOD (JR >26) 

IF (JRS.^£a,.0) JRS=2 «l . . 

C CURVES IDENTIFIED WITH SEQUENTIAL CHARACTERS OF ALPHABET IF NC=1 

IF (NC ,EQ.i) L I NE(NP.OS) =ALPHA.(JRS ) 

427 CONTINUE 

. . IF (MODUSUfl U- 25 j<g4,£5 

24 MM = MM-1 

WRITE (6,52) YGRAD(MM) , (LINE(M) »M=16#126) 

GO TO 401 

25 IF ( JS . NE . 25 . ANp , JS , NE » 26 ) GO TO 36 
IF (JS.EQ.26) GO TO 27 

DO. 33 1=1 i 13 

33 LINE(I) = LY1(I) 

GO TO 36 

27 DO 34 1=1,15 

34 LINE(I) = LY2( I) 

36 WRITE (6,55) (LINE(M) ,M=1,126) 

55 FORMAT ( IX t 15A1» 1,X f 1UA1 1 
401 CONTINUE 

NNN = NNN+1 
GO TO 9 

C PRINT X-AXIS GRADATION MARKINGS 
712 WRITE (6,53) (XGRAD(M) ,KBLNK ,M=1 , 11 ) , XGRAD ( 12) 

53 FORMAT (14X# 111F6„3» A4 ) t F6.3) 

WRITE (6,56) (LXKI) ,1=1,4) 

56 FORMAT (»0* 56X, 4Ab ) 

WRITE (6,58) (LX2(I) *1=1,8) ,NV, (TZ(I) ,I=1,NV) 

58 FORMAT (»+* 81X, 8A6, 11 / 86X, 5F8.U) 

IF (LXSF.EQ.O) GO TO 57 
WRITE (6,55) LXSF 

55 FORMAT ( *0* 5QX *X-AXIS SCALING: X VALUES*10** * , 12 ) 

37 RETURN 

END 
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o n o 


suBroOt inTTmPRnT Tt M ,twt,ni,n2) 

ROUTINE TO PRINT WALL TEMPERATURE AND AXIAL INTERFACE VS TIME PROFILES 
PARAMETER NAX=70» NT 1=40 

COMMON /C82/ Z(NAX) , ZSSL(NT I ) oZsLL(NTI) ,ZL LV (NTI) , ZL W(NTI) 

dimension g t i > , p fi ) P tm < i ) , twt a ) 

WRITE (6,10) 

10 FORMAT (*1* //// 43X »WALL TEMPERATURE, INTERFACE VARIATIONS VS TI 
XME* // 15X 'TIME * WALL TEMP SOLID/SOLID-LIQUID SOLID 
X-LIQUID/LIQUIO LIQUID/LXQUID-VAPOR LlQUXD-VAPOR/VAPOR* / 

X • 15X SEC.) . . * (DEGREES. R) ». iftX „• CEIJ V 3U0.X JJFXllL_-ZL.-g5-X 

X »*»> 

DO 30 I=N1,N2 

30 WRITE (6,40) TM ( I ) » TWT ( I ) , ZSSL ( I } , ZSLL ( I ) , ZLLV ( I ) , ZLVV ( I ) 

40 FORMAT ( 13X, F3.4, * * E12.5, 5X» El2.5f 3(lQX» £12.5) ) 

. RETURN 

ENTRY TMPRP (TMrGfP.rNI *N2) 

WRITE (6,20) 

20 FORMAT ('O' // 44X 'PRESSURE, MASS FLOW RATE VARIATIONS VS TIME' / 

X / 43X 'TIME * PRESSURE* 8X ’'MASS FLOW RATE' / 43X 

X '(SEC) * (PSIA)' 9X ' (LB/FT**2-SEC) • / 51X •*» ) 

DO SO I=N1,N2 

50 WRITE (6,60) TMCI > r P ( I ) »&.( I) 

60 FORMAT ( 41X, F8.4, 2X, ♦*', E12.5, 7X» £12.5 ) 

RETURN 

END 


96 



80 

til 

90 

91 

92 

93 

94 

95 

9o 

97 

98 

99 


SUBROUTINE DEBUGS l TIME* IT) 
REAL MH2iM02 j MH20j-MH£iMU 
X MU2VST »MU3VST ,MU4VST 

HAB A METER. _ N A X= 7.Q. t _NT.I=4i}_. . 


PARAMETER NPP=30 

X . ' T(NAX)7P(NAX)'»PT(NAX) ,H(NaX) »CIH2 (NAX) »CI02(NAX)#CIH20(NAX) , 
X _ CIH20S (NAX) »CIH20U(NAX) , C.JH2Q V.INAX )_t_C IHE ( MAX ) » DELTXNT.I).,. 

X PPT(NAX) 

COMMON /CB2/ Z(NAX) jZSSLJNTl) ,ZSLk(NTI) »2LUV(NTI) »ZLVV(NT.I1 

COMMON /CB3/ ALPHA , GAMMA » KAPPA , HC » KP » A , DELTT » KCH2 , KC02 , KCH20 , TCH2 , 


X .. . . JC02jT,CH2Q.i.P.CH2li£C02,i;PXH2Q.»MH2-f.MQ2<J4Jd2:0.».MH£j-HRj±jAj±lAlAHA2.f 

X AW»CW,MW,DC»TA»DOH2»0002,DOH20,ALPHAK,AP,DELTA,G»MU»HS»HL, 

X . H V , T W » RHOM , WH2 , WQ2 » WH20S ,4slH2GL» WH20 V.» WHE , CFBAR « MB AR » CICP.l , ... 

X ClCP2,ClCP3»RH0,RH0V ,DELHC,TITPS,DTDZ»DHDZ,TS»TL,DPDZ,DC1DZ, 

X PC 2DZ , DC3SDZ f DC 3LDZ , D.C3VD2 , 0C4DZ , RHOP , CP , DELTP . DPTPS , SD.TIME * 

X . PF , TSS , AC , VC » PKEST , GSS , MBARSS 

COMMON. /CB4/. NOFZ f NPASS NAXIAL f JLNPERP.f.MODPC./ MQDNQ aEUNCHxLTSIEEil.^ 

X JDP , ISHUT , IHOT 

COMMON /TABLES/ AVSZ(40) ,APVSZ(4G> ,DELVSZ(4Q) ,HRVST(36) ,T.C1(22), 

X • TC2(20> ,TC3V(2G> ,TC4(18) ,MUlVST<34> »MU2VST(34) ,MU3VST<32) , 

X MU4VST ( 34 ) , CF1 VST { 34 ) , CF2VST (34 ) , CF3VST (34 ) , CF3SVT (20 ) » 

X . CF4VST(2Q) ,VPVST(60) ,TVSVP(48) »DHCVST{24> 

NAXIAT = NAXIAL 
NAXIAL = NOFZ 

WRITE (6,90) (CIH2(I) ,1=1, NAXIAL) 

WRITE (6,91) (CI02(I) ,1=1, NAXIAL) 


WHITE (6,92) (CIHHOS(I) ,1=1, NAXIAL) 

WRITE (6,93) (GIH20UI), 1=1, NAXIAL) 

WRITE (6,94) (CIH2Q.V (I), 1=1, NAXIAL) 

WRITE (6,95) (CIHE(I) ,1=1, NAXIAL) 

WRITE (6,96) (T(I)f 1=1, NAXIAL) 

WRITE (6,81) (PT ( I ) ,1=1, NAXIAL) 

WRITE (6,97) (H(I) ,1=1, NAXIAL) 

WRITE (6,98) (P(I) ,1=1, NAXIAL) 

WHITE (6,80). (ZTI) 1 1=1, NAXIAL) 

WRITE (6,99) TIMEtTw - 

FORMAT (’0* 57X ’ AXIAL POSITIONS* / ( IX , 10E13, 5 ) ) 

FORMAT l'0» SOX ’PARTICLE SURFACE TEMPERATURES’ / ( IX ,1 0E13 . 5) ) 
FORMAT ( *0’//48x» INTERSTITIAL CONCENTRATIONS OF H2’ / (1X,10£13.5) 


X ) 

format (*o» 48x ’Interstitial concentrations of 02 * / (ix,ioei3 .5) 


x ) 

format (*o’ 47x ’Interstitial concentrations of h 2 q(S>» / ux, 

X 10E13.5) ) 

FORMAT (’O’ 47X ’INTERSTITIAL CONCENTRATIONS OF H20(L)’ / (IX, 

X . 10E13.5) ) 

FORMAT (’O’ 47X ’INTERSTITIAL CONCENTRATIONS OF H20(V.)_» ./ IlX, _ - 
X 10E13.5) ) 

FORMAT (’O’ 48X ’INTERSTITIAL CONCENTRATIONS OF HE’ / (1X,10E13.5) 


X ) 

FORMAT C’O* 52X ’INTERSTITIAL TEMPERATURES* / (1X,10E13.5) ) 
FORMAT (’0* 60X ’ENTHALPIES’ / (1X,10£13.5> ) 

FORMAT (’O' 53X ’INTERSTITIAL PRESSURES’ / (1X,10E13.5) ) . 

FORMAT (’O’ 22X ’TIME =*,E11.5, 3X ’WALL TEMP =’,E11.5) 

NAXIAL = NAXIAT 

RETURN 

END 
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;o o o r> o 


SUBROUTINE UNBAR (T, ikVXIN, YIN,ZZ,KK> 


0 


THIS ROUTINE WILL DO LINEAR, QUADRATIC, OR CUBIC INTERPOLATION' BETWEEN 

tabulated points in table t — interpolation may be either univariate 

OR BIVARIATE 


DIMENSION A { 6 ) ,X(6) ,Y(6f,TU> 
II = IK+1 
N = 3 
N2 a 2 

IP (T( II )-3» ) 700,701,702 
..7M .1F..1TAIU) „6Q.f7.pijm 

704 IF (T(in-2.) 705,706,701 ~ 

705 N a 1 

GO TO 707 

706 N a 2 

707 N2 a 1 

.701 II = lit 1. 

702 N1 = N+l 

DO 50 L=II,II 
IF (T (L) ) 60,60,51 
60 KK = -1 

zz - o, 

GO TO 9999 

51 NX = T(L) 

IF (T(L+1) ) 60,52,50 

52 NY a 0 
60 TO 53 

50 NY = TCL+1) 

53 CONTINUE 
KK = 0 
KY = 0 
XX a XIN 
YY = YIN 
J1 a I 1+2 

J2 = NX+II+i 

IF (XX-T ( Jl ) ) 301 » 306 ,400 
400 DO 302 J=J1,J2 

IF (XX-T(J)) 304,304,302 

302 CONTINUE 
309 KK = 2 

XX = T ( J2) 

308 JX1 = J2-N 
60 TO 305 
301 KK = 1 

XX = T(J1) 

306 JX1 = Jl 
60 TO 305 

304 IF (J-Jl-1) 301,306,307 

307 IF (J-J2) 303,308,309 

303 JX1 = J-N2 

305 CONTINUE 
XINT — XX 

IF (NY) 1500,1500,3000 
1500 DO 1599 L=i,Nl 
X (L) a T(JX1) 

LY = JX1+NX 
Y (L) = T (LY) 


10 

_ -.20 ___ 

30 

fid _ 

50 

6fi 

70 
.60, 
90 
10 Q. 
110 

120 

130 
14Q 
150 
160 . 
170 



190 
200 
210 
220 
230 
24Q..-. 
250 
260 
270 
280 
290 
30Q 
310 
320 
330 
340 
350 
360 - 
370 
380 
390 
400 
410 
420, . 
430 
440 
450 
460 
470 
480 
490 
500 
310 . 
520 
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1599 JXl = UXi+1 

. i=.i 

GO TO 54 
3Q QQ J1 = Jl+NX 
J2 = J2+NY 

. I F ... .( YY-TJ jJ 1.)_ 3lL#_ 3_16_# 40_1 
401 DO 312 J=Ji » J2 

. .IF. (YY-T(J) l 314,314,31 2 __ 

312 CONTINUE 

„3.12.__KX_=j6 

YY = T(J2) 

316 <J_Yl. = J2-N 
GO TO 315 

311 KY = 3 

YY s T(J1) 

3.16 31SJ.JS sJfJL 

GO TO 315 

314 IF .(J-Jl-1) 311,316 #317 

317 IF (J-J2) 313»318#319 

313 JY1 5 J~N2 

315 CONTINUE 
JX2 m JXl 

LY = JY1+NY ♦ ( JX2-I 1-1 ) 

LYl = LY 
DO 3099 L=l#Ni 
X(L) = T(JX2) 

Y (L) = T(LY1) 

LYl = .LY1+NY 
3099 JX2 = JX2+1 
1=0 
GO TO 54 
3098 YU) = ZZ 

DC 4400 1=1 #N 
1Y1 = LY+I 
Y (1+1) = 0. 

DO 4050 MM=l,Nl 
y ( i+i ) = Ya+i)+T(Lri)*x(MM) 
4050 LYl = LY1+NY 
4400 CONTINUE 

DO 4199 L=i#Nl 
X (L) = T(JY1) 

4199 JYl = JY1+1 
XINT = YY 
I = 1 

54 D = 1, 

X(N+2) =. X(l) 

XCN+3) = X<2) 

DO 55 J=l#Ni 
A (J+l) = X(J+1)-X(J) 

TPAL1 = XINT-X(U) 

IF (TPAL1) 57,56, 57 
58 a. - Y(J) 

X(l) = o. 

X(2) = 0. 

X (3) = 0. 

X (4 ) = 0. 

X(.J) = 1, 

GO TO 59 

57 D = D*TPAL1 

GO TO (711 #712 #713 ) t N 
711 X(J) = TPALl/A (J+l) 

GO TO 55 


530 

540 

550 

s $Q 

570 
5 8Q. . 
590 
600 
610 
-.62 9 . 
630 

. 6_4Q 

650 

66J5L _ 

670 

6fl.Q — 

690 

..JZOO 

710 

_7.2Q. _ 

730 

740 

750 

760 

770 

7.30. . 

790 

. 800 - 
810 
820 
830 
840 
850 
-6611.--. 
370 
680 
390 
900 
910 

- _92a... 

930 

940 

950 

960 

970 

380 — 

990 
1000 
1010 
1020 
.1030 
1 OML .. 
1050 
1060 
1070 
1080 
1090 
uacL. 
1110 
1120 
1130 
1140 
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712 X(J) s -TPALi 1150 

GO TO 55 1160- 

713 X(j> = ( X ( J+2 ) -X ( J ) ) +TPAL 1 ' " 1170 

55 CONTINUE 1180 

A(ll 5 A(N+2) 1190 

ZZ - 0, _ JL2CL0_ 

DO 56 J=1»N1 1210 

XjU) = D/<A(J>*A{J+1)*X<J> > 1220 

ZZ s ZZ+Y(J>*X<J) 1230 

5.6 CONTINUE ___ 1240 

59 IF (I) 3098 #3096 ,9999 1250 

9999 KK = KK+KY 1260 

RETURN 1270 

END 1280 
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n c. o r n or o r r o o h n n n c n coo > o r> 


BLOCK DATA 


REAL MUiySTiMU.SVST ,My3_VS.T,HU4_VST 

COMMON /TABLES/ AVSZ(40) , APVSZ140) »OELVSZ(40) »HRVST(36) ,TC1 (22) , 

X MU4VST (34) ,CF1VST(34) ,CF2VST(34) ,CF3VST(34) ,CF3SVT(20 ) t 

X CF4VST(20),VPVST<60)fTVSVP<48)rDHCVST(24) . . 


HEAT. OF REACTION OF 02 VS TEMPERATURE JPEG R) 


- DATA .(HRVST1 1 ) >.1=1 r 

36) 

/_o, t i.. 

t. .16 ». 1 

Q. 




temperatures 









X . .0. r 

180. 

, 

360. t 

720. 

I 

1080. 

, 1440. 

» 

X 1800, , 

2160. 

» 

2520. , 

2860. 

9 

3240. 

, 3600, 

9 

X 3960. f 

4320, 

, 

4680. , 

5040. 

t 




HEATS OF REACTION 








. . X.-. -J6423... i 

-6461. 

r 

t64.78^ . , 

-6529.. 

I 

-6581. 

» . —6626.. . 

* 

X -6665. » 

-6696 

, 

-6720 . , 

-6740. 

9 

-6755. 

, —6766. 

9 

X . -6777. , 

-6786. 

t 

-6794. t 

-6801. 

/ 





FOUR VISCOSITY TABLES FOLLOW ... 


H2 VISCOSITY (LB/FT-SEC) VS TEMPERATURE (PEG R) 


DATA (MUIVST(I) ,1=1,34) 
TEMPERATURES 
X 160. # 360* * 

X 2X60. , 2520. » 

X 4320. » 46B0* * 

H2 VISCOSITIES 
X 2,54 E-6 * 4,47 fc-6, 

X 14.97 L-6 , 16,54 E-6, 

X 23,45 E-6, 24.66 E-6, 


/ 0 . , 1 , , 15 , , 0 , , 

720, , i080. * 1440. , 
2880. t 324q, , 36QQ, , 
5040 , , 

7.30 E-6, 9,54 E-6, 11.48 E-6, 
16,04 E-6 , 19.47 E-6, 20.84 E-6, 
25,89 E-o / 


1800. , 
3960. » 


13.29 E-6, 
22.16 E-6 r 


02 VISCOSITY (LB/FT-SEC) VS TEMPERATURE (DEG R) 


DATA (MU2VSTU) ,1=1,34)/ 0, , 1. , 15. , 0, , 

TEMPERATURES 


X 180, » 

3t>0 • , 

720. 

, 

1080. , 

1440. , 

1800. , 

X 2160, t 

2520. f 

2880. 

r 

3240. t 

3600. t 

3960. ...j. 

X 4320. , 

4680, , 

5040, 

, 




02 VISCOSITIES 






X 5.15 E-6, 

9.94 E-6, 

17,24 E 

-6, 

22,94 E-6, 

27.87 E-6, 

32.30 E-6, 

X 36,33 E-6, 

40,10 L-6» 

43,74 E' 

-6, 

47.29 E-6, 

50,65 E-6, 

53.90 E-6, 

X 57.03 E-6, 

60,08 L-6, 

63,03 E< 

-b 

/ 



H20 VAPOR VISCOSITY (LB/FT-SEC) 

vs 

temperature 

(DEG R) 


DATA (MU3V5TU) ,1=1,32)/ 0. , 

1. 

, 14. , 0. 

9 


temperatures 







X 360. r 

720. » 

1080, 

, 

1440. , 

1800. , 

2160, , 

X 2520, t 

2680, , 

3240, 

f 

3600. t 

3960. t 

4320. , 

X 4680. t 

5040. » 






H20 VAPOR VISCOSITIES 






X 5,18 E-6, 

9.62 E-6, 

14.44 E- 

-6, 

19.35 E— 6, 

24.10 E-6, 

28.67 E-6, 

X 32.99 E-6, 

37,07 E-o, 

40,92 E' 

-6 , 

44,60 E-6, 

48.U E-6 , 

51.49 E-6, 
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X 54.73 E-6, 57.87 £-6 / 

C 

C HE -VISCOSITY (LB/FT-SEC) VS TEMPERATURE (DEG R) 

C. _ __ 

DATA~CMG&VST<"l)»X=r»34)/ 0. *■ 1. 15, , 6. , 

c temperatures 

X 180, » 360. r 720. , 1080. » 1440. , 1800. , 

X 2160, » 2520. » 2880, , 3240, , 3600, , 3960, r 

X 4320, » 4680, , 5040. , 

.£ HE.. „V_i SCP5.IXIIS 

X 6.71 £-6, 10.54 E-6, 16.50 E-6, 21,43 E-6, 25,80 E-6, 29.79 E-6, 

X 33.50 E-6, 37.01 E“6, 40,34 E-6, 43,52 £-6, 46,58 E-6, 49,54 E-6, 

X 52.40 E-6, 55,17 E-6, 57.87 E-6 / 

C 

L FIVE SPECIFIC HEAT TABLES FOLLOW 

C . 

c 

L H2 SPECIFIC HEAT (BTU/LB-DEG R) VS TEMPERATURE (DEG R) 

C 

DATA (CFIVST(I) ,1=1,34)/ 0, , 1. , 15, , 0. , 

C TEMPERATURES 


X 

180, 

9 

360 * 

9 

720, 

, 

1060. 

, 

1440. , 

1.8. QQ, 

J, 

X 

2160. 

9 

2520, 

9 

2880* 

, 

3240. 

, 

3600. , 

3960. 

9 

X 

. 4320. 

9 

4680, 

9 

5040. 

, 






H2 

SPECIFIC 

HEATS 









X 

£.675 

, 

3.233 

9 

3,460 

, 

3.474 

, 

3.509 , 

3.577 

t 

X 

3,670 

, 

3.772 

9 

3.872 

, 

3,966 

, 

4,050 , 

4.126 

9 

X 

4,19.3 

r 

4,253 

9 

4*306 

/ 







c, 

C 02 SPECIFIC HEATS (BTU/LB-DEG R) VS TEMPERATURE (JEG R) 
C 

DATA (CF2VSTU) ,1=1,34)/ 0. , 1. , 15. , 0, , 

C TEMPERATURES 


X 

180. 

, 360, 

, 

720. 

F 

1060. i 

r 1440. , 

1600. 

! 

X 

2160. 

, 2520, 

, 

2880, 

9 

3240. i 

• 3600. , 

3960. 

9 

X 

4320. 

, 4680. 

, 

5040, 

9 





02 

SPECIFIC 

HEATS 








X 

,217 

, .217 

, 

.224 

9 

,239 . 

. .251 , 

.260 

9 

X 

,266 

, .270 

, 

.274 

9 

.278 i 

- .281 , 

.285 

9 

X 

,288 

, ,291 

t 

,294 

/ 






C 

L H20 VAPOR SPECIFIC HEATS (BTU/LB-uEG R> VS TEMPERATURE (DtG R) 
C 

DATA (CF3VST ( I ) , 1=1,34)/ 0, , 1. , 15. , 0. , 

C TEMPERATURES 



X 

180, 

9 

360 • 

9 

720, 

9 

1080. 

t 

1440. 

f 

1800. 

F 


X 

2160, 

9 

2520. 

9 

2680 , 

9 

3240. 

9 

3600 . 

9 

3960. 

F 


X 

4320. 

9 

4680, 

9 

504C. 

9 







c 

H20 

VAPOR 

SPECIFIC HEATS 










X 

,4406 

9 

,4411 

9 

.453 

9 

.480 

t 

.512 

9 

.546 

9 


X 

.d79 

9 

.609 

9 

,636 

9 

.659 

9 

.678 

9 

.694 

9 

c 

c 

X 

.709 

9 

.721 

9 

.731 

/ 







H2G 

SOLID 

SPECIFIC HEATS 

(BTO/LB-dEG 

H) VS 

TEMPERATURE 

(QhG 

R) 



w 

data (CF3SVrU> ,1=1,16)/ o. , 1. , 6. , 0. , 

c temperatures 

X 200. , 30u. t 350, , 400. , 450. , 491.4 » 

L H20 SOLID (ICE) SPECIFIC HEATS 

X .227 , ,320 , ,367 , ,41b , .462 , .502 / 

C 

C HE SPECIFIC HEAT (BTU/LB-DEG R) VS TEMPERATURE (DEG R) 
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c 

DATA. (CF4VS.Ta).»l=JU16)/ 0^. t JU t 6, t 0, . 

c temperatures 

. .. JL..MCU l__2M£L , — 5MSL. j 

C HE SPECIFIC HEATS 

X 1.24 , 1 ,2ft.. t 1.24 . , .1.24 i 1.24 , X*24.._. .. 

C 

C VAPOR PRESSURE (PSIA) vs .TEMPERATURE (DEG R) ... 

C 


DATA (VP.VST.lI}.fI=li60) / Q.»_ CU*. . 


c 

temperatures 








X 

210. 

t 

281. t 

301. » 

331 • , 

360, 

1 

380, 


X 

401. 

, 

420. 1 

440. , 

460, » 

480, 

9 

492. 


X 

513. 

1 

562. r 

586. , 

622. , 

653. 

9 

688. 


X 

727. 

, 

768. , 

768, , 

818. , 

842, 

9 

877, 


X 

9JQ5. 

t 

927. t 

1005. , 

, 1155. , 


m 

„ 

c 

pressures 









X 

..17 E-il 

t 

,441 E-8 , 

.588 E-7 , 

.147 £-5 , 

,232 E-4 

9 

.00012 


X 

.00056 

t 

.0019 , 

.0062 , 

,0185 , 

.0505 

9 

.0885 


X 

0.2 

( 

1.0 t 

2.0 , 

5.0 1 

10. 

9 

20. 


X 

40. 

, 

75, , 

100. , 

150. , 

200 . 

9 

300. 

t 

X 

. 400. 

1 

500. » 

1000, , 

3000* / 





V 

c 

temperature 

(DEG R) VS VAPOR PRESSURE (PSIA) 





c 

UATA (TVSVP(I) fl=l»48) / 0. , 1. » 22. , 0. , 


C PRESSURES 


X .20 1 

.50 

f 

1.0 

1 2.0 

f 

3.5 

1 

5*0 

-- 4 

X 7.5 , 

10. 

9 

15. 

, 20. 

f 

25. 

9 

30. 

9 

X 40. 

50. 

9 

75. 

, 100. 

9 

150. 

9 

2Q0, 

f 

X 300, 

400. 

9 

500. 

, 1000. 

f 





temperatures 










X 513. 

540, 

9 

562, 

» 586, 

t 

608. 

9 

622. 

9 

X 640. 

653. 

9 

673. 

, 688. 

f 

70 U. 

9 

•710, 

9 

X 727. 

741. 

9 

768. 

, 788. 

9 

818. 

9 

842. 

9 

X 877. 

905. 

9 

927. 

, 1005. \ 

/ 






C 

C DELTA HEAT OF CONDENSATION VS TEMPERATURE 
C 

DATA (DHCVST (I)., 1=1 *24)/ 0. * 1. » 10. , 0. , 

L TEMPERATURES 



X 513, , 

562. 

, 566, , 

622, 

9 

653. 

9 

688, , 


X 710. » 

741. 

, 768. , 

788. 

t 




c 

DELTA HEATS OF 

CONDENSATION 







X 1064. , 

1036. 

, 1022. , 

luOl, 

9 

932, 

» 

960 , , 

c 

c 

X 945. , 

924. 

, 904, , 

889, 

/ 




TAbLE OF THERMAL 

CONDUCTIVITY OF H2 vs 

TEMPERATURE 

: (DEG 

R) 



data (tckd, 

1-1,22) 

/ 0», 1., 9 . , 

0 • , 





V. 



TEMPERATURES 






X 132. , 

312. 

, 492. , 

672. 

t 

852. 

9 

1032. , 


A 1200. t 

1500, 

t 2000, , 





. 


C ‘ THERMAL. CONDUCTIVITIES (H2) 

X .oi4 E-5 # 1,011 E-5, 2.663 E-5, 3.444 E-5, 4,122 E-5, 4.736 E-5, 
X 5.292 E-5, 6.194 E-5, 7,583 E-b/ 


C 

c table of thermal conductivity of 02 vs temperature (deg ri 
data (TC2(I) ,1=1,20) / 0. , l,f 8,» 0., 
c temperatures 

X 132. , 312. , 492* , 582. , 672. » 300. , 

X 1000. , 2000, , 

C THERMAL CONDUCTIVITIES (02) 
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X * 106 E-5, .253 E-5, ,394 E-5, .461 E-5, .522 E-5, .600 E-5, 

X .714 E-5, 1.222 E-5/ 

TABLE OF THERMAL CONDUCTIVITY OF H20 VAPOR VS TEMPERATURE _( PEG R) 

DATA (TC3V (i ) , I=T»ToT/ 0 . » IT, 8., Q.» 

TEMPERATURES 

X' 672. , 852, , 1032, , 1212. » 1392. , 1600. , 

X 1800, » 2000, , 

THERMAL CONDUCTIVITIES (H20 VAPOR) 

X .378 E-5, ,506 E-5, .639 E-5 p ,775 E-5, .911 E-5, 1.061 E-5, 

X 1.208 £-5, 1.356 E-5/ 

TABLE OF THERMAL CONDUCTIVITY OF HELIUM VS TEMPERATURE (DEG R) 

DATA (TC4(I) ,1=1,18) / 0., 1,, 7., C., 

TEMPERATURES 

X 132, , 312. , 492. , 672. , !Q0Q, . #„.A500j 

X 2000, , 

THERMAL CONDUCTIVITIES (HELIUM) 

X .939 £-5, 1.700 £-5, 2.272 E-5, 2,744 E-5, 3.458 E-5, 4.347 E-5, 
X 5.139 E-5/ 


END 
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